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Absorptiometric  titrations  of  a  model  drug  compound, 
3-aminoacridinium  cation,  with  salts  of  calf  thymus  DNA 
were  carried  out  at  pH  5.90  in  a  series  of  solutions  in 
which  the  supporting  electrolyte  was  the  dihydrogen 
phosphate  salt  of  one  of  the  following:  lithium,  sodium, 
potassium,  rubidium,  cesium,  or  tetramethylammonium  ion. 
Titrations  were  also  conducted  in  solutions  of  the  acetate 
salts  of  magnesium,  calcium,  strontium,  and  barium.   The 
ionic  strengths  of  the  solutions  were  varied  between  0.25 
and  0.00063 (M)  for  each  alkali  metal  electrolyte  and  between 
0.036  and  0.00010  for  the  alkaline  earth  electrolytes. 
Apparent  equilibrium  association  constants  for  the  surface 


xiv 


binding  and  for  the  intercalative  binding  of  3-amino- 
acridinium  cation  to  DNA  were  computed  for  each  concen- 
tration of  each  alkali  metal   electrolyte  solution  and 
are  reported.   Apparent  constants  for  the  surface  binding 
of  3-aminoacridinium  in  the  presence  of  alkaline  earths 
are  also  presented.   Values  for  the  intercalative  binding 
mode  in  the  presence  of  group  IIA  metal  cations  could  not 
be  determined. 

The  extent  of  binding  of  the  model  drug  to  DNA  was 
dependent  on  the  nature  and  concentration  of  the  electrolyte. 
Essentially  complete  binding  was  attainable  for  all  concen- 
trations of  all  alkali  metals  while  as  little  as  70%  binding 
was  possible  in  the  presence  of  0.010N  magnesium.   Spectral 
characteristics  of  the  drug-DNA  complex  were  also  dependent 
upon  the  nature  of  the  background  electrolyte  and  are  dis- 
cussed. 

The  reaction  systems  did  not  behave  in  accordance  with 
limiting  Debye-Hiickel  relationships  over  the  range  of  ionic 
strength  investigated.   After  applying  extended  Debye- 
Hiickel  corrections,  the  relationships  between  the  surface- 
binding  equilibrium  association  constants  and  the  ionic 
strength  indicated  that  specific-ion  competition  between 
metal  and  drug  cations,  for  DNA  phosphate  sites,  was  signifi- 
cant.  A  simple  1:1  competitive  model  was  evaluated. 

At  ionic  strengths  less  than  about  0.010,  experimental 
precision  was  poor,  presumably  due  to  alterations  in  the 
conformation  of  DNA  at  low  ionicities,  suggesting  that 


evaluations  of  interactions  between  small  molecules  and 
DNA  be  conducted  in  solutions  of  ionic  strength  in  excess 
of  0.010. 

Thermodynamic  parameters  for  surface  and  intercalative 
binding  of  3-aminoacridinium  cation  and  its  benzolog,  7- 
aminoquinolinium,  were  evaluated  from  equilibrium  association 
constants  at  15.0°C,  25.0°C,  and  35.0°C.   The  standard 
enthalpy  change,  AH°,  for  the  3-aminoacridinium  binding 
was  less  negative  than  AH°  for  the  7-aminoquinolinium  process, 
while  the  standard  entropy  change,  AS°,  for  the  former  was 
more  positive  than  for  the  latter.   Thus,  the  greater  degree 
of  disordering  resulting  from  the  relaxation  of  the  solvent 
around  the  tricyclic  3-aminoacridinium  cation  as  it  binds  to 
DNA  appears  responsible  for  the  more  negative  standard  free 
energy  change  for  the  3-aminoacridinium  binding  process. 
Values  of  AS°  for  the  intercalative  binding  of  the  two  drugs 
suggest  that  significant  ordering  of  the  polymer  helix  occurs 
upon  intercalation. 
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CHAPTER  I 

INTRODUCTION 

The  interaction  of  small  molecules  with  DNA  and 
biopolymers,  in  general,  has  been  the  subject  of  intense 
investigation  during  the  past  several  decades.   Interest 
in  this  area  has  been  due,  in  large  part,  to  the  early 
realization  that  such  interactions  may  result  in  muta- 
genesis (1)  and  carcinogenesis  (2) ,  especially  when  the 
small  molecules  are  cationic  aromatics  such  as  substituted 
acridines.   Numerous  publications  deal  with  various  aspects 
of  the  interactions  of  cationic  aromatics  with  DNA,  some 
of  which  will  be  briefly  reviewed  here.   The  specific 
systems  to  be  considered  in  this  work  include  3-amino- 
acridine(I).    Hence,  a  summary  of  representative  and 
pertinent  research  involving  aminoacridines  is  included. 
Moreover,  we  have  investigated  the  effects  of  metal  ions 
on  the  interactions  of  3-aminoacridine  with  DNA  and, 
therefore,  shall  review  the  past  research  efforts  in  this 
area.   It  is  emphasized  that  this  review  is,  by  no  means, 


Roman  numerals  in  parentheses  after  chemical  names  refer  to 
their  listings  in  Figure  1  of  appendix  I. 


exhaustive,  but  rather,  is  intended  to  provide  perspective 
and  an  introduction  to  the  work  that  is  included  in  this 
report. 

Binding  of  Small  Molecules  to  Biopolymers 

Probably  the  first  acknowledgement  of  a  significant 
interaction  between  a  member  of  the  acridine  chemical 
family  and  a  biological  system  was  reported  in  1913  by 
Ehrlich  and  Benda  (3)  when  they  reported  their  observation 
on  the  antibacterial  activity  of  3 , 6-diamino-10-methyl- 
acridinium  chloride (II)  against  Trypanosonoma  brucei. 
Subsequently,  evidence  accumulated  that  the  antibacterial 
activity  of  substituted  aminoacridines  is  dependent  on 
their  capacity  to  compete  with  hydronium  ion  for  a  vital 
anionic  site  in  a  bacterium  (4-6,  7:  Ch.  2)  and  that 
antibacterial  activity  and  mammalian  toxicity  were  not 
directly  related  (8,9).   The  antibacterial  action  of 
prof lavine (III)  is  reported  to  be  due  to  its  ability  to 
inhibit  DNA-dependent  RNA  synthesis  (10)  which  suggests 
that  DNA  may  contain  the  anionic  site.   However,  it  is 
inadvisable,  at  this  time,  to  make  a  direct  correlation 
since  proflavine  also  binds  strongly  to  a  variety  of  other 
cellular  structures,  any  one  of  which  may  contain  the 
vital  anionic  site  (11) . 

Numerous  amino-substituted  acridines  are  known  to  be 
carcinogenic  while  few  simple  acridines  are  cancer  inducing 
(12)  which  may  be  due  to  the  former's  greater  ability  to 


intercalate  into  DNA  (13) .   It  is  worthy  of  note  that 
benzacridines (IV)  and  dibenzacridines (V) ,  first  exten- 
sively studied  in  1935,  are  nearly  all  highly  carcino- 
genic (2) .   Two  possible  modes  of  action  have  been 
proposed.   The  benzacridine  may  be  metabolically  activated 
to  a  cationic  form  which  then  reacts  as  an  electrophile 
toward  nucleophilic  centers  of  DNA,  RNA,  and  protein  (14) , 
or  the  unmetabolized  species  may  interact  physically  with 
the  DNA  (15) ,  perhaps  by  intercalating  into  the  double 
stranded  helix. 

Mutagenesis,  presumably  arising  from  the  intercalation 
of  acridines  with  nucleic  acids,  was  first  reported  by 
DeMars  in  1953  (  1  )  in  relation  to  his  work  on  proflavine 
and  T2  bacteriophage  (Escherichia  coli) .   The  mutagenic 
activity  of  proflavine  does  not  involve  its  integral 
incorporation  into  the  newly  synthesized  polynucleotide 
chain  but  rather,  results  from  the  acridine's  capacity  to 
cause  deletions  and  insertions  of  nucleotides  into  the 
DNA  polynucleotide  (16,17).   In  fact,  proflavines'  capacity, 
in  this  regard,  was  used  in  1961  to  discover  the  triplet 
nature  of  the  genetic  code. 

In  searching  for  effective  anticancer  and  antitumor 
agents,  advantage  is  taken  of  the  fact  that  compounds 
which  may  interfere  with  cell  replication  at  the  molecular 
level  may  have  medicinal'  utility.   As  a  result,  a  large 
number  of  acridine  derivatives  have  been  investigated,  only 
to  find  that,  while  a  select  few  do  have  antineoplastic 


properties,  their  associated  toxicities  to  the  host  make 
them  less  desirable  than  other  available  compounds.   On 
the  other  hand,  aminoacridines  have  enjoyed  success  as 
antimalarial  agents  and  are  still  in  limited  use  today. 
Investigations  were  begun  in  the  192  0' s  for  a  substitute 
for  quinine,  as  an  antimalarial,  which  led  scientists 
first  to  methylene  blue,  a  phenothiazine  dye,  then  to 
pamaquine,  both  of  which  had  low  effectiveness  against 
acute  stages  of  malaria.   Structural  similarities  in  the 
above  two  compounds,  with  substituted  acridines,  eventually 
led  investigators  to  quinacrine (VI)  which  became  the  drug 
of  choice  throughout  World  War  II.   Later,  the  more  effec- 
tive and  less  toxic  chloroquine (VII)  was  discovered. 
Presently,  considerable  synthetic  research  is  devoted  to 
acridine-derived  antimalarials  to  supplant  chloroquine 
and  to  combat  chloroquine-resistant  P.  falciparum  — 
discovered  in  1961  (18) .   To  date,  the  most  promising 
compounds  of  this  type  have  a  diamine  side  chain  in  the  9- 
position (VIII ) .   Structure-activity  relationships  for  a 
variety  of  quinacrine  analogues  having  side  chain  variants, 
terminal  amine  variants,  and  aromatic  substituent  variants, 
as  well  as  quinine  related  acridine  antimalarials,  are 
tabulated  by  Henry  (19)  . 

In  light  of  the  broad  biological  significance  of  the 
binding  of  aminoacridines  to  DNA  and  other  biopolymers, 
in  vitro  and  in  vivo  studies  aimed  at  understanding  the 
exact  nature  of  the  interactions  have  been  undertaken. 


The  following  paragraphs  will  briefly  review  several  of 
the  more  important  articles  germane  to  our  investigation. 
Those  dealing  with  the  modes  of  binding  will  be  outlined 
first,  followed  by  studies  of  the  effects  of  the  ionic 
medium  (including  possibly  competing  ions) ,  with  a  review 
of  papers  involving  theoretical  discussions,  third. 

An  excellent  comprehensive  article  on  the  interactions 
of  aminoacridines  with  nucleic  acids  and  the  methods  by 
which  they  may  be  studied  was  published  in  1956  by  A.R. 
Peacocke  and  J.N.H.  Skerrett  (20).   Another,  in  1968,  by 
Blake  and  Peacocke  (21)  considered  developments  in  the 
field  after  1956. 

Aminoacridines  have  long  been  known  to  bind  to 
double-stranded  DNA  via  two  distinct  modes  —  termed  type 
I  and  type  II  binding.   Type  I  is  the  stronger  and  was 
first  proposed  by  Lerman  (22) ,  to  be  an  internal  or 
intercalative  process.   Insertion  of  a  molecular  ion  is 
accomplished  by  transient  rotation  of  base  pairs  allowing 
the  flat,  lipophilic,  aromatic  portion  of  the  molecule  to 
become  situated  between  them.   The  charged  portion  of  the 
molecule  was  considered  to  remain  external  to  the  helix. 
H.J.  Li  and  D.M.  Crothers  (23)  corroborated  this  model 
when  they  presented  indirect  evidence  that  the  base  pairs 
remain  perpendicular  to  the  helix  axis  upon  intercalation 
of  the  molecule.   Another  model  for  type  I  binding  was 
advanced  by  N.J.  Pritchard,  A.  Blake,  and  A.R.  Peacocke  (24) 
in  which  the  small  molecule  is  intercalated  between  adjacent 


bases  on  the  same  polynucleotide  chain.   In  both  of  these 
models,  the  cationic  heteroatom  is  associated  with  an 
indeterminate  number  of  anionic  phosphate  moieties  along 
the  DNA  backbone.   The  viscosity  of  heat-denatured  DNA 
is  essentially  unaffected  by  binding  of  acridines  (25) 
whereas  the  viscosity  of  native  double-stranded  DNA  is 
increased.   An  increase  in  viscosity  is  indicative  of  a 
lengthening  of  the  polymeric  chain  such  as  would  result 
from  intercalation  between  base  pairs  or  between  adjacent 
bases.   Since  the  single-stranded  DNA  undergoes  no  vis- 
cosity change,  the  Lerman  model  appears  more  plausible. 

Some  specific  geometrical  parameters  of  substituted 
acridine-DNA  intercalative  complexes  were  established  by 
G.R.  Kelly  and  T.  Kurucsev  (26)  through  the  evaluation  of 
linear  dichroic  spectra  of  stretched  DNA  films.   An  advan- 
tage of  films  over  conventional  methods  involving  aqueous 
solutions  is  that  they  do  not  require  large  ratios  of  drug 
to  DNA  phosphate.   Taking  advantage  of  the  dichroism  of 
two  mutually  perpendicular  polarized  transitions  of  the 
9-aminoacridine  nucleus,  the  workers  were  able  to  assess 
the  angles  of  tilt  and  twist  of  the  plane  of  the  dye 
relative  to  the  polymer  helix.   The  long  axis  of  acri- 
flavine(II)  was  essentially  perpendicular  to  the  helix 
axis  (85  +  5°),  prof lavine (III)  was  tilted  slightly  (77  + 
2  1/2  ),  while  9-aminoacridine (IX)  was  tilted  significantly 
(67  +  3°). 


Type  II,  surface,  or  external  binding  has  come  to 
be  generally  accepted  as  an  association  of  the  small 
molecules  with  the  exterior  surface  of  the  DNA  polyanion. 
This  interaction  is,  consequently,  more  strongly  affected 
by  alterations  of  the  solvent  system  than  is  type  I 
binding.   Li  and  Crothers  (23)  determined,  on  the  basis 
of  temperature- jump  kinetic  data,  that  proflavine  binds 
both  internally  and  externally  to  calf  thymus  DNA  and 
only  about  7%  of  the  total  amount  bound  binds  externally 
in  the  presence  of  0.2M  Na   (sodium  phosphate  -  sodium 
nitrate) .   Decreasing  the  sodium  ion  concentration  to 
0.02M  results  in  about  30%  of  the  total  being  externally 
bound . 

A  careful  study  of  the  intercalative  binding  of  a 
series  of  aminoacridines  was  conducted  by  Drummond  and 
co-workers  (27) .   Binding  decreases  in  affinity  in  the 
order:  acranil (X) >neomonacrin (XI) >atebrin(VI) >9-amino- 
acridine (IX) >prof lavine (III )  9-amino-l ,2,3, 4-tetra- 
hydroacridine (XII) .   This  sequence  shows  that  side  chain, 
and  ring,  substituents  are  not  predominant  in  influencing 
binding  affinity.   The  effect  of  greatly  decreasing  the 
planar  area  of  9-aminoacridine  by  hydrogenation  is 
evident  upon  comparing  its  position  with  that  of  XII.   The 
authors  discuss  various  geometries  of  binding  and  conclude 
it  is  not  necessary  that  there  be  exact  and  complete 
intercalation  for  type  I  (strong)  binding.   Rather,  a 
modified  Lerman  model  in  which  the  charged  heterocyclic 
ring  nitrogen  interacts  with  a  DNA  phosphate  while  the 


rings  only  partially  interact  with  the  bases  is  more  in 
accord  with  observed  behavior. 

It  has  been  shown  that  the  degree  of  intercalation 
may  be  dependent  on  the  base  composition  of  the  DNA. 
Heterogeneity  of  intercalative  sites  was  demonstrated 
by  J.C.  Thomes,  G.  Weill,  and  M.  Daune  (28)  using 
fluorescence  quenching  techniques.   They  demonstrated 
that  between  2%  and  3%  of  the  base  pairs  form  sites  where 
proflavine  is  strongly  bound  while  the  remainder  have 
binding  constants  3  to  4  times  weaker.   Proflavine 
fluorescence  is  relatively  unaffected  in  the  former  case, 
but  is  totally  quenched  in  the  latter.   The  stronger 
sites  correspond  to  adenine-thymine  rich  regions  of  the 
DNA  indicating  proflavine' s  specificity  for  the  base 
pair.   Similarly,  Ramstein  and  Leng  (29)  showed  that  the 
location  of  bound  proflavine,  within  the  DNA,  was  dependent 
upon  the  base  composition.   The  absorptiometric  and  fluori- 

metric  titrations  of  DNA  and  partially  methylated  DNA 

7 
(  18%,  mainly  on  N   of  the  guanine  residues)  indicated 

acridine  was  intercalated  in  both  the  native  and  methylated 

sites.   However,  greater  apparent  equilibrium  constants  for 

the  methylated  DNA  than  those  for  the  native  type  were  ob- 

4  -1 
served  at  two  ionic  strengths:  14.25  x  10  M   and  4  3.5  x 

104M_1  in  0.1M  NaCl  and  40  x  104M_1  and  83.5  x  104M_1  in 

0.01M  NaCl,  respectively.   These  differences  may  be  due  to 

an  increase  in  the  distance  between  base  pairs  arising  from 

the  presence  of  the  methyl  groups  which  may,  in  turn,  result 


in  the  shift  in  the  equilibrium  concentrations.   S. 
Georghiou  (30)  obtained  much  the  same  results  as  Ramstein 
and  Leng  based  on  the  fluorescence  decay  of  proflavine 
when  it  was  bound  to  DNA  at  a  phosphate  to  drug  ratio  of 
420.   Gabbay  and  co-workers  (31)  invoked  the  "Bookmark" 
model  of  intercalation  to  postulate  ten  possible  distinct 
sites  of  binding  of  N-substituted-N-methylphenanthrolinium 
cations (XIII)  to  nucleic  acids.   They  also  investigated 
the  topography  of  some  nucleic  acids  in  solution  using 
PMR  (32) . 

Binding  models  for  the  interactions  of  acridine 
orange (XIV)  and  prof lavine (III)  with  DNA  at  ionic  strengths 
of  0.002M,  0.020M,  and  0.20M  (tetramethylammonium  caco- 
dylate  buffer,  pH  6.5,  22°C)  were  proposed  based  on 
thermodynamic,  spectroscopic, and  hydrodynamic  properties 
of  the  systems  (3  3) .   The  methods  precluded  consideration 
of  surface  binding.   Binding  isotherms  were  generated  from 
equilibrium  association  constants  which  were,  in  turn, 
obtained  from  theoretical  mass  action  relationships. 
Additional  information,  useful  in  developing  binding 
models,  was  afforded  by  Goswami,  Das,  and  Das  Gupta  (3  4) 
who  measured  the  decrease  in  the  static  dielectric  con- 
stant of  a  prof lavine-DNA  complex  while  decreasing  the 
DNA  phosphate/prof lavine  (P/D)  ratio.   The  authors  con- 
firmed that:  the  bound  dye  neutralizes  the  charge  on  the 
DNA  polyanion;  at  high  P/D  values,  intercalation  predomi- 
nates and;  that  when  P/D  becomes  small  ($10) ,  a  cooperative 
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increase  in  the  electrostatic  binding  of  phosphate  groups 
brings  about  a  rapid  decrease  in  the  total  surface  charge 
of  the  polymer.   In  brief,  at  P/D  values  between  10  and 
100  most  of  the  cations  are  intercalated  with  a  concomi- 
tant elongation  of  the  polyanionic  DNA  helix,  whereas 
below  P/D  of  10,  surface  binding  is  predominant.   In  the 
latter  case,  the  spatial  extension  of  the  helix  is  de- 
creased due  to  charge  neutralization  with  no  accompanying 
expansion  due  to  intercalation.   These  results  indicate 
that  alterations  in  the  dielectric  constant  of  a  medium 
containing  such  a  system  would  manifest  dissimilar  effects 
on  the  two  equilibrium  association  constants. 

The  specific  effects  of  cations,  other  than  the  dye 
species,  on  the  reaction  between  DNA  and  a  dye  is  an 
integral  part  of  this  investigation.   Therefore,  a  brief 
consideration  of  articles  dealing  with  the  binding  of 
metal  cations  to  DNA  and  related  topics  is  in  order.   A 
tabulation  of  sites  and  thermodynamic  quantities  associated 
with  proton  and  metal  ion  interactions  with  DNA,  RNA,  and 
their  constituent  bases,  nucleosides,  and  nucleotides  has 
been  published  by  Izatt  et  al . ,  (35)..  H.  Sigel  and  D.B. 
McCormick  (36)  provide  an  introduction  to  the  fundamental 
theories  concerning  metal-ligand  interactions  in  biological 
systems,  including  discussions  of  the  Irving-Williams 
series,  competition  of  metals  and  protons  for  sites,  the 
formation  of  ternary  complexes,  and  effects  of  solvent 
polarity  and  ionic  strength.   The  Pearson  theory  (37)  of 
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hard  and  soft  acids  and  bases  is  applied  to  metals  and 
polymers  in  relation  to  the  biological  roles  of  Na  ,  K  , 
Ca   ,  and  Mg    and  to  the  carcinogenicity  of  heavy  metals 
(38)  . 

The  relative  effectiveness  of  ions  in  altering  the 
conformation  of  macromolecules  via  specific  ion  inter- 
actions is  termed  the  Hofmeister,  or  lyotropic,  effect 
(39).   Generalizations  based  on  observations  of  simple, 
ideal  polymers  may  be  confidently  extrapolated  to  more 
complex  polymers  since  it  has  been  shown  that  the 
Hofmeister  effect  of  a  given  ion  is  remarkably  independent 
of  the  nature  of  the  macromolecule.   Moreover,  while 
changes  in  the  entropies  of  polymer  segments  are  affected 
by  both  anions  and  cations,  as  are  changes  in  the  partial 
molal  volume,  the  partial  molal  internal  energy  of  a  segment 
is  dependent   only  on  cations.   Considering  the  above, 
M.J.  Hey,  J.M.  Clough,  and  D.J.  Taylor  (39)  reported  that 
the  ion  binding  strength  with  protein  is  (in  decreasing 
order  of  strength)  NH4+<K+<Na+  and  Na+,  Li+<H20<Ca++. 
This  series  may  be  useful  in  separating  solvent  effects 
from  ion  binding  effects  where  the  above  cations  are 
present  along  with  3-aminoacridinium  cation  and  DNA. 

Stabilizing  effects  of  various  ions  on  four  biopoly- 
mers,  including  DNA,  were  evaluated  by  Hippel  and  Schleich 
(40) .   The  effects  were  discussed  in  terms  of  free  energy 
of  transfer  for  model  compounds,  activity  coefficient 
variations,  and  correlations  between  ion  effects  on 
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macromolecules  and  on  water  structure.   The  relative 

effectiveness  of  the  ions  studied,  in  increasing  the 

stability  of  the  native  configuration  of  DNA,  paralleled 

the  Hofmeister  series:  (CH3)4N+>K  >Na  >Li   and  Cl~,  Br~> 

CH-.COO- >C10  ~>CNS~.   It  was  found  for  a  series  of 
3        4 

tetraalkylammonium  salts,  the  greater  the  length  of 
the  alkyl  chains,  the  lower  the  stabilization  of  the 
native  form  of  DNA. 

Hen  and  Strauss  (41)  employed  equilibrium  dialysis, 
dilatometry,  and  viscometry  to  assess  the  counterion 
binding  of  a  series  of  cations  with  poly (vinylsulfonate) . 

The  Hofmeister  series  was  maintained  (in  order  of  de- 

+   +    +   +    +        + 
creasing  binding  affinity):  Ag  >K  >Na  >H  ,Li  >(CH.J.N 

and  Ba   >Mg   .   All  of  the  cations  except  H   and  Li 
are  capable  of  cross-linking  the  polyanion  through 
simultaneous  site  binding  at  two  sulfonate  groups. 

The  specific  interactions  of  univalent  and  di- 
valent cations  with  calf  thymus  DNA  can  be  correlated 
with  their  respective  unhydrated  ionic  radii  (42,43). 
The  order  of  most  weakly  bound  to  most  strongly  bound 
is:  (CH  )  N+<K+<Na+<Li+<Ca++<Mg++<<Mn++.   Comparing 
binding  parameters  with  DNA  with  polyphosphate  as  polyanion, 
the  authors  showed  that  the  sites  on  the  biopolymer  were 
the  phosphates  rather  than  bases  or  base  pairs.   These 
results  were  in  agreement  with  those  of  other  workers 
(44,45) . 
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Conductometric  titrations  of  DNA  with  divalent 
cations  were  first  performed  by  J.  Shack,  R.J.  Jenkins, 
and  J.M.  Thompsett  (4  6)  who  found  a  sharp  end  point  cor- 
responding to  0.8  equilvalents  of  Mg    and  Ca    per  mole 
of  polymer  phosphate.   Similar  titration  behavior  was  seen 
for  poly  A  and  poly  U  by  Felsenfeld  and  Huang  (47)  who 
developed  an  improved  binding  model  which  showed  that  the 
apparent  stoichiometry  of  0.8  was  in  error  and  that  the 

true  equivalence  point  corresponded  to  1.0  eq  of  Mg    or 

++ 

Ca   per  mole  of  DNA  phosphate.   The  authors  also  disproved 

the  previously  held  belief  (48)  that  dications  are  more 
strongly  bound  to  purines  than  pyrimidines  in  a  poly- 
nucleotide. 

It  has  been  reported  (49)  that  there  is  binding  to 
base  pairs  of  the  DNA  helix  as  well  as  to  phosphates.   The 
investigators  reported  that  the  effects  of  adding  Mg   to 
heat-denatured  calf  thymus  DNA  was  to  produce  an  absorption 
spectrum  similar  to  that  of  denatured  DNA  alone,  except 
for  reduced  absorbance.   However,  when  Mg   was  added  to 
native  DNA,  which  was  then  denatured,  a  spectral  shift  to 
longer  wavelength  occurred.   The  authors  concluded,  on 
these  grounds,  that  the  metal  was  binding  only  to  phos- 
phates in  the  former  case  while  in  the  latter,  there  was 
interaction  with  the  aromatic  nitrogens  of  the  bases.   The 
bonding  between  Mg    and  the  nitrogens  is  presumably  not 
able  to  occur  unless  the  hydrogen  bonds  between  the  bases 
are  broken  in  the  presence  of  the  metal. 
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Mathematical  relationships  incorporating  changes  in 
activities  of  ions  in  solution  and  deviations  in  activity 
coefficients  were  developed  by  Lyons  and  Kotin  (45)  to 
qualitatively  distinguish  between  specific  (site)  binding 
and  non-specific  (diffuse  electrostatic)  binding  of 
metals  to  polyanions.   The  rule  of  additivity  of  the 
activities  of  countercations  of  solutions  of  polyelectro- 
lytes  and  simple  salts  was  re-examined.   It  was  found 
that,  under  certain  conditions,  there  is  significant 
exchange  between  the  counterions  associated  with  the 
polyanion  and  those  of  the  bulk  solution.   A  decrease  in 
the  activity  coefficient  of  sodium  ion,  upon  dilution  of 
the  system,  was  observed,  indicating  increased  binding  at 
lower  concentrations.   In  the  case  of  magnesium  the  results 
were  not  as  simple  but  suggested  that  at  high  magnesium 
ion  concentrations,  specific  site  binding  was  predominant 
while  at  low  concentrations,  non-specific  binding  obtained. 
The  differences  between  Na   and  Mg   may  be  due,  in  part, 
to  the  ability  of  the  latter  to  bridge  between  two  adjacent 
anionic  sites  on  the  polymer  whereas  the  monovalent  cation 

cannot.   Note  that  the  distance  between  adjacent  phosphates 

o 
along  the  DNA  backbone  is  approximately  7  A.   C.Sander  and 

and  O.P.  Ts'o  (50)  believe  that  the  binding  of  Mg++  to  DNA, 

RNA,  poly  A,  poly  A-U,  poly  I,  poly  I-C,  and  denatured  DNA 

can  be  rationalized  using  a  linear  Scatchard  relationship 

(vide  infra) .   Their  specific  site  model  was  applicable  for 

up  to  about  70%  saturation. 


15 


The  degree  of  binding  of  acridine  orange  to  native 
and  denatured  DNA  was  determined  by  equilibrium  dialysis 
in  0.1M  and  0.001M  NaCl  at  20°C  (51).   Three  successive 
stages  of  binding  were  observed:  the  first  corresponding 
to  intercalation  of  monomers;  the  second,  formation  of 
bound  dimers;  and  the  third,  external  binding  of  aggre-  • 
gates.   The  appearance  of  the  three  modes  of  binding 
compared  to  two  for  most  aminoacridines  is  attributed 
to  acridine  orange's  proclivity  to  dimerize.   The  dialysis 
data  was  augmented  with  absorption  and  fluorescence 
spectra,  electric  dichroism,  electric  birefringence,  and 
circular  dichroic  spectra. 

Scruggs  and  Ross  (42)  found  that  the  intrinsic 
viscosity  of  three  types  of  DNA  (salmon  sperm,  calf  thymus, 
and  T4  phage)  decreases  in  the  presence  of  univalent  cations 
as  the  ionic  strength  increases,  attaining  separate  limiting 
values  for  each.   Prior  to  their  work,  the  apparent  insen- 
sitivity  of  the  viscosity  of  DNA  to  variations  in  ionic 
strength  was  considered  anomalous  because  polyelectrolyte 
viscosities  are  usually  highly  dependent  on  salt  concentra- 
tion.  The  workers  asserted  that  previous  invar iance  was 
due  to  the  presence  of  trace  amounts  of  basic  proteins, 
polypeptides,  and  polyamines  and  that  DNA  was,  in  fact, 
similar  to  all  other  polyelectrolytes . 

S.  Bram  (52)  investigated  the  secondary  structure  of 
DNA  in  solutions  of  varying  concentrations  of  Li  ,  Na  , 
Rb  ,  and  Cs   using  X-ray  diffraction  spectroscopy. 
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Differences  between  scattering  patterns  by  sodium  DNA 
in  solution  and  those  expected  from  theoretical  calcu- 
lations based  on  an  unperturbed  B  form  of  DNA  were  due 
to  structural  changes  in  the  B  form  arising  from  the 
polymer's  interaction  with  the  species  in  the  medium. 
This  is  in  concert  with  others' (53-55)  conclusions  that 
the  number  of  base  pairs  per  turn  of  the  DNA  helix  is 
a  function  of  the  ionic  environment.   For  cesium  DNA 

the  extrapolated  lengths  for  10  nucleotide  pairs  are 

o  o 

43.5  A  at  zero  ionic  strength  and  34.6  A  at  finite  CsCl 

o 
concentration,  whereas  in  NaCl ,  the  values  are  3  3.5  A 

o 
and  32.0  A,  respectively.   (A  novel  DNA  conformation 

model  which  differs  markedly  from  the  classical  Watson- 
Crick  B  form  is  proposed  by  Stig  Er lander  (56)  who  argues 
that  the  differences  between  the  values  given  above  cannot 
be  rationalized  using  the  Watson-Crick  model.)   A  theo- 
retical treatment  of  cooperative  binding  was  advanced  by 
G.  Schwarz  (57)  in  which  nearest  neighbors  along  a  linear 
polymer  interact  with  one  another.   Two  types  of  intrinsic 
internal  binding  are  assumed:  (1)  that  of  an  isolated 
ligand  (nucleation) ;  (2)  that  of  a  ligand  to  a  site 
immediately  adjacent  to  one  already  occupied.   The  model 
also  includes  considerations  for  dimer  formation  and 
competitive  (non-cooperative)  binding  at  adjacent  sites 
by  unrelated  ligands.   Matrix  methods  of  calculating 
species  concentrations  under  equilibrium  conditions  are 
presented. 
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An  understanding  of  the  effects  of  ionic  strength 
on  the  binding  of  small  molecules  to  DNA  requires  at 
least  some  knowledge  of  the  charge  distribution  on  the 
surface  of  the  polyanion.   Electrophoresis  and  membrane 
equilibrium  experiments  allow  the  determination  of  the 
apparent  fractional  charge  per  DNA  phosphorous,  i,  which 
is  an  index  of  net  charge  (58) .   The  value  of  i  increases 
from  0.250  to  0.34  on  going  from  a  solution  containing 
0.005m  NaCl  to  one  having  0.05m  NaCl,  showing  that  the 
net  charge  of  the  polyanion  increases  with  increasing 
ionic  strength.   Ross  (58)  considered  the  above  information 
in  light  of  Gorin's  model  of  a  long  rigid  cylinder  (59) 
and  applied  it  to  DNA,  which  is  intermediate  between  a 
rigid  cylinder  and  a  random-coil  polyelectrolyte .   Ross1 
data  substantiates  Gorin's  assumption  that  a  rigid  cylinder 
(and,  therefore,  DNA)  obeys  Debye-Hiickel  approximations. 

Manning  (60)  postulated  that,  as  the  ionic  strength 
approaches  zero,  the  effective  charge  of  a  polyelectrolyte 
is  maintained  at  a  critical  constant  by  condensation  of 
counterions  along  its  surface  as  long  as  the  formal  charge 
does  not  exceed  a  critical  value.   His  predictions  were 
borne  out  by  independent  investigators  (33)  who  found 
that  at  low  ionic  strength  (0.002M)  the  electrophoretic 
mobility  of  DNA  was  unaffected  by  bound  acridine  orange. 

Lerman's  model  for  intercalation  of  small  molecules 
into  DNA  satisfies  many  experimental  observations  but 
leaves  some  uncertainties  unresolved  arising  from 
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consideration  of  the  energetics  involved  (61) .   One  of 
the  most  perplexing  problems  is  the  strong  dependence 
of  the  intercalative  mode  of  binding  on  the  ionic  strength 
(albeit  not  as  dependent  as  surface  binding) .   Also,  the 
existence  of  a  maximum  degree  of  association  as  a  function 
of  ionic  strength  requires  explanation.   Attempts  have 
been  made  to  resolve  these  problems  by  invoking  "stacking 
energies"  rationalized  in  terms  of  short  range,  ionic 
strength  independent  interactions  between  DNA  base  pairs 
and  the  aromatic  portions  of  the  acridines.   These  ratio- 
nalizations, however,  are  incompatible  with  the  fact  that 
aromatic  hydrocarbons  intercalatively  bind  only  very  weakly. 
M.  Gilbert  and  P.  Claverie  (61)  point  out  that  the  solution 
energies  of  the  intercalating  compounds  parallel  their 
binding  affinities  and,  therefore,  may  be  at  least  as 
important  as  the  stacking  energies.   Moreover,  the  limiting 
of  intercalation  processes  long  before  all  available  sites 
have  been  occupied  mandates  a  limitation  mechanism  inde- 
pendent of  the  heterogeneity  of  sites.   The  critical  effect 
of  transferring  positive  charges  from  solvent  "cavities," 
having  dimensions  on  the  order  of  the  small  molecules', 
into  large  cavities  interior  to  the  lipophilic  DNA  must 
be  accounted  for  in  a  model  of  the  dye-DNA  complex  —  as 
it  is  in  their  model. 

T.  Herskovits  (62)  studied  the  relationship  between 
electrolyte  concentration  and  solvent  denaturation  transi- 
tions in  various  media.   Both  increasing  chain  length  and 
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increasing  hydrocarbon  content  increase  the  effectiveness 
of  the  denaturant.   It  was  concluded  that  these  observations 
demonstrate  the  importance  of  hydrophobic  forces  in  main- 
taining DNA's  aqueous  configuration.   In  most  cases,  the 
denaturation  was  reversible  with  added  electrolyte. 

It  is  Bradley's  opinion  (63)  that  by  using  a  combina- 
tion of  statistical  and  quantum  mechanical  methods,  it 
should  be  possible  to  relate  all  absorption,  emission, 
equilibrium  constant,  and  optical  rotation  observations 
in  a  unified  way.   With  this  aspiration  in  mind,  he 
presented  a  statistical  model  of  dye-polymer  binding  as 
a  step  toward  that  goal.   He  assumed  that  specific  polymer 
sites  were  in  a  linear  array  and  dye  molecules  were  bound 
to  them  by  simple  electrostatic  forces,  with  possible 
enhancement  by  neighboring  sites.   Furthermore,  the  sites 
were  sufficiently  close  to  allow  bound  dye-bound  dye 
interactions  with  concomitant  changes  in  their  absorption, 
fluorescence,  and  phosphorescence  spectra.   The  distribution 
of  N  consecutive  sites  occupied  by  dye  molecules,  in  their 
varying  states  of  aggregation,  can  be  computed  theoretically 
in  terms  of  either  the  ground  state  free  energy  of  inter- 
action of  a  pair  of  neighboring  dye  molecules,  AF,  or  in 
terms  of  a  stacking  coefficient,  K,  related  to  AF, 

K  =  exp(-AF/kT)    .  (1-1) 

The  AF  and  K  terms  may  be  equated  to  an  experimentally 
obtainable  ratio  of  P  sites  per  D  dye  species, 
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P/D  =  (1-F1/2)  1(1+(K-1)F)  (1-2) 


Values  of  K  range  from  1.25  for  DNA  to  826  for  poly- 
phosphate (K  =  6.2  for  denatured  DNA) .   There  is 
evidence  that  the  value  of  K  is  directly  related  to 
the  conformation  of  the  polymer.   The  fact  that  K  is 
equal  to  1.25  and  not  to  1.00  for  DNA  —  as  it  would  be 
for  a  completely  random  site  occupancy  —  indicates  some 
non-randomness  which  could  arise  from  dye-dye  interaction 
along  the  polymer.   A  simplified  quantum  mechanical 
treatment  is  presented  to  allow  computation  of  the  magni- 
tudes of  the  relationships  between  the  dye  molecules  of 
filled  units  with  empty  sites. 

Kinetics  of  intercalation  of  aminoacridines  into  DNA 
were  first  studied  in  depth  by  Li  and  Crothers  (23)  in  1969. 
Employing  data  from  temperature- jump  relaxation  methods, 
they  proposed  a  relatively  simple  two-step  process  for 
intercalative  binding.   It  is  worth  noting  that  static 
experiments  will  not  allow  a  separation  of  these  two  steps 
of  the  intercalative  mode.   The  data  fit  the  mechanism 

k12  k23 

P  +  DNA    ■   -*-z>»    (P)  ZJ  v    (P)  .        (1-3) 

s  .   "     out   <*r-r in 

K21  32 

where  P  is  proflavine  and  (P)   ,  and  (P) .   are  surface 
r  out        in 

bound  and  intercalatively  bound  dye,  respectively.   The 

overall  insertion  reaction  occurs  in  the  millisecond  time 

range  and  is  first  order  from  the  (P)   ,  state.   The  two- 
out 

step  hypothesis  is  substantiated  by  comparing  the  large 
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relative  thermodynamic  and  spectral  changes  of  the 
exterior  binding  with  the  small  changes  for  the  interior- 
bound  species  that  occur  upon  glucosylation  of  the  DNA 
(vide  infra) .   Two  limiting  mechanisms  are  advanced  and 
alternate  mechanistic  pathways  discussed. 

Passero,  Gabbay,  and  others  (64)  used  a  simplistic 
model  which  included  competitive  effects  of  metal  ions 
in  solution  to  arrive  at  apparent  association  constants 
of  two  reporter  molecules  (Xllla,  XHIb)  to  DNA.   Electro- 
static potentials  and  interactions  between  near  neighbors 
were  ignored,  as  were  any  forms  of  surface  binding  by  the 
reporter  molecules.   The  overall  binding  process  was 
expressed  as  two  separate  equilibria 

R  +  P  — ►  RP  (1-4) 

and 

M  +  P  — >  MP  (1-5) 

where  R,  M,  and  P  are  free  reporter,  metal,  and  DNA 
phosphate,  respectively,  and  RP  and  MP  the  corresponding 
bound  species.   An  apparent  binding  constant,  K'_  may  be 
defined  relative  to  the  metal  ion  concentration 

K'R  =  [RP]  /  [R]  [P]   .  (1-6) 

Finally,  assuming  the  binding  of  R  or  M  to  the  polyanion 

results  in  only  one  phosphate  being  affected  per  R  or  P, 

the  apparent  association  constant  may  be  expressed  in  terms 

of  association  constants  for  reporter  binding,  K  ,  and 

R 
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metal   binding,    K 

K'R  =   KR  /    (1   +   KM    [M+])  .  (1-7) 

The  authors  maintain  that  their  assumptions  are  valid 
by  virtue  of  the  fact  that  they  obtain  reasonably 
constant  values  of  K   and  K   under  a  variety  of  conditions. 
However,  if  the  reporter  molecules  are  intercalatively 
bound,  as  the  authors  assert  (and  they  are  probably 
correct)  then  it  is  not  valid  to  compare  the  reporter 
and  metal  binding  because  two  entirely  different  types 
of  sites  are  being  occupied.   The  highly  hydrophilic, 
symmetrically  charged  metal  cations  are  almost  certainly 
surface  bound.   Thus,  the  constancy  of  the  results  are 
probably  fortuitous  or,  at  best,  important  deviations 
are  being  masked  by  experimental  error. 

Relaxation  kinetics  of  ethidium  bromide (XV)  binding 
to  DNA  revealed  three  types  of  bound  species  at  equilib- 
rium (65) .   Aside  from  the  conventional  surface-bound  and 
intercalated  species  ascribed  to  type  II  and  type  I  binding, 
a  third  species  was  found  which  involved  an  ethidium  bromide 
transiently  attached  to  sites  of  two  separate  DNA  strands. 
The  bimolecular  rate  constant  for  transferring  an  ethidium 
bromide  from  one  DNA  phosphate  to  another  is  3  to  6  times 
larger  than  the  rate  of  intercalation  of  free  molecules. 
The  more  rapid  bimolecular  mechanism  is  dependent  on  the 
small  molecule's  ability  to  bind  to  two  sites  simultaneously. 
In  the  case  of  ethidium  bromide,  the  tricyclic  region  of  the 
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molecule  provides  one  site  while  the  phenyl  moiety  is 

believed  to  provide  the  second,  less  favorable,  site. 

Divalent-cation-specific  electrodes  were  used  to 

study  the  binding  of  Mg++  to  DNA  (50)  in  a  solvent 

.  .         -3 
containing  5  x  10   M  pH  7  phosphate  buffer.   Results 

were  treated  in  terms  of  site-binding  to  phosphates  in 

a  multiple  equilibrium  process,  employing  a  Scatchard 

treatment,  which  showed  linearity  over  a  range  of  30% 

to  7  0%  saturation  of  sites.   Approximations  to  the 

Scatchard  method  yielded  an  apparent  intrinsic  binding 

constant  of   6  x  10   (AG°  =5.2  Kcal/mole) .   The  use  of 

cation-specific  electrodes  mandated  DNA  phosphate  concen- 

—  3 
trations  in  excess  of  5  x  10   M,  causing  aggregation  of 

DNA  and  precipitation  of  MgDNA.   Thus,  the  values  of  the 
constants  may  include  significant  error. 

Three  methods  for  evaluating  equilibrium  association 
constants  for  the  binding  of  small  molecules  to  DNA  will 
be  outlined.   They  are,  the  classical  and  most  commonly 
employed  Scatchard  method,  a  modified  Scatchard  model, 
and  a  mass-action  approach  which  will  be  used  for  this 
investigation.   Scatchard 's  original  treatment  was  designed 
for  binding  of  small  molecules  to  proteins.   In  Scatchard' s 
words: 

If  the  various  groups  on  a  protein  molecule 
act  independently,  we  can  apply  the  Law  of 
Mass  Action  as  though  each  group  were  on  a 
separate  molecule  and  the  strength  of  binding 
can  be  expressed  as  the  constant  for  each 
group.   Often  a  single  constant  will  express 
the  behavior  of  several  groups.  (66:660) 

If  the  groups  do  act  independently  of  one  another,  changes 
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in  free  energy  for  the  reaction  of  the  protein  with 
small  molecules  is  made  up  of  the  statistical  entropy 
terms  for  v  plus  a  term  proportional  to  v.   when  the 
total  number  of  sites  (or  groups  of  sites)  is  large,  the 
calculations  become  tedious  —  if  the  total  number  is 
not  known  the  calculations  cannot  be  performed.   Such 
a  situation  does  arise  if  one  attempts  to  apply  the 
equations  without  modification  to  nucleic  acid  binding. 
Since  DNA  polyanions  commonly  have  molecular  weights 
in  excess  of  one  million  and  may  contain  as  many  as  10 
phosphates  (67)  Peacocke  and  Skerrett  (20)  redefined  the 
parameters  of  the  Scatchard  equations  to  allow  their 
application  to  nucleic  acids.   Their  assignments  and  an 
abbreviated  derivation  of  the  Scatchard  equations  are  as 
follows. 

Let  R  represent  the  amount  of  small  molecules,  D, 
bound  per  mole  of  nucleic  acid  phosphorus  (generally  a 
fraction) ,  and  C  represent  the  molar  concentration  of 
free  D.   Consider  a  number,  n.,  of  P  classes  of  binding 
sites  for  D  on  the  DNA  polyanion.   If  each  class  has 
associated  with  it  an  intrinsic  binding  constant  K.,  then 
we  may  write  for  any  general  case 


J=P      n.K.C 

R  =  .A       rHb  (1~8) 

j=i  j 


This  may  be  simplified  when  all  sites  are  identical 

R/C  =  Kn  -  Kr  (1-9) 


* 
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If  there  are  two  distinct  types  of  binding  sites,  I 
and  II,  then 

n  K  C        n   K   C 

R  °  1  f  KlC     +    1  +  k"C  d"10) 

Figures  2a  and  2b  are  representative  plots  of  R/C  vs  C 
for  one  and  two  types  of  sites,  respectively.   Curvature 
may  result  in  the  case  of  only  one  type  of  site  from 
variations  in  the  electrostatic  free  energy  term  which 
is  dependent  on  R.   Such  non-linearity  can  become  critical 
when  the  intercalating  species  bear  opposite  charges  but 
may  be  partially  obviated  by  replacing  K  with  the  term 

K'j  exp  (AG°  /(RT)  (1-11) 

where  K'   is  independent  of  R,  and  AG°   is  an  electro- 
static free  energy  term  dependent  only  on  the  total  effec- 
tive potential.   The  case  in  which  there  are  two  types  of 
binding  sites  is  of  primary  interest  in  this  work.   The 
stronger,  intercalative  binding  corresponds  to  type  I  of 
figure  2b  whereas  surface  binding  is  represented  by 
region  II.   With  two  types  of  sites,  it  is  sometimes 
possible  to  discern  distinct  linear  portions  in  the  plots 
of  R/C  vs  R  as  suggested  in  the  idealized  figure.   Only  if 
the  magnitude  of  one  of  the  constants  is  about  10  to  100 
times  greater  than  the  other,  can  reliable  values  of  the 
individual  association  constants  be  obtained  by  extrapola- 
tion of  the  data  in  each  of  the  linear  regions.   Even  then, 
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potentially  erroneous  assumptions  regarding  changes  in 
the  electrostatic  free  energy  must  be  considered.   More- 
over, as  previously  stated,  the  proposed  model  does  not 
allow  for  any  form  of  interaction  between  the  two  types 
of  binding  sites.   This  is  not  generally  serious  in 
protein  binding  where  the  sites  are  often  noninteracting 
but  can  be  a  problem  in  the  case  of  the  nucleic  acids. 
Factors  such  as  cooperative  binding  and  near-neighbor 
exclusion  will  introduce  sigmoidal  characteristics  into 
the  plots  of  R/C  vs  R  (21) .   Another  shortcoming  of  the 
Scatchard  method  is  that  extrapolation  of  the  legs  of 
the  hyperbolic  plot  to  obtain  slopes  and  intercepts 
often  places  an  inordinately  high  degree  of  reliance 
on  the  data  at  the  extrema  of  titrations.   It  is  at  the 
extrema  that  the  data,  regardless  of  how  it  was  obtained, 
will  be  the  least  precise. 

For  aminoacridines,  intercalative  binding  pre- 
dominates up  to  a  value  of  R£0.20  corresponding  to  4  or 
more  phosphates  per  bound  cation.   Surface  binding,  on 
the  other  hand,  is  predominant  for  values  of  R  between 
approximately  0.5  and  1  for  compounds  which  do  not  dimerize 
and  up  to  2  or  more  for  those  that  are  known  to  aggregate 
as  free  species  at  moderate  concentrations. 

Numerous  investigators  have  suggested  improvements 
of  the  basic  Scatchard  method  as  modified  by  Peacocke  and 
Skerrett.   Included  are  corrections  for  electrostatic 
effects  (27,31,61,63),  competition  of  metallic  ions  in 
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solution  (45,46,67),  near-neighbor  effects  (23,50,57, 
63,65),  small  molecule  aggregation  (51,57),  and  configura- 
tional  effects  due  to  solvent  and  ionic  factors  (23,25,28, 
39,44,52,62,68).   One  of  these  improved  models,  proposed 
by  Armstrong,  Kurucsev,  and  Strauss  (3  3)  incorporates 
some  novel  features  including:  (1)  the  total  concentration 
of  intercalation  sites  available  is  fixed,  a  priori,  to 
include  every  slot  between  successive  DNA  base  pairs;  (2) 
an  intercalated  dye  inhibits  intercalation  at  its  two 
immediately  adjacent  slots;  and  (3)  inclusion  of  a  term 
to  allow  for  an  intercalated  monomer  forming  a  dimer  with 
a  free  monomer  resulting  in  a  spectroscopically  distinct 
species.   The  authors  found  (3)  became  important  when  the 
degree  of  binding,  B,  exceeded  about  0.2  mole  of  dye  per 
mole  of  DNA  phosphate  (dyes  were  acridine  orange  and  pro- 
flavine) .   Their  expressions  for  the  association  constants, 
based  upon  theoretical  mass  action  relationships  are,  for 
intercalation 


2(3n-3?)  (1-23,) 
K  =  ± f -±-  (1-12) 

CM(1-401> 

and  for  surface  binding 

e2 

K9    C  (8  -B  )  (1-13) 

1         ^MVP1  *V 

where  B,  =  the  number  of  intercalated  molecules  per  DNA 
phosphate. 
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32  =  the  number  of  externally  bound  diraers  per 

DNA  phosphate. 
C  =  the  molar  concentration  of  free  dye  monomer. 
K,  =  the  apparent  binding  constant  for  the  monomer 

species  intercalating. 
K„  =  the  apparent  binding  constant  for  surface 
binding. 
These  relationships  imply  that  the  electrostatic  potential 
of  the  DNA  is  unaffected  by  binding.   Values  of  3,  and  g„ 
may  be  spectroscopically  evaluated  knowing  the  molar 
absorptivities  of  the  bound  monomer,  e  ,  and  the  bound 
dimer,  e  , 

£T  =  (£M3M  +  eDBD)/3  d-14) 

where  3M  =  \~^2'    3D  =  2fV  and  £T  is  the  totai  "m°lar 
absorptivity"  of  all  bound  species.   The  model  was  tested 
by  determining  the  binding  of  acridine  orange  and  pro- 
flavine to  DNA  (PH  6.5,  22°C)  at  ionic  strengths  of  0.002m, 
0.02  0m,  and  0.2  00m.   The  above  method  is  an  improvement 
over  the  Scatchard  procedure  but  still  demands  a  priori 
selection  of  an  exact  mode  of  reaction.   That  is,  though 
the  overall  stoichiometries  are  experimentally  evaluated, 
the  thermodynamic  manner  in  which  the  free  species  react 
is  not  defined.   For  example,  in  the  intercalative  case, 
where  4  phosphates  are  involved  for  each  3-aminoacridinium 
binding  site,  are  the  4  phosphates  to  be  considered  as  4 
separate  entities,  as  2  separately  interacting  moieties, 
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each  containing  2  phosphates,  or  as  one  entity  containing 
all  4  phosphates?   In  each  case,  the  overall  stoichiometry 
is  4  to  1,  but  the  equilibrium  expression  of  the  first 
includes  a  fourth  power  term  in  DNA  phosphate  concentration, 
the  second  would  have  a  squared  phosphate  concentration 
term,  while  the  third  would  be  linear  in  phosphate. 

To  resolve  this  problem,  Capomacchia  (69)  has  proposed 
a  mass-action  approach  based  on  best-fit  analyses  of 
spectroscopically  derived  concentrations  of  free  and  bound 
species.   Assuming  that  the  binding  process  obeys  the  Law 
of  Mass  Action  and  total  concentrations  of  all  species  are 
sufficiently  low  that  aggregation  of  like  species  is  in- 
significant, we  may  write 

BH  +  qSu       N   BH(Su)  (1-15) 

SI 

where  BH  is  the  free  cationic  form  of  the  small  molecule, 
q  is  the  number  of  DNA  binding  sites  that  complex  one 
small  molecule  and  Su  denotes  the  unbound  DNA  binding 
sites.   Each  site,  Su,  acts  as  one  distinct  entity  in 
the  equilibrium  expression  and  may  contain  any  number 
(1,2,3,.  .  .  m)  of  phosphates.   Note  that  throughout  this 
text  binding  sites  will  be  discussed  in  terms  of  DNA 
phosphates,  simply  for  convenience,  because  the  total 
concentration  of  DNA  phosphate  is  easily  obtained  ex- 
perimentally.  This  terminology  does  not  imply,  for 
instance,  that  phosphates  are  the  binding  sites  for 
intercalatively  bound  drug.   The  association  constant, 
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on   a  molar  basis,  for  the  process  is 

K  -  [BHP]/[BH] [Su]q  (1-16) 

where  [BHP]  is  the  molar  concentration  of  bound  species, 
[BH]  is  the  concentration  of  free  small  molecule,  and 
[Su]  is  the  concentration  of  the  uncomplexed  DNA  sites. 
Unbound  site  concentration  may  be  determined  from  the 
total  site  concentration,  S.f  and  bound  site  concentration, 

S 

c 

.   Su  =  St  "  Sc  f1"17) 

Now,  since 

Sc/q  -  [BHP]  (1-18) 

and  st  =  [Pt/m]  (1-19) 

where  P   is  total  DNA  phosphate  concentration,  we  may 
write 

Su  =  [pt/m]~  qfBHp]    •  (1-20) 

Substituting  equation  (1-20)  into  (1-16)  yields 

k. urn •  (1.21) 

[BH]  ([Pt/m]-  q[BHP])q 

Taking  the  log  of  both  sides  of  (1-21)  and  rearranging, 
we  obtain 

log  (Tfif^)  =  log  K  +  qlo<?  (tpt/m]-  q[BHP])q      (1-22) 
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Molar  concentrations  of  free  and  bound  small  molecule 
species  may  be  determined  from  any  of  a  variety  of 
experimental  techniques.   The  overall  stoichiometry  of 
the  process  may  be  obtained  from  a  Scatchard  treatment  or 
a  Job's  plot.   Knowing  the  overall  stoichiometry,  initial 
values  of  m  and  q  may  be  selected  which  satisfy  the  simple 
relationship  mq  =  n,  where  n  is  the  total  number  of  DNA 
phosphates  per  bound  small  molecule.   A  rapid  determination 
of  the  most  probable  combinations  of  values  may  be  made  by 
plotting  log  ([Pt/m]-  q[BHP])  vs.  log  (  [BHP]/'[BH]  )  .   That 
pair  which  yields  the  most  consistent  straight  line  having 
a  slope  of  q  may  then  be  used  to  compute  values  of  K  for  a 
series  of  ratios  of  [BHP]  to  [BH] .   Again,  the  pair  which 
provides  the  most  consistent  calculated  equilibrium  associa- 
tion constants  for  all  data  points,  corresponding  to  either 
the  surface  or  intercalative  mode,  are  taken  as  the  exponents 
in  the  binding  expression.   For  interactions  between  3- 
aminoacridinium  cation  and  calf  thymus  DNA,  best-fit  values 
of  m  and  q  are  3  and  1, respectively  for  surface  binding  and 
2  and  2  for  the  intercalative  mode  (69)  .   To  summarize,  the 
method  of  Capomacchia  allows  for  the  determination  of 
equilibrium  association  constants  based  on  mechanistic 
considerations.   Furthermore,  by  computing  values  of  K  from 
data  in  the  midregion  of  a  titration,  undue  reliance  on  the 
least  trustworthy  data  at  the  extrema  of  a  titration  is 
avoided. 

Many  workers  have  been  satisfied  with  reporting  associ- 
ation constants  calculated  from  molar  quantities  rather  than 
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from  the  more  correct  activities  of  the  various  species. 
Doing  so  implies  two  things:  activities  of  the  various 
species  in  the  expression  are  equal  to  molarities,  and 
the  molarities  are  equal  to  molalities.   It  is  unlikely 
that  approximating  molalities  with  molarities  introduces 
significant  error  throughout  the  usual  concentration 
ranges  employed.   But  activities,  on  the  other  hand, 
may  differ  markedly  from  molarities.   For  this  reason, 
evaluation  of  association  constants  in  terms  of  activities 
is  recommended.   Molality  is  related  to  activity  by 


a.  =  a.m.  (1-23) 

1     11 

where  a.  is  the  activity,  a.  is  the  activity  coefficient, 
l  u         i 

and  m.  is  the  molality  of  species  i.   Note  that  the 
activity  equals  the  molality  only  when  a  =  1.   For  a 
process  in  solution 

A+  +  nB"  ==^   AB~(n+l)  (1-24) 

^       n 

the  equilibrium  constant  is  rigourously  defined  as 

aAB"(n+1)             aAB  MAB 
K  =  -_" =   _ 5 5 .         (1_25) 

(AA+)(V}         (V*A+*(VV> 

It  is  extremely  difficult  to  determine  individual  activities 
since  their  is  no  way  of  separating  effects  due  to  positive 
ions  from  those  due  to  the  accompanying  negative  ions.   We 
may,  however,  define  a  mean  activity,  a—,  which  is  the 
geometric  mean  of  the  individual  activities 
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A-  =     (A+-A")1/2  (1-26) 


Similarly, 

a-  -  (a+-a~)1/2   .  (1-27) 

For  ideal  dilute  solutions,  the  mean  ionic  activity 
coefficient  will  be  equal  to  1.   If  deviations  from 
ideality  are  caused  entirely  by  electrical  interactions 
of  point  charges,  it  may  be  shown  that,  for  aqueous 
solutions  at  25°C  (70) 

log  a-  =  -.509/Z+Z_/l1//2  (1-28) 

where  Z   and  Z_  are  the  charges  of  interacting  species 
and  I  is  the  total  ionic  strength  of  the  system 

1=1/2  Z      .m.Z2   .  (1-29) 

ill 

Equation  (1-28)  is  an  expression  of  the  Debye-Hiickel 
limiting  law  for  activity  coefficients  and  is  applicable 
only  to  dilute  solutions  containing  point  charges.   As  the 
ionic  strength  becomes  large  (e.g.,  >0.01m)  and/or  the 
ions  become  less  ideal  (become  larger,  more  polarizable, 
etc.)  the  relationship  becomes  increasingly  approximate. 
Under  these  conditions,  an  extended  Debye-Hiickel  equation 
may  be  invoked  which  introduces  a  correction  for  the  finite 
sizes  of  ions: 

-0.509/Z  Z  7l1//2 

log  a-  =  s-75 

1  +  Bdl  '  (1-30) 
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B  is  a  constant  for  a  given  solvent  and  includes  such 
terms  as  solvent  viscosity,  dielectric,  and  temperature 
while  d  is  an  average  effective  diameter  of  the  ions.   A 
factor  not  considered  in  equation  (1-30)  is  the  effect  of 
ions  on  solvent  molecules  and  solvent  structure.   This 
may  become  especially  important  for  solutions  containing 
large  polyanions  such  as  DNA. 

Lyons  and  Kotin  (44,45)  report  a  method  of  assessing 
values  of  a   and  a   independently  of  one  another.   They 
also  determined  a—  for  Na   and  Mg   when  binding  to  DNA 
and  three  synthetic  polyanions.   Their  somewhat  tedious 
procedure  involving  equilibrium  dialysis  allowed  them  to 
distinguish  site-specific  and  non-specific  forms  of  surface 
binding.   The  activity  coefficient  of  Na  was  shown  to 
decrease  with  increasing  dilution  of  metal  ion  in  the 
presence  of  DNA  suggesting  an  increased  degree  of  binding. 

Thermodynamics  of  Binding  Reactions 

A  second  topic  of  this  work  is  concerned  with  the 
thermodynamic  parameters  associated  with  the  binding  of 
small  molecules  to  DNA.   Hence,  a  brief  review  of  several 
pertinent  articles  is  presented  here. 

M.  Gilbert  and  P.  Claverie  (61)  measured  total  energies 
for  an  intercalation  process  and  attributed  them  to  three 
types  of  electrostatic  interactions.   They  were   (1)  the 
attractions  between  the  bound  cations  and  the  DNA  phosphates, 
(2)  the  repulsive  forces  between  cations,  and  (3)  solvent 
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response  to  the  electric  field  generated  by  the  charged 
species.   At  high  ionic  strength  the  phosphates  are 
relatively  neutralized  causing  (1)  to  become  insignificant 
and,  thereby,  causing  a  decrease  in  both  the  strength  and 
extent  of  binding.   Some  limiting  of  the  binding  process 
occurs  even  at  very  low  ionicities  because  cation-cation 
repulsions  are  still  important  with  respect  to  intercalated 
species.   The  third  type  of  interaction  serves  to  reduce 
all  forms  of  electrostatic  energies  and  will  be  most  ef- 
fective for  those  species  most  intimately  in  contact  with 
the  solvent.   The  authors  have  shown  that  phosphates  and 
complexed  acridines  are  more  solvated  than  free  acridines 
so  that,  as  the  ionic  strength  is  raised,  their  stabilizing 
interactions  are  more  reduced  than  are  acridine-acridine 
interactions.   In  that  the  forces  arising  from  the  latter 
are  destabilizing,  the  result  of  the  countervailing  trends 
is  a  net  destabilization  of  the  complex  with  increasing 
ionic  strength.   The  model  used  here  can  also  be  applied 
to  other  cationic  small  molecules,  including  protons. 
Similar  results  were  reported  for  the  binding  of  proflavine 
to  DNA  via  f luorimetry  (71) .   It  was  found  that  the  electro- 
static contribution  to  the  total  free  energy  decreased  with 
increasing  ionic  strength  due  to  screening  of  the  potential 
by  an  ionic  atmosphere. 

Thermodynamic  parameters  accompanying  the  intercalation 
of  calf  thymus  DNA  were  evaluated  by  Chambron,  Daune,  and  Sadron 
(72)  using  equilibrium  dialysis.   Temperatures  ranged  from 
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0°C  to  70°C  and  ionic  strengths  of  0.01M,  0.1M,  and  1.0M 
were  maintained  using  acetate  buffer  at  pH  5.9  in  the 
presence  of  EDTA.   The  number  of  sites  available  for 
binding  decreased  upon  increasing  either  the  temperature 
or  the  ionic  strength.   The  enthalpy  of  reaction  decreased 
with  decreasing  temperature;  AH°=  -19  Kcal/mole  at  70°C 
and  -4.5  Kcal/mole  at  0°C  (I  =  0.01).   This  was  ascribed 
to  a  thermally  reversible  change  of  state  of  DNA  which 
occurs  at  approximately  4  0°C.   Total  standard  free 
energy  was  considered  to  be  a  sum  of  three  terms.   The 
first  of  these,  AG  ,  characterizes  the  dipole-induced 
dipole  interaction  between  base  pairs  and  cations  while 
the  second,  AG  ,  is  the  electrostatic  free  energy.   The 
third,  AG°,  is  the  solvation  free  energy  corresponding 
to  desolvation  of  reactants  and  solvation  of  the  bound 
complex.   At  high  ionic  strengths,  AG°  and  AG°  are  very 
small,  hence,  the  total  free  energy  of  6.4  Kcal/mole  is 
due  to  AG°.   The  authors'  value  of  AS°=  11  Kcal/mole 
degree  is  in  agreement  with  theoretical  predictions. 

Thermodynamic  energies  of  reaction  between  proflavine 
and  two  forms  of  DNA  determined  from  temperature- jump 
relaxation  methods  are  summarized  below  (23) . 


E 
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„+   • 

jj  is  the  activation  energy  for  transformation  from  state  i 


37 


to  j  (see  equation  (1-13)  and  AH?,  and  AS?,  are  the  thermo- 
dynamic enthalpy  and  entropy  changes  for  the  same  reaction. 
Measurements  were  done  over  the  temperature  range  10  C  to 
25°C,  0.20M  Na+  and  pH  6 . 9 .   Of  particular  note  here  is 
that  the  enthalpy  change  for  outside  binding  is  much  more 
strongly  affected  by  glucosylation  of  the  DNA  than  it  is 
for  the  intercalated  form  (T2  DNA  is  glucosylated) .   For 
example,  AH°   changes  from  -9.8  to  -3.4  Kcal/mole  dye  on 
going  from  calf  thymus  to  T2  DNA  but  for  the  overall  inter- 
calation process  (AH°   +  AH°   =  AH°3)  tne  change  is  only 
from  -7.8  to  -7.3  Kcal/mole  dye.   It  is  unlikely  that  small 
differences  in  base  composition  between  the  two  forms  of 
DNA  can  account  for  the  large  differences  in  the  thermo- 
dynamic values  of  the  binding  reaction. 

Entropy  and  enthalpy  changes  during  reaction  of 
ethidium  bromide  with  calf  thymus  DNA  were  computed  from 
temperature -jump  relaxation  data  by  Bresloff  and  Crothers 
(65)  over  the  temperature  range  15°C  to  32°C  in  1.0M  Na  . 
Values  for  the  formation  of  a  surface  bound  species  were: 
AH°  =1.1  Kcal/mole  and  AS°(19°C)  =5.6  Kcal/mole,  while 
for  one  of  a  pair  of  intercalatively  bound  species:  AH  = 
-7.8  Kcal/mole  and  TAS° (19°C)  =  -2.3  Kcal/mole.   Values 
for  a  second  intercalated  form  could  not  be  obtained  due 
to  the  rapidity  of  its  binding  process. 

Bradley  (63)  determined  the  ground  state  free  energy 
of  intercalation  of  a  pair  of  neighboring  acridine  orange 
molecules  bound  to  DNA  to  be  -0.065  Kcal/mole  dye.   By 
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comparing  this  to  the  spectral  shift  accompanying  aggrega- 
tion (504  nm  to  464  nm)  corresponding  to  4.8  Kcal/mole  dye 
separation  in  energy  between  the  ground  and  excited  states, 
he  concluded  that  the  spectral  shift  results  from  an 
increase  in  energy  of  the  excited  state,  rather  than  a 
decrease  in  the  ground  state. 

Enthalpy  and  entropy  changes,  AH  and  AS  respectively, 
of  reactions  may  be  determined  from  the  relationship  be- 
tween the  equilibrium  association  constants  for  the  processes 
at  various  temperatures  and  the  temperature  as  set  forth 
by  the  van't  Hoff  relationship. 

-AH°    1 
log  K  =  „  3Q3R  m  +  Constant  (1-31) 

This  important  equation  will  be  derived  below  (73) . 

For  ideal  systems,  any  change  in  the  total  free 
energy,  G,  can  be  said  to  be  due  to  changes  in  either 
temperature  or  pressure,  since  the  state  function,  G, 
can  be  defined  in  terms  of  pressure,  volume,  and  temperature. 

dG  =  (If)   dP  +  (f§)   dT  (1-32) 

T  P 

We  define 

G  =  H  -  TS  (1-33) 

and 

H  =  E  +  PV  (1-34) 

where  H  is  the  enthalpy,  S  is  the  entropy,  and  E  is  the 


39 


internal  energy  for  the  system  under  consideration.   Recall 
that,  the  enthalpy  changes  are  commonly  termed  "heats  of 
reaction"  for  processes,  entropy  may  be  used  as  a  measure 
of  the  degree  of  order  or  disorder,  and  that  the  internal 
energy  is  unavailable  for  work  of  any  form.   The  change  in 
the  free  energy,  dG,  may  now  be  expressed  as 

dG  =  dE  +  PdV  +  VdP  -  TdS  -  SdT   .  (1-35) 

For  a  reversible  process  in  which  no  work  other  than 
expansion  occurs, 

PdV  =  dw  (1-36) 

and 

dE  -  dq  +  dw  =  0   .  (1-37) 

The  terms  q  and  w  refer  to  the  heat  absorbed,  and  work 
done,  by  the  system  respectively.   Equation  (1-37)  is  a 
direct  consequence  of  the  First  Law  of  Thermodynamics  and 
states,  in  effect,  that  the  total  change  in  the  internal 
energy  of  a  system  is  the  sum  (or  difference)  of  the  amount 
of  heat  absorbed,  dq,  by  the  system  and  the  amount  of  work 
done,  dw,  by  the  system.   Notice  that  heat  and  work  are  not 
functions  of  state  (which  are  conventionally  assigned  upper 
case  symbols) .   Equation  (1-35)  now  becomes 

dG  =  VdP  -  SdT  (1-38) 
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which  upon  equating  with  equation  (1-32)  yields 

(IS)   -  "S  d-39) 

dT   p 

and 

<f§)   "  V   '  (1_40) 

Equation  (1-39)  suggests  an  experimentally  tractable 
relationship  between  free  energy  (via  K)  and  temperature, 
when  pressure  is  held  constant.   Moreover,  we  may  extend 
the  relationship  to  include  the  changes  in  free  energy  as 
a  function  of  temperature  for  a  process  involving  reactants 
going  to  products.   We  may  now  equate  (1-38)  with  equation 
(1-32) ,  rearrange,  and  solve  for  the  partial  differentials 
to  obtain 

(9AG/9T)p  =  -AS  (1-41) 

where  AG  =  (G     ,   ,   -  G      .   .  )  and   S  =  (S     , 

products     reactants  products 

S  reactants^  *   The  entroPy  term  of  equation  (1-41)  may  be 
eliminated  by  recalling  that  AS  =  (-AH  +  AG)/T  and  substi- 
tuting 

(9G/9T)p  =  (-AH  +  AG)/T   .  (1-42) 

Assume,  henceforth,  that  the  pressure  will  be  held  constant. 
Rearranging  (1-42)  gives 
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dAG/dT  -  AG/T  =  -AH/T   .  (1-43) 


Equation  (1-43)  may  be  expressed  as 

T  |^  (AG/T)  -  -AH/T  (1-44) 

since,  upon  differentiation,  it  would  become 


_  T (dAG/dT) -AG    dAG/dT    AG  n  ,n 

T  ~ =  -  .  (l-4b) 

T  T      T 


Equation  (1-45)  may  be  written  in  terms  of  standard  states, 
and,  in  that 

AG°  =  -RT  In  K  (1-46) 

it  can  be  stated  as 

d  In  K/dT  =  AH°/RT2  (1-47) 

which  is  equivalent  to 

d  log  K     -AH°   ,  _    .   .  ,.  ... 

d(l/T)   =  2T303R  +  Constant  (1-48) 

Assuming  that  AH°  is  temperature  independent,  equation  (1-48) 
can  be  integrated  to  yield  equation  (1-3.1)  .   The  above 
relationship  may  be  applied  to  reactions  in  solution  whose 
equilibrium  association  constants  are  based  on  the  activities 
of  the  species  involved  [see  equations  (1-23)  through  (1-30)] 
Plots  of  log  K  vs.  1/T  should  yield  straight  lines  having 
slopes  equal  to  -AH°/2.303R.   If  the  relationship  between 
log  K  and  1/T  is  not  linear  then  the  enthalpy  of  binding  is 
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not  constant  with  respect  to  temperature.   For  systems 
containing  DNA,  enthalpies  of  binding  determined  from  van ' t 
Hoff  plots  are  often  not  constant  with  temperature,  but 
decrease  with  increased  temperature,  presumably  due  to 
changes  in  DNA's  conformation  (72).   As  a  result,  it  should 
not  be  surprising  that  binding  curves  obtained  spectro- 
photometrically  and  by  dialysis  equilibrium  can  agree  at 
one  temperature  and  disagree  at  another  (21) . 

Rationale  for  the  Selection  of  the  Model  System 
and  Experimental  Design 

As  this  study  is  concerned  with  the  effects  of  com- 
peting metal  ions  and  of  alterations  of  the  ionic  strength 
on  the  binding  of  small  molecules  to  biopolymers,  the  choice 
of  DNA  as  the  model  polymer  is  evident.   The  selection  of 
3-aminoacridine  as  the  small  molecule  is  based  on  several 
factors,  two  of  which  are  especially  important.   The  first 
criterion  that  the  small  molecule  must  meet  is  that  it  bind 
to  DNA.   To  do  so,  it  is  necessary  that  the  species  exist 
as  a  cation  since  neutral  compounds  do  not  bind  to  DNA  — 
or  if  they  do,  the  binding  is  too  weak  for  our  purposes. 
(It  should  be  evident  that  electrostatic  repulsive  forces 
prevent  anions  from  interacting  with  DNA.)   Furthermore, 
the  compound  must  be  positively  charged  in  solutions  whose 
pH  is  between  4  and  11  so  that  denaturation  of  the  polymer 
is  minimal  (74).   Hence,  if  the  compound  is  acidic  its  pK 
should  be  at  least  5.   The  pK   of  3-aminoacridine  is  8.04 

a. 

(7:Ch  4)  meaning  it  is  essentially  totally  protonated  in  a 
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solution  of  pH  6.   Moreover,  3-aminoacridine  allows  con- 
venient evaluation  of  both  of  the  principal  modes  of  binding 
of  small  molecules  to  DNA,  under  a  given  set  of  conditions, 
in  a  single  titration.   That  is,  the  range  of  the  ratio  of 
total  drug  to  total  DNA  concentration  over  which  the  two 
modes  of  binding  can  be  experimentally  evaluated  is  about 
0.1  to  10.   Obviously,  it  is  mandatory  that  each  type  of 
binding  occur  in  separate  regions  of  the  titration. 

The  second  major  criterion  that  the  model  compound 
must  satisfy  is  that  it  allow  an  expedient  means  of  measur- 
ing the  extent  of  binding  to  the  polyanion.   Of  the  various 
methods  of  determining  equilibrium  constants  one  of  the  best 
is  absorptiometric  spectroscopy.   The  3-aminoacridine  cation 
is  particularly  amenable  to  this  method  for  the  following 
reasons.   (1)  It  absorbs  strongly  in  the  visible  region  of 
the  spectrum  as  the  free  cation  having  molar  absorptivities 
of  1.413  x  10  M~1cm~1  at  365  nm  and  1.259  x  104M_1cm_1  at 
454  nm  (75) .   Such  magnitudes  allow  good  sensitivity  at  low 
concentrations  of  the  compound.   The  molar  absorptivities 
of  the  bound  drug  in  this  region  are  8.494_  x  10  M~  cm   at 
368  nm  and  9.071  x  10  M~  cm-1  at  463  nm.   The  large  dif- 
ference between  the  absorbances  of  the  band  maxima  of  the 
free  and  bound  drug  provides  good  sensitivity  in  following 
a  titration.   The  longest  wavelength  absorption  maximum 
for  DNA  is  at  approximately  260  nm  meaning  that  there  is 
little  spectroscopic  interference  from  the  polymer.   (2) 
Both  the  free  and  bound  species  have  two  distinct  absorption 
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envelopes  in  the  visible  region  corresponding  to   L,  and 
L  bands  [Piatt's  nomenclature  (76);  a  and  para  bands, 

3. 

respectively,  in  the  Clar  nomenclature  (77)  ] .  These  arise 
from  transitions  from  the   A  ground  state  to  the   L, 
(-365  nm)  ,  and  to  the   L.   (=460  nm)  excited  states.   The 

a 

transition  moment  of  the  former  is  parallel  to  the  long 
axis  of  the  acridine  ring  system  while  the  moment  of  the 
latter  is  orthogonal  to  the   L,  and  parallel  to  the  short 
axis  of  the  molecule.   By  observing  shifts  in  the  two 
maxima,  along  with  any  changes  in  vibronic  character, 
during  the  course  of  a  titration,  one  may  be  able  to 
deduce  binding  geometries  (78)  .   The  3-aminoacridine 
molecule  also  fluoresces  in  its  free  cationic  form  and, 
less  intensely,  as  the  bound  species,  thereby  providing 
another  spectroscopic  method  with  which  to  study  its 
interaction  with  DNA. 

The  more  thoroughly  studied  aminoacridines,  proflavine 
and  acridine  orange,  meet  the  criteria  mentioned  above  in 

regard  to  3-aminoacridine.   Both  of  them,  however,  begin  to 

-5 
aggregate  at  concentrations  as  low  as  10   M,  thereby 

introducing  undesireable  complications.   Also,  the  mono- 
substituted  3-aminoacridine,  being  a  simpler  molecule  than 
either  proflavine  or  acridine  orange,  is  preferable  as  a 
model  compound.   It  is  hoped  that  information  gained  using 
the  simple  compound  may  be  applied  to  more  complex  systems. 

Many  references  to  the  effects  of  ionic  strength  on  the 
binding  of  small  molecules  to  biopolymers  can  be  found  in 
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the  literature.   Unfortunately,  the  majority  are  ancillary 
to  other  studies  and  are  not  sufficiently  complete  in  them- 
selves to  provide  much  specific  information.   Too,  even  for 
a  given  molecule  binding  to  a  given  polymer,  the  wide 
variety  of  experimental  variables  under  which  the  separate 
investigations  were  done  precludes  any  valid  comparisons  of 
results.   The  present  work  was  undertaken  as  a  result  of 
this  lack  of  a  comprehensive  investigation.   Our  pharmaceu- 
tical interest  led  us  specifically  to  DNA  as  the  polyanion, 
though  other,  synthetic  polymers  may  have  provided  a  less 
complex  system. 

Some  of  the  immediate  questions  which  come  to  mind  and 
which  may  be  answered  by  an  investigation  of  the  effects  of 
changing  the  ionic  constitution  of  a  drug-DNA  solution  are: 
are  the  effects  predictable  on  the  basis  of  electrostatic 
theory?   That  is,  may  the  effects  be  explained  in  terms  of 
the  simple  Debye-Huckel  limiting  law  at  low  ionic  strengths 
and  by  extended  Debye-Huckel  equations  at  higher  ionicities? 
If  the  above  questions  cannot  be  answered  affirmatively, 
might  deviations  at  low  ionic  strengths  be  due  to  the 
polymeric  nature  of  the  anion?   Are  deviations  at  high 
ionic  strengths  due  to  specific  ion  interactions  and,  if  so, 
can  they  be  rationalized  on  the  grounds  of  simple  competition 
for  binding  sites  or  must  a  less  simplified,  mixed  complex 
model  be  introduced?   Other  questions  are:  will  altering  the 
ionic  medium  affect  surface  and  intercalative  binding  in  the 
same  manner?   If  not,  how  can  the  differences  be  reconciled? 
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How  will  changing  the  ionic  radii,  ionic  charge,  and 
charge  density  of  the  metallic  ions  affect  the  equilibria? 
The  last  question  is  particularly  germane  in  in  vivo 
biological  binding  studies  because  of  the  wide  variations 

_i_     _i_      +  -+■  4*4- 

in  concentrations  of  Na  ,  K  ,  Mg   ,  and  Ca    in  a  living 
organism  (79) . 

To  try  to  answer  these  questions,  equilibrium  associa- 
tion constants  for  the  surface  and  intercalative  binding 
of  3-aminoacridinium  cation  to  DNA  in  a  variety  of  ionic 

media  were  determined.   Background  electrolytes  included 

-4 
concentrations  between  0.15M  and  6.3  x  10   M  of  the  phosphate 

salts  of  Li  ,  Na  ,  K  ,  Rb  ,  Cs  ,  and  (CHO.N  .   Also  con- 
sidered were  the  acetates  of  the  alkaline  earths:  Mg   ,  Ca   , 
Sr   ,  and  Ba    at  concentrations  ranging  from  0.025N  to  1  x 
10~5N.   The  effect  of  H„PO^  as  opposed  to  CH3COO~  as  the 
counteranion  of  the  supporting  electrolyte  was  evaluated. 
The  pH  for  all  of  the  titrations  was  5.90  and  the  tempera- 
ture was  25.0°C.   A  minimum  of  two  titrations  were  done  for 
each  set  of  conditions. 

Thermodynamic  parameters  were  evaluated  for  systems 
containing  3-aminoacridine  and  DNA  with  a  background 
electrolyte  of  0.010M  KH^PO.  at  pH  5.90.   Similar  investiga- 
tions were  conducted  using  7-aminoacridinium  cation (XVI)  in 
place  of  3-aminoacridinium.   By  studying  the  binding  of  the 
linear  benzalog  of  3-aminoacridine  and  comparing  the  thermo- 
dynamic parameters  of  the  two,  the  relative  importance  of 
enthalpy  and  entropy  contributions  to  overall  free  energies 
may  be  assessed. 


CHAPTER  II 


EXPERIMENTAL 


Materials  and  Procedures 


All  chemicals  were,  unless  otherwise  stated,  either 
reagent  or  analytical  grade.   Water  was  multiply  distilled 
in  a  tin-lined  still  (Barnstead,  Sybron  Corp.).   It  was 
redistilled  until  its  specific  conductance  was  less  than 
1.0  megaohm  cm-1  at  25.0°C.   Stock  solutions  of  supporting 
electrolytes  were  prepared  within  three  weeks  of  use  and 
were  refrigerated  during  storage  to  minimize  biological 
growth.   Approximately  100  to  200  ml  of  a  0.15M  solution 
of  each  electrolyte  were  prepared  as  outlined  below  and 
titrated  with  0.72M  H3P04 ,  1M  CH3C02H,  or  1M  metal  hydroxide 
solution  to  pH  5.9,  in  all  cases.   Concentrations  of  the 
hydroxide  solutions  of  the  various  cations  were  determined 
by  titration  with  standardized  1.000N  sulfuric  acid. 

A  stock  solution  of  tetramethylammonium  phosphate  was 
made  by  titrating  13.16  ml  of  1.14M  tetramethylammonium 
hydroxide  (Mallinckrodt  Organic  Reagent,  10%  in  water, 
lot:  15)  in  70  ml  of  water  with  0.72M  phosphoric  acid 
(Fisher,  lot:  755459)  and  adjusting  the  final  volume  to 
100  ml.   Lithium  phosphate  stock  solution  was  prepared  by 
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titration  of  0.476M  lithium  hydroxide  with  0.72M  phosphoric 
acid.   The  hydroxide  solution  was  made  by  reacting  stoi- 
chiometric amounts  of  anhydrous  lithium  sulfate  (PCR,  Inc., 
Gainesville,  Florida,  lot:  10465)  and  barium  hydroxide 
heptahydrate  (Mallinckrodt,  lot:  RTZ) .   Sodium  and  potas- 
sium phosphate  stock  solutions  were  prepared  using  sodium 
phosphate  monohydrate  (Mallinckrodt,  lot:  ABK)  and  mono- 
basic potassium  phosphate  (Fisher  Certified,  Reagent,  lot: 
20175),  respectively.   The  acid  level  of  each  solution  was 
adjusted  with  0.72M  phosphoric  acid  to  pH  5.9.   The  stock 
solution  of  rubidium  phosphate  was  prepared  by  titrating 
an  aliquot  of  0.4  3M  rubidium  hydroxide  (Pfaltz  and  Bauer, 
no  lot  number)  with  0.72M  phosphoric  acid.   Similarly, 
cesium  phosphate  was  prepared  from  0.201M  cesium  hydroxide 
(Pfaltz  and  Bauer,  no  lot  number) .   Stock  solutions  of 
potassium  acetate  and  tetramethylammonium  acetate  were 
prepared  by  titrating  standardized  aliquots  of  the  respec- 
tive hydroxides  with  1M  acetic  acid  (Baker  Chemical,  lot: 
41397) .   Magnesium  acetate  solutions  were  made  from 
magnesium  acetate  tetrahydrate  (Mallinckrodt,  control: 
NSB) .   Solid  calcium  acetate  was  made  by  dissolving  calcium 
carbonate  (Mallinckrodt,  no  lot  number)  in  excess  acetic 
acid.   The  calcium  acetate  was  recrystallized  twice  from 
aqueous  ethanol,  washed  repeatedly  with  95%  ethanol,  and 
dried  for  twenty-four  hours  at  110°C  to  yield  the  mono- 
hydrate  of  the  salt.   The  sample  was  found  to  contain  97.3% 
Ca(C2H302) .H20,  by  dissolving  a  known  amount  of  the  acetate 
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in  water,  precipitating  the  calcium  as  calcium  carbonate 
and  titrating  the  carbonate  with  hydrochloric  acid. 
Strontium  acetate  was  prepared  from  its  carbonate  (Matheson, 
Coleman  &  Bell,  lot:  19)  and  assayed  in  a  manner  similar 
to  that  for  calcium  acetate.   It  was  determined  to  contain 
95.0%  Sr(C2H_0„).   A  recrystallized  sample  of  barium 
acetate  was  dried  for  twenty-four  hours  at  170  C  and 
assayed  by  forming  barium  carbonate.   The  solid  carbonate 
was  reacted  with  excess  l.OON  sulfuric  acid,  carbonic  acid 
removed  by  gentle  heating,  and  the  excess  sulfuric  acid 
determined  by  titration  with  standard  sodium  hydroxide 
solution.   The  assay  of  the  barium  acetate  was  99.6% 
Ba(C2H302) . 

Solutions  of  the  supporting  electrolytes  to  be  used 
as  solvents  were  made  from  the  stock  solutions.   Aliquots 
of  the  stock  solutions  were  diluted  to  within  90%  of  their 
calculated  final  volumes  and  their  acidities  adjusted  to 
pH  5.90  +  0.05,  if  necessary,  using  the  appropriate  metal 
hydroxide,  phosphoric  acid,  or  acetic  acid.   They  were 
then  diluted  to  their  exact  desired  final  volumes  and  their 
ionic  strengths  computed,  taking  into  account  any  additional 
ions  added  to  adjust  the  acid  levels. 

A  commercially  available  sodium  salt  of  DNA  was 
dialyzed  against  a  solution  of  the  desired  metal  phosphate 
or  acetate  to  obtain  the  DNA  salts  having  the  same  cation  as 
that  of  the  supporting  electrolyte.   Solutions  containing 
approximately  0.01  mole  of  DNA  phosphate  per  liter  were  made 
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by  placing  3.6  mg  of  DNA,  sodium  salt  (Calbiochem  calf 
thymus  DNA,  A  grade,  lot:  900007:  8.04%  P;  12.21%  N; 

moisture,  14%;  E^fiOnm  =  192)  per  milliliter  of  solvent 
into  twenty  milliliter  volumetric  flasks.   Between  5  and 
20  ml  of  solvent  were  then  added.   The  solvent  was  a  0.01M 
solution  (pH  5.90)  of  the  metal  phosphate  or  acetate.   The 
flask  and  its  contents  were  placed  on  a  vertical  rotating 
mixer  (Scientific  Products  model  150-V)  operating  at  6  rpm 
for  about  one  hour,  after  which  the  sample  was  placed  in 
an  ultrasonic  mixer  for  two  to  three  minutes.   Rotation 
and  sonication  were  repeated  for  a  total  of  at  least  five 
hours  or  well  beyond  the  time  at  which  the  moderately 
viscous  sample  appeared  homogeneous.   The  sodium  DNA  solu- 
tion was  then  transferred  to  dialysis  tubing  (Union  Carbide 
Films  Packaging  Division,  7 . 3  mm  dia.).   The  tubing  and  its 
contents  were  immersed  in  a  volume  of  0.01M  solution  of  the 
desired  salt  which  was  at  least  ten  times  greater  than  the 
volume  of  DNA  solution  being  dialyzed.   The  external  solu- 
tion was  agitated  for  two  hours  and  then  replaced.   A 
total  of  six  volumes  of  external  solution  were  used  for 
each  dialysis  procedure.   The  dialyzed  DNA  was  again  mixed 
for  one  hour  on  the  rotating  mixer  and  submitted  to  ultra- 

sonication  for  two  to  three  minutes.   Then,  the  sample  was 

p 
divided  among  Pyrex    culture  tubes  (6x50  mm)  each 

receiving  approximately  one  half  milliliter  of  sample. 

These  were  immediately  stored  at  -5  C. 


51 


Determinations  of  the  concentrations  of  DNA  after 
dissolution  and  after  dialysis  yielded  values  which  were 
identical,  within  experimental  error  and  also  within  +  5% 
of  the  calculated  concentration  based  on  the  manufacturer's 
data.   The  agreement  of  these  data  were  taken  as  evidence 
that  any  degradation  of  the  biopolymer  during  the  pro- 
cedures was  minimal. 

The  efficiency  of  the  dialysis  method,  used  here,  to 
replace  the  sodium  cations  of  the  original  DNA  salt  with 
the  cations  of  the  external  solution  was  evaluated  by 
atomic  absorption  spectroscopy.   A  Perkin-Elmer  model  290B 
atomic  absorption  spectrophotometer  outfitted  with  a  combina- 
tion sodium/potassium  lamp  was  used  to  measure  the  residual 
concentration  of  sodium  ion  in  a  sample  of  sodium  DNA  which 
had  been  dialyzed  against  a  potassium  phosphate  solution. 
Instrument  readings  were  taken  using  the  contents  of  the 
dialysis  bag  and  compared  to  those  obtained  for  standard 
solutions  in  which  the  ratio  of  potassium  to  sodium  ions 
was  approximately  one  thousand.   Note  that  in  this  instance 
high  relative  concentrations  of  potassium  ion  must  be  included 
in  the  standard  solution.   Less  than  0.1%  of  the  original 
sodium  remained  within  the  dialysis  bag  when  the  above 
procedure  was  employed. 

Crystalline  3-aminoacridine  was  prepared  by  Dr.  Timothy 
Roy  using  a  standard  procedure  (80)  except  for  additional 
purification  of  the  product  by  vacuum  sublimation.   Thin 
layer  chromatographic  analysis  of  the  purified  product 
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using  a  variety  of  solvents  and  solid  supports  failed  to 
indicate  the  presence  of  any  impurities.   The  spectra  of 
solutions  of  3-aminoacridine  at  pH  11.0  and  at  pH  2.5  had 
the  same  relative  absorptivities  as  those  reported  by 
Albert  (75).   Reported  (pH  11.0)  *™J  (log  e)  :  237  (4.46), 
262  (4.83),  321  (3.35),  337  (3.65),  353  (3.92),  410  (3.79). 
(pH  2.5)  233  (4.62),  274  (4.65),  349  (4.03),  365  (4.15), 
454  (4.10).   The  melting  point  agreed  with  that  reported 

(81):  218°C,  uncorr.,  223. 5°C,  corr. 

-3 
Stock  solutions  of  1  x  10   M  3-aminoacridine  were 

prepared  by  dissolving  2  mg  of  the  solid  in  7  ml  of  absolute 
ethanol.   The  solutions  were  stored  at  -5  C  in  glass.   One 
hundred  microliter  aliquots  of  the  ethanolic  solutions  were 
used  for  each  titration.   Periodic  tests  of  the  stock 
solutions  using  thin  layer  chromatography  showed  no  observ- 
able degradation  during  the  time  periods  in  which  the  samples 
were  being  employed. 

Absorptiometric  titrations  were  conducted  by  successive 
additions  of  microliter  aliquots  of  DNA  stock  solution  to 
10  ml  solutions  of  3-aminoacridinium  cation.   After  each 
addition,  the  absorbance  spectrum  between  550  nm  and  320  nm 
of  the  solution  was  recorded.   To  prepare  the  solutions,  10  ml 
of  supporting  electrolyte  at  a  given  concentration  were 
measured  in  a  10  ml  volumetric  flask.   This  volume  was  de- 
termined when  the  liquid  was  at  the  same  temperature  as  it 
was  to  be  during  the  titrations  (e.g.,  15.0,  25.0,  or  35.0  C) . 
Immediately  prior  to  a  titration,  0.1  ml  of  stock  3-amino- 
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acridine  in  ethanol  was  added  to  the  flask,  the  contents 
mixed  and  then  transferred  to  a  4.00  cm  quartz  cell  (Bolab 
Inc.,  Derry,  New  Hampshire).   The  cell  and  its  contents 
were  placed  in  the  thermostated  compartment  of  the 
spectrophotometer  (Beckman  model  DB-GT)  and  allowed  to 
become  thermally  equilibrated.   The  temperature  of  the 
instrument's  cell  compartment  was  maintained  at  15.0  + 
0.2°C,  25.0  +  0.05°C  or  35.0  +  0 . 2°C  by  circulating  thermo- 
stated water  (Aminco  refrigerated  water  bath)  through  its 
external  jacket. 

The  spectrum  of  the  sample  was  recorded  after  which 
an  aliquot  of  stock  DNA  solution  was  added.   Syringes  (Uni- 
metrics,  Inc.)  were  employed  to  dispense  volumes  of  DNA 
solution  ranging  from  4  yl  to  100  ul  while  volumes  greater 
than  100  Ul  (generally  added  toward  the  ends  of  the 
titrations)  were  measured  with  a  100  microliter  micro- 
pipette  (Centaur  Chemical,  Danbury,  Connecticut).   The 
individual  culture  tube  of  the  stock  DNA  solution  from 
which  aliquots  were  being  withdrawn  for  a  given  titration 
was  maintained  at  25°C,  regardless  of  the  temperature  at 
which  titration  was  being  conducted.   After  adding  the  DNA 
to  the  3-aminoacridinium  solution  in  the  cell,  the  solution 
was  mixed  by  withdrawing  about  5  ml  of  the  cell's  contents 
into  a  serological  pipet,  then  rapidly  forcing  them  back 
into  the  cell.   This  was  done  three  or  four  times  after 
each  addition  of  DNA. 
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Concentrations  of  the  stock  DNA  solutions  were 
determined  by  measuring  the  absorbance  of  a  solution 
containing  50  ml  of  stock  DNA  in  10.00  ml  of  0.010M 
supporting  electrolyte  at  pH  5.90  (1.00  cm  pathlength) . 
The  concentrations  of  at  least  ten  samples,  each  taken 

from  separate  tubes,  were  calculated  using  the  molar 

4-1   -1 
absorptivity  at  260  nra  of  1.413  x  10  H   cm   (75).   All 

glassware  was  cleaned  by  immersion  into  concentrated 
sulfuric  acid  saturated  with  potassium  dichromate.   Spec- 
trophotometer cells  were  rinsed  with  water  followed  by 
ethanol  and  were  stored  in  50%  concentrated  hydrochloric 
acid  in  ethanol  to  prevent  accumulative  absorption  of  the 
dye  or  of  DNA  to  their  surfaces. 

Solid  7-aminoquinoline  was  a  gift  from  Dr.  D.  Jackson, 
Texas  Tech.  University,  Lubbock,  Texas,  and  was  determined 
to  be  at  least  98%  pure  according  to  molar  absorptivity 
values  (82-84).   Reported  (pH  12)  X™    (log  e)  :  240  (4.56), 
275s  (3.58),  335  (3.68),  346  (3.64).   (pH  2.0):  260  (4.44), 
285s  (3.57),  392  (3.90).   The  reported  pKa  of  the  compound 
is  6.65  (84).   A  stock  3.4  x  10   M  solution  in  absolute 
ethanol  was  prepared  and  30  VI   aliquots  added  to  8.00  ml 
of  0.10M  KH„P04  buffer  solution  at  pH  5.9.   The  7-amino- 
quinolinium  solutions  were  titrated  at  15.0  +  0.2°C,  25.0  + 
0.05°C,  and  35.0  +  0.2°C  using  the  same  equipment  and  pro- 
cedures as  for  the  3-aminoacridinium  titrations,  vide 
supra. 
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Calculations 

The  molar  absorptivity  of  the  3-aminoacridinium-DNA 
complex,  £"_,„„,  at  368  nm  was  calculated  from  absorbance 
spectra  of  solutions  in  which  essentially  all  of  the 
compound  was  in  the  bound  form.   Complete  binding  (>  99% 
of  all  drug  present)  was  assumed  when  at  least  three  suc- 
cessive additions  of  excess  DNA  to  the  solution  resulted 
in  no  change  in  the  total  absorbance,  corrected  for 
dilution.   It  was  further  presumed  that,  under  these 
conditions,  all  absorbance  at  368  nm  was  due  to  the  complex 
and  no  3-aminoacridine  existed,  either  as  the  neutral  or 
monocationic,  unbound  form.   Also,  only  solutions  containing 
monovalent  metals  (and  tertamethylammonium  cation)  at  low 
ionicities  were  used  since,  for  other  systems,  the  large 
amounts  of  DNA  needed  to  drive  the  equilibrium  toward 
products  caused  serious  light  scattering  and,  in  some 
instances,  precipitation  of  DNA  salts.   Calculated  values 
of  £„„„,  in  the  presence  of  each  type  of  electrolyte  are 
presented  in  Table  1.   The  averaged  value  of  8494_  M  cm 
was  used  for  e_„_.  in  the  presence  of  all  supporting  elec- 
trolytes  at  all  concentrations.   A  similar  procedure  was 
used  to  evaluate  the  molar  absorptivity  of  the  7-amino- 
quinolinium-DNA  complex,  £,„„,  at  its  band  maximum  (405  nm) . 
At  pH  5.9  the  ratio  of  neutral  to  protonated  7-amino- 
quinoline  is  0.178.   It  was  assumed  that  only  the  cationic 
form  of  the  drug  binds  to  DNA.   The  averaged  molar  absorp- 
tivity of  the  7-aminoquinolinium-DNA  complex  was  evaluated 
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using  three  simultaneous  equations  and  was  found  to  be 
5332  M~1cm~1  at  405  nm. 

Concentrations  of  free  3-aminoacridinium  cation, 
[BH] ,  and  bound  monocation,  [BHP] ,  were  calculated  from 
absorptiometric  data  assuming  simple  additivity  of  the 
absorbances  of  the  two  species  at  the  analytical  wave- 
length 

A  =  eBH£[BH]  +  eBHp£[BHP]  (2-1) 

and  considering  the  mass  balance  expression 

Cfc  =  [BH]  +  [BHP]  (2-2) 

where 

A  =  the  total  absorbance  of  the  solution  at  a  given 
wavelength. 

£BH  =  the  molar  absorptivity  of  the  free  species, 
1.413  x  103  M~1cm~1. 

eBHP  =  the  molar  absorptivity  of  the  bound  complex, 
8.494  x  103  M~1cm~1. 
i   =   pathlength  of  light  in  the  cell. 
Ct  =  the  total  concentration  of  free  and  bound  species. 
[BH]  =  the  molar  concentration  of  free  3-aminoacridinium 
cation. 
[BHP]  =  the  molar  concentration  of  bound  3-aminoacridinium- 
DNA  complex. 
Solving  equation  (2-2)  for  [BH] ,  substituting  into  equation 
(2-1) ,  and  rearranging  yields 
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tBHPl  -  e        T-e      I      ■  <2"3> 

BHP    BH 


Before  any  DNA  has  been  added,  the  total  absorbance  of 
the  solution  will  be  due  entirely  to  free  BH  and,  since 
BH  at  this  point  will  be  equal  to  C  ,  the  total  absorbance 
before  any  titration,  A  ,  may  be  expressed  as 

Ao  =  eBHV   '  .        •  (2"4) 

Similarly,  when  all  of  the  compound  has  been  bound,  the 
final  absorbance,  A  ,  is  equal  to 


A-  =  £BHPCt£  <2"5> 


and  for  any  point  in  the  titration,  where  the  total  absor- 
bance is  A  , 

[BHP]  =  [(At-Ao)/(Aoo-Ao)]Ct  (2-6) 


and 


[BH]  =  [l-(At-Ao)/(Aoo-AQ)]Ct   .  (2-7) 


Since  an  appreciable  amount  of  the  total  7-amino- 
quinoline  exists  in  the  solution  in  the  neutral  form  at 
pH  5.9,  the  equations  representing  the  protolytic  equilibrium, 


K   -  £H+]  [A]  ,,  Q> 

a [AH]~  (2"8> 
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the  conservation  of  mass  for  7-aminoquinoline, 

Cfc   =     [A]    +     [AH]    +     [AHP]  ,  (2-9) 

and  the  photometric  absorbance  at  any  point  in  the 
titration  of  the  small  molecule  with  DNA, 

At  =  £A[A]£  +  eAR[AH]£  +  ?AHp[AHP]£  (2-10) 

must  be  employed  to  determine  the  concentrations  of  the 
species  involved  in  the  binding.   Concentrations  of  drug 
species  are,  explicitly, 

(£A-T7+  eAH+)  V  "V1  +  Ka/[H+]}. 
[AHP]  =  LE_J (2-11) 

{£A^-+  £AH  "  ?AHP(1  +  Ka/[H+])U 


At    £AHP  Ct* 
[AH]   =  (2-12) 

{£A-T7+    £AH   -   FAHP(1    +   Ka/[H+])U 


(At   "    £AHP    Ct£)    Ka[H+] 
[A]    =   — 4  (2-13) 

K  K 

{EA— r   +    £AH    "    £AHP(1    +   -T">>* 
[H    ]  Ati  AHF  [H    ] 
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(2-14) 


It  was  also  useful  to  directly  calculate  the  ratio, 

tAHP]    ^^  +  £AH+)  Ct"  "  AU  +  Ka/[H+]} 
[AH]  Afc  -  FAHp  Ctl 

where 


K  =  dissociation  constant  for  the  drug  or  probe, 
2.24  x  10~7  M-1. 
[H  ]  =  hydronium  ion  concentration,  1.26  x  10   M. 
£AHP  =  molar  absorptivity  of  the  bound  monocation  at 
the  specified  analytical  wavelength,  5.332  x 
103  M~1cm~1  at  420  nm. 
eA„  =  molar  absorptivity  of  the  free  monocation  at 
the  specified  analytical  wavelength,  7.943_  x 
103  M~1cm"1  at  420  nm. 

e     -   molar  absorptivity  of  the  free  base  at  the 

2   -1 
specified  analytical  wavelength,  3.47  x  10  M 

cm   at  392  nm. 
I   =  absorption  cell  path,  4.00  cm. 
A  =  the  absorbance  at  any  point  during  the  titration 

at  the  specified  analytical  wavelength. 
C   =  total  concentration  of  drug  or  probe. 

For  each  point  in  a  titration  the  total  absorbance  was 
determined  along  with  appropriate  constants  which  were 
substituted  into  the  necessary  relationships  to  yield  the 
concentrations  of  free  and  bound  species.   Raw  data  for 
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all  titrations  are  presented  in  Tables  2  through  17.   Note 
that  volume  corrections  for  the  addition  of  DNA  solution 
to  the  sample  were  made  in  all  cases. 

The  equilibrium  association  constants  for  the 
surface  and  intercalative  modes  of  binding  were  calculated 
for  each  point  in  the  absorptiometric  titrations  of  the 
small  molecules  with  DNA,  using  the  relationships  developed 
by  A.C.  Capomacchia  and  S.G.  Schulman  (85) . 


Ks  =  [BHP]/[BH]  ([Pt/3]-[BHP])  (2-15) 

and 

Kj  =  [BHP]/[BH]  ([Pt/2]-2[BHP])2  (2-16) 


where  K   and  KT  are  the  apparent  association  constants 
si  rsr 

for  the  surface  and  intercalative  modes,  respectively,  and 

P   is  the  total  molar  DNA  phosphate  concentration.   The 

equations  for  the  binding  of  7-aminoquinolinium  to  DNA 

are  the  same  except  [AHP]  and  [AH]  are  substituted  for 

[BHP]  and  [BH] ,  respectively.   The  extended  Debye-Huckel  equation 

(1-28)  was  used  to  calculate  activity  coefficients  at  the 

various  ionic  strengths,  for  charged  species  involved  in 

the  association  equilibrium. 

Apparent  association  constants  were  corrected  for 
deviations  of  molar  concentration  from  activity,  by 


aRHP 
K'  =  K  (2-17) 

n 

aBHap 
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where  K'  is  the  corrected  apparent  association  constant f 

anur./  ar,tJ/  and  a   are  activity  coefficients  for  bound 
BHP    BH        p 

species,  free  species,  and  DNA  phosphate,  respectively. 
The  value  of  n  is  1  for  surface  binding  and  2  for  inter- 
calative  binding.   For  the  binding  of  7-aminoquinolinium 
to  DNA,  the  subscripts  BHP  and  BH  are  replaced  with  AHP 
and  AH  corresponding  to  bound  and  free  7-aminoquinolinium, 
respectively.   These,  and  other  calculations,  were  done 
using  a  Litton-Monroe  model  Beta  326  Scientist  computer 
equipped  with  a  tape  cassette  drive  (Litton-Monroe, 
Orange,  New  Jersey) .   The  simultaneous  solution  of  (2-8) , 
(2-9) ,  and  (2-10)  to  obtain  concentrations  of  A,  AH,  and 
AHP  was  accomplished  using  a  program  supplied  by  the 
computer  manufacturer.   Programs  for  specific  calculations 
used  here  were  designed  and  are  presented  in  Appendix  III. 

Thermodynamic  parameters  for  the  binding  of  3-amino- 
acridinium  and  7-aminoquinolinium  to  DNA  were  evaluated 
using  equation  (1-31).   Values  of  log  K  at  15°C,  25°C,  and 
35°C  for  surface  and  for  intercalative  binding  of  each 
species  were  determined.   These  values  were  then  plotted 
vs.  1/T  (K) .   The  slopes  of  the  lines,  equal  to  AH°/2.303R, 
were  used  to  calculate  the  standard  entalpies  of  binding. 
Values  for  the  standard  entropies  of  the  reactions  were 
evaluated  using 

AS°  =  (AH°  -  AG°)/T  (2-18) 


CHAPTER  III 


RESULTS  AND  DISCUSSION 


General  Spectral  and  Titration  Characteristics 

A  number  of  interesting  spectral  features  and 
titration  characteristics  are  common  to  all  of  the 
systems  investigated.   Space  limitations  preclude  inclu- 
sion of  all  of  the  absorptiometry  titration  spectra, 
though  four  which  are  representative  are  presented 
(Figures  3  through  6).   These  are  spectra  of  titrations 
done  in  the  presence  of  the  background  cation  having  the 
lowest  charge  density,  Cs  ,  and  the  greatest  charge  density, 
Mg   ,  at  the  highest  and  lowest  concentrations  of  each. 
As  a  rule,  the  lower  the  charge   and  charge  density,  of 
the  background  cation,  the  less  effect  it  appears  to  have 
on  the  two  modes  of  binding.   The  highest  concentration  of 
the  smallest  divalent  cation  represents  the  other  extreme, 
as  it  severely  inhibits  binding  processes.   The  rest  of 
the  systems  are  intermediate  between  these  two.   In  each 
of  the  figures,  only  7  curves  are  drawn  for  the  sake  of 
clarity,  though  usually  between  18  and  22  spectra  were 
recorded  for  any  given  titration.   For  the  curves  ommitted 

from  figures  3  through  6,  the  absorbances  at  A362nm  0f  the 

max 
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L,  bands  may  be  obtained  from  the  appropriate  tables  in 
Appendix  II. 

In  all  cases,  addition  of  DNA  to  solutions  of  3- 

aminoacridinium  cation  resulted  in  decreases  (hypochromisms) 

.  .  1       1 

in  the  intensities  of  both  the   L,  and   L   bands.   It  is 

b       a 

noteworthy  that  even  small  amounts  of  DNA  (<  1  mole  DNA 
per  mole  of  drug)  added  to  3-aminoacridinium  solutions 
produced  only  hypochromic  shifts.   However,  addition  of 
small  amounts  of  DNA  to  other  compounds  such  as  7-amino- 
quinolinium  causes  an  initial  increase  in  the  total  intensity 
(hyperchromism)  of  the  drug's  spectrum.   Hypochromism  sub- 
sequently occurs  when  DNA  in  excess  of  about  1  mole  per  mole 
of  drug  is  added.   Presently,  it  is  not  known  what  factors 
are  responsible  for  the  initial  hyperchromism  of  some  species' 
spectra  or  why  some  exhibit  initial  hyperchromism  while 
others  exhibit  only  hypochromism  throughout  their  whole 
titration,  though  significant  variations  in  binding  behavior 
appear  highly  probable. 

It  may  be  seen  from  Figures  3  through  6  that  addition 
of  increasing  increments  of  DNA  to  3-aminoacridinium  results 
in  red  (bathochromic)  shifts  of  both  bands  of  the  drug's 

spectrum.   The  bathochromism  of  the  1L,'^ 1A  and  1Ll_^ ■'"A 

d  b 

bands  are  a  result  of  decreases  in  the  energy  separations 
between  the  ground  and  excited  states.   This  may  be  due  to 
a  lowering  of  the   L  states,  a  raising  of  the  ground,  state, 
or  a  combination  of  both.   If  binding  of  small  molecules  to 
DNA  occurs  as  the  reaction  proceeds,  then  it  is  reasonable 
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to  hypothesize  that  the  cationic  drug  is  moving  from  a 

highly  polar,  high  dielectric  medium  (water)  to  a  less 

polar,  more  lipophilic  solvent  such  as  the  region  in,  and 

around,  the  DNA  helix.   An  environmental  change  of  this 

type  would  tend  to  favor  a  decrease  in  charge  separation 

* 
within  a  molecule  and,  thus,  favor  it  — >  it   transitions 

in  which  the  formal  positive  charge  of  the  heterocyclic 
ring  nitrogen  is  decreased.   Such  a  case  would  mean  that 
a  lowering  of  the  energies  of  the  excited  states  is  the 
predominant  factor  in  causing  the  red  shift  upon  binding 
if,  indeed,  such  "solvent  effect"  arguments  are  valid  when 
considering  binding  processes.   Caution  must  be  exercised 
in  this  regard  as  Capomacchia  (69)  has  found  that,  for  a 
number  of  compounds,  spectral  shifts  upon  binding  to  DNA 
are  not  in  qualitative  agreement  with  shifts  observed  when 
the  compounds  were  transferred  from  polar  aqueous  solvents 
to  less  polar  organic  solvents.   Maximum  shifts  of  the  two 
bands  are  from  362  nm  to  368  nm  for  the   Lfa  band  and  from 
4  54  nm  to  463  nm  for  the  """L  band.   These  differences  (6  nm 
and  9  nm,  respectively)  correspond  to  decreases  in  the 
transition  energies  of  450  cm-1  (1.29  Kcal/mole  dye)  and 
428  cm-1  (1.23  Kcal/mole  dye).   The  similar  magnitudes  show 
that,  with  excess  DNA  at  low  total  ionic  strength,  the  two 
transitions  are  equally  affected  from  an  energetic  stand- 
point in  the  presence  of  excess  DNA.   However,  a  large 
number  of  titrations  have  substantiated  the  observation 

that  the  1L,_  moves  only  slightly,  from  362  nm  to  =  364  nm, 
fc> 
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over  most  of  the  titration  and  then  abruptly  shifts  from 

s  364  nm  to  368  nm  while  the  shift  of  the   L  band  is 

smooth  throughout  the  whole  titration.   The  abrupt  shift 

of  the  1L1  takes  place  at  approximately  the  point  in  the 
b 

titration  where  any  additional  DNA  results  in  no  further 
decrease  in  1L  band  intensity.   In  fact,  this  phenomena 

D 

has  been  accepted  as  evidence  that  all  of  the  protonated 
3-aminoacridine  has  been  bound.   Moreover,  the  rapid  shift 
of  the  higher  energy  band  always  occurs  significantly 
beyond  the  point  in  the  titration  where  very  little  change 
is  noted  in  band  intensity  or  position  of  the  maximum  of 
the  longer  wavelength   Lb  band.   These  observations  suggest 
that  the  environment  appears  to  affect  both  transitions 
equally  in  the  presence  of  excess  DNA,  but  they  are  affected 
differently  during  the  course  of  a  titration.   A  cursory 
consideration  of  the  data  may  lead  one  to  postulate  that 
the  abrupt  shift  is  due  to  intercalation  of  the  molecule. 
This  is  plausible  since  internal  binding  would  affect  the 
1L,  more  dramatically  than  the  Lfl/  especially  if  the  charged 
ring  nitrogen  remained  exterior  to  the  helix.   Further 
consideration,  however,  reveals  substantial  evidence,  both 
in  this  work  and  in  the  literature,  which  shows  that  inter- 
calation occurs  much  earlier  in  the  titration  and  that  it  is 
much  less  abrupt  a  process  than  the  bathochromic  shift 
considered  here  would  require.   The  abrupt  change  might 
arise  from  some  form  of  near-neighbor  interaction  of  bound 
drug  molecules  (63) .   During  most  of  the  titration,  where 
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the  total  DNA  phosphate  to  total  3-aminoacridine  ratio, 
P/D,  is  low,  the  bound  molecules  can  interact  with  one 
another.   But  as  the  ratio  increases,  bound  species  become 
isolated  from  each  other  to  the  point  of  no  longer  being 
able  to  interact.   The  hypochromic  shift  is  not  incon- 
sistent with  this  proposal  as  both  red  and  blue  shifts 
are  possible  upon  exciton  formation.   Another  factor  to 
be  considered  is  that  as  the  DNA  concentration  is  increased 
aggregation  of  the  DNA,  itself,  may  begin  to  occur.   At 
some  critical  level,  the  bound  dye  may  experience  a  marked 
decrease  in  the  hydrophilicity  of  its  environment,  perhaps 
as  a  result  of  the  neighboring  polymers  excluding  water 
between  themselves.   Under  these  circumstances,  red  shifts 
would  be  expected,  though  why  they  would  be  seen  only  for 

the   L,  and  not  for  the   L   transition  is  unanswered. 
b  a 

The  onset  of  the  rapid  red  shift  appears  to  be  depen- 
dent on  the  total  ionic  strength  of  the  medium  as  well  as 
the  type  of  cations  in  solution.   In  the  presence  of 
large,  diffusely  charged,  monovalent  cesium  a  distinct 
bathochromism  occurred,  even  in  the  highest  concentration 

(0.15M)  at  the  moderate  P/D  ratio  of  65.   On  the  other 

-4     ++ 
hand,  6.3  x  10   N  Mg   was  sufficient  to  prevent  the  shift 

to  368  nm  even  when  phosphate  is  in  excess  of  drug  by  a 

factor  of  65.   In  0.0025M  CsH2PO. ,  a  P/D  ratio  of  only  14 

results  in  the  spectral  shift  (in  fact,  it  probably  occurs 

at  about  P/D  =  11) .   Charge  density  of  the  countercations 

influenced  the  point  at  which  the  abrupt  red  shift  was  seen 
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in  a  manner  parallel  to  the  degree  of  completion  of  a 
titration  for  a  given  P/D  ratio.   The  overall  spectral 
behavior  of  the  systems  in  the  presence  of  the  alkaline 
earths  was  the  same  as  for  the  alkalies  except  the 
maximum  concentrations  of  M   that  allowed  at  least  90% 
binding  at  P/D  <  100  was  much  lower  than  for  alkali 
metals.   Concentrations  of  group  Ila  cations  in  excess 
of  about  0.02N  resulted  in  precipitates  forming  in  the 
samples.   These  were,  presumably,  the  metal  salts  of  DNA 
and  appeared  long  before  sufficient  polymer  was  added  to 
appreciably  bind  the  drug. 

Throughout  this  discussion,  the  absorbances  at  the 
band  maxima  have  been  assumed  to  be  directly  related  to 
the  areas  under  the  curves.   Since  it  is  the  areas,  and 
not  peak  heights  at  the  maxima  which  are  true  measures  of 
transition  probabilities,  care  must  be  taken  in  equating 
peak  heights  with  relative  transition  probabilities.   If 
the  overall  geometry  of  the  L,  band  were  to  change  during 
the  course  of  a  titration,  peak  heights  would  not  be  valid 
representations  of  the  progress  of  the  titration.   In  our 
studies,  inspection  of  the  bands  showed  that  their  geome- 
tries were  constant  during  the  titrations,  hence,  we  feel 
approximating  peak  heights  with  peak  areas  was  valid. 

In  addition,  we  have  assumed  that  the  0-0  vibronic 

band  maximum  of  the   L,  < A  transition  is  coincident 

b 

with  the  maximum  of  the  total  transition  envelope.   This 
may  be  misleading  if  the  maximum  of  the  envelope  arises 


■ 
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from  the  sums  of  the  intensities  of  vibronic  transitions 
which  are  close  enough  together  to  appreciably  overlap 
at  the  maximum  wavelength.   If  the  vibronic  energies  are 
altered  in  dissimilar  ways  during  a  titration,  then  the 
maximal  envelope  intensity  will  reflect  this  dissimilar 
change  along  with  any  real  shift  in  the  0-0  band.   It  is 
possible,  therefore,  that  shifts  of  the  1L  envelope 
maximum  may  be  due  to  either  a  disproportionate  change  in 
the  magnitudes  of  the  band  maxima  of  the  0-0  and  0-1  vibronic 
bands  or  a  bona  fide,  pure  0-0  shift. 

Variations  in  the  spectral  characteristics  of  the  3- 
aminoacridinium-DNA  systems  in  the  presence  of  metal  phos- 
phates and  metal  acetates  were  not  due  to  the  different 
counteranions,  as  evidenced  by  the  results  of  titrations 
using  KH2P04  and  K02CCH3  as  background  electrolytes.   Re- 
sults obtained  with  the  two  buffer  ions  were  the  same  in 
all  respects  except  at  very  low  ionic  strengths.   Deviations 
in  this  region  are  probably  due  to  lack  of  sufficient  pH 
control  or  alterations  in  the  conformation  of  DNA  possibly 
arising  from  mild  denaturation. 

The  Lfa  envelope  of  the  free  3-aminoacridinium  cation 
consists  of  two  distinct  vibronic  bands  having  maxima  at 
348  nm  and  362  nm,  corresponding  to  the  0-1  and  0-0  vibronic 
transitions,  respectively.   As  a  titration  proceeds,  the  fine 
structure  of  the  envelope  is  lost  as  a  result  of  the  348  nm 
peak  becoming  increasingly  less  pronounced,  to  the  point  of 
becoming  a  mere  perturbation  of  the  overall  envelope.   It  is 
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difficult  to  determine  the  exact  position  of  the  shoulder 

when  the  titration  is  nearly  complete,  but  it  appears  that 

its  shift  is  essentially  the  same  as  that  of  the  362  nm 

band.   In  cases  where  the  abrupt  364  nm  to  368  nm  shift 

of  the  v __  vibronic  band  of  the   L,  occurs,  there  is  an 
0-0  b 

accompanying  shift  of  the  shoulder  maximum  from  a  350  nm  to 
about  354  nm  with  a  slight,  but  definite  increase  in  the 
resolution  of  the  two  bands.   The  degree  of  fine  structure 
of  a  compound's  spectrum  provides  qualitative  information 
regarding  the  environment  of  the  species.   Loss  of  fine 
structure  may  be  rationalized  on  the  basis  of  increasing 
the  degrees  of  vibrational  freedom  of  the  molecule.   If, 
upon  binding,  the  energy  separation  between  the  v  =  0  and 
v  =  1  vibrational  levels  of  the  ground  state  were  to  be 
reduced,  a  coalescing  of  the  two  peaks  would  be  anticipated. 
Such  a  situation  may  arise  when  a  molecule  moves  from  an 
environment  in  which  its  vibrational  motions  are  restricted 
to  one  in  which  they  are  less  so.   For  this  to  maintain  in 
our  systems,  the  bound  species  would  have  to  have  greater 
vibrational  freedom  than  the  wholly  water-solvated  free 
cation.   This  is  feasible  considering  the  strong  degree  of 
interaction  between  water  molecules  and  any  dipole  in  their 
midst  compared  to  the  lipophilic,  noninteracting  "solvent" 
of  the  surface  and  interior  of  DNA. 

The  increase  in  fine  structure  which  accompanies  the 
dramatic  red  shift  of  the  0-0  band  of  the  L,  envelope  at 
high  DNA  concentrations  is  not  well  understood  but  serves 
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to  corroborate  evidence  for  a  dramatic  environmental  change 

of  bound  drug  at  these  very  high  levels. 

To  summarize,  titrations  of  3-aminoacridinium  cation 

with  DNA  salts  in  the  presence  of  monovalent  and  divalent 

metal  ions  (including  tetramethylammonium  ion)  of  varying 

concentrations  result  in  hypochromism  only.   In  all  cases, 

both  the  L  and  L,  envelopes  of  the  free  drug's  spectrum 

red  shift,  at  least  to  some  extent,  during  titration  — 

the   L   from  454  nm  to  463  nm  and  the   L,  from  362  nm  to 
a  b 

about  364  nm.   In  the  presence  of  all  concentrations  of 
supporting  electrolytes  having  large,  monovalent  cations, 
a  further,  sharp  red  shift  from  364  nm  to  368  nm  takes  place 
when  a  P/D  ratio  of  10  to  15  is  reached.   The  shift  may  be 
effected  in  the  presence  of  smaller  alkalimetal  ions  at 
P/D  ratios  approaching  70.   Regardless  of  the  concentration 
or  charge  density  of  divalent  cations,  the  dramatic   L, 
shift  did  not  occur.   Loss  of  vibronic  fine  structure  of 
the  L,  envelope  was  observed  in  all  cases  as  the  titrations 
progressed,  up  to  the  point  at  which  the  dramatic  batho- 
chromic  shift  occurred.   Concomitant  with  the  rapid  shift 
was  a  reemergence  of  the  0-1  vibronic  band. 

Using  relationships  developed  in  the  Experimental 
section,  concentrations  of  free  and  bound  3-aminoacridinium 
cation  were  computed  from  the  total  absorbances  of  solutions 
which  contained  varying  amounts  of  DNA.   Data  for  all  of 
the  titrations  are  not  included  here,  though  some,  which 
include  the  extremes  of  experimental  conditions,  are  presented. 
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Tables  18  through  21  contain  tabulations  of  [BH]  and 
[BHP]  along  with  log  ([BHP]/[BH])  values  for  reactions  in 
0.15M  and  0.0025M  CsH2PC>4  and  in  0.010N  and  6.3  x  10~4N 
Mg(02CCH3)2.   In  lower  concentrations  of  alkali  metals 
and  tetramethylammonium  ion,  the  addition  of  as  little 
as  5  pi  of  9  x  10   M  DNA  phosphate  solution  results  in  a 
significant  degree  of  binding  of  the  drug.   For  instance, 
nearly  10%  of  the  total  3-aminoacridinium  present  in 
0.0025N  CsH2P04  is  bound  after  addition  of  only  5.0  pi  of 
DNA.   Since  volumes  of  DNA  less  than  5  yl  cannot  be 
measured  with  much  reliability  due  to  its  viscosity,  the 
ratio  of  bound  to  free  drug  of  0.1  should  be  considered 
as  the  effective  minimum  for  titrations  of  this  sort.   When 
lower  ratios  were  possible,  such  as  in  solutions  containing 
divalent  cations  or  >  0.010M  monocations,  random  experi- 
mental error  was  sufficiently  high  to  make  the  data  unreliable, 
As  a  result,  our  studies  indicate  that  absorptiometry 
titrations  of  small  molecules  with  DNA,  under  conditions 
similar  to  ours,  are  probably  invalid  whenever  the  ratio  of 
bound  to  total  drug  is  less  than  0.1.   It  should  be  recalled 
that  our  use  of  4.0  0  cm  cells  maximized  sensitivity,  thereby 
allowing  minimal  increments  of  DNA  solution  to  be  used. 
Smaller  amounts  of  DNA  could  be  added  in  5  yl  increments  by 
employing  a  more  dilute  stock  solution,  but  the  errors  arising 
from  other  sources  would  likely  be  unimproved.   At  the  other 
titration  extreme,  ratios  of  bound  to  free  drug  in  excess  of 
about  0.95  also  result  in  unacceptable  randomness  in  most 
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cases.   Our  procedures  are  designed  to  minimize  the  impact 
of  these  sources  of  imprecision  on  the  calculated  association 
constants. 

The  maximum  extent  of  total  binding  which  could  be 
achieved  varied  with  the  character  and  concentration  of 
the  supporting  electrolyte.   Table  22  lists  the  percent 
of  total  drug  bound  to  DNA  after  the  final  increment  of 
DNA  was  added  in  titrations  done  in  various  media.   The 
P/D  values  listed  are  the  maximum  that  could  be  obtained 
due  to  solubility  and  pH-control  limitations.   No  effort 
was  made  to  determine  how  the  drug  was  bound  or  what  per- 
centage of  the  total  was  surface  bound  compared  to  inter- 
calatively  bound.   As  discussed  previously,  for  those  titra- 
tions in  which  an  abrupt  bathochromic  shift  of  the   L,  band 
was  observed,  the  binding  was  considered  to  be  complete 
(100  +  2%).   The  averaged  molar  absorptivity,  £„„„,  determined 
from  these  titrations  was  used  to  compute  concentrations  of 
free  and  bound  species  for  all  systems.   In  light  of  the 
experimental  error  incurred  in  evaluating  c  and  the 

simple  fact  that  thermodynamic  equilibrium  demands  that 
there  cannot  be  total  binding,  percentages  greater  than 
about  96%  are  taken  as  indicative  of  complete  binding.   From 
Table  22,  it  may  be  seen  that  at  least  97%  of  the  total  3- 
aminoacridinium  present  is  bound  in  solutions  of  <  0.010M 
MH-PO..   To  attain  essentially  complete  binding  in  these 
systems  only  about  10  to  20  times  as  much  DNA  phosphate  as 
total  drug  need  be  present  (see  P/D  column  of  Table  22) .   As 
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expected,  larger  and  larger  ratios  of  DNA  to  drug  are 
necessary  for  complete  binding  as  the  concentration  of 
a  given  electrolyte  is  increased.   Similarly,  it  becomes 
increasingly  more  difficult  to  shift  the  equilibrium  wholly 
toward  bound  product  in  the  order  (of  electrolyte  cations) : 
(CH3)4N+  <  Cs+<Rb+<K+<Na+.   A  14-fold  excess  of  DNA  was 
required  for  97%  binding  in  0.0025N  Cs+;  17-fold  was  needed 
for  97%  binding  in  0.010M  Cs;  while,  even  with  a  65-fold 
excess  of  DNA,  only  90%  of  the  drug  was  bound  in  0.15M 
CsH2P04.   It  was  relatively  easy  to  totally  bind  the  3- 
aminoacridinium  cation  in  low  concentrations  of  monovalent 
cations  but  very  difficult  to  attain  greater  than  90% 

binding  in  even  the  lowest  concentrations  of  the  divalent 

-4    ++ 
electrolytes.   In  6.3  x  10   N  Ba    (the  dication  having 

the  lowest  charge  density  of  those  investigated)  only  90% 
of  the  drug  was  bound  at  a  P/D  ratio  of  47  and,  for  Mg++ 
at  the  same  concentration,  90%  was  bound  at  P/D  =  65. 
Similarly,  in  0.010N  Ca++  and  0.010N  Mg++  the  maximum 
degrees  of  binding  were  77%  and  70%,  respectively,  with 
P/D  ratios  of  103  and  104.   The  choices  of  the  upper  limits 
of  total  phosphate  concentrations  reported  here  were  de- 
pendent on  the  solubility  of  species  in  the  concentrated 
electrolytes  and  on  pH  considerations  for  the  dilute  systems. 
Visible  precipitation  or  erratic  spectral  behavior  were 
evident  in  solutions  having  P/D  values  greater  than  those 
listed  in  Table  2  2  for  barium  and  magnesium.   The  maximum 
concentration  of  the  divalent  electrolytes  was  set  at  0.010N, 
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since  obvious  precipitation  occurred  in  solutions  of 

2+ 
0.025N  M   well  before  the  titrations  were  complete  (i.e., 

at  30%  to  50%  binding) .   Acetates  of  the  alkaline  earths 
were  used  instead  of  their  phosphates  due  to  the  latter' s 
very  low  solubilities.   Despite  dialysis  of  DNA  solutions 
against  pH  5.90  buffer,  their  pH  levels  were  still  ap- 
parently close  to  7.   Thus,  addition  of  moderate  volumes 
of  DNA  to  the  very  poor  acetate  buffer  systems  caused 
excessive  increases  in  the  pH  levels  of  the  dilute  buffer 

solutions.   Consequently,  acetate  solutions  having  more 

-4 
than  approximately  5  x  10   M  DNA  phosphate  had  pH  values 

considerably  in  excess  of  5.90. 

Upon  perusing  the  second  column  of  Table  22,  it  becomes 
quickly  evident  that  the  extent  of  binding  is  not  predomi- 
nately and  directly  related  to  ionic  strength.   Total  ionic 
strengths  were  computed  for  the  systems  having  the  P/D 
values  listed  using  equation  (1-29)  and  Program  4  of 
Appendix  III.   Notice  that  while  the  ionic  strengths  of 
0.0025M  CsH2P0.  and  6.3  x  10~4N  Ba(02CCH3)2  are  the  same, 
the  binding  is  complete  in  the  former  and  only  90%  complete 
in  the  latter.   This  shows  that  the  effect  of  the  formal 
charge  on  the  cations  of  the  electrolyte  is  greater  than 
predicted  on  the  basis  of  the  simple  electrostatic  effects 
incorporated  in  the  limiting  Debye-Ruckel  equation.   That 
the  contribution  of  excess  DNA  to  the  total  ionic  strength 
is  not  responsible  for  these  large  differences  may  be  seen 
by  comparing  0.010N  barium  and  0.010N  magnesium  solutions 
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which  both  contain  the  same  amount  of  DNA  and  have,  there- 
fore, the  same  total  ionic  strengths.   The  extent  of 
binding  in  the  former  is  77%  while  in  the  latter  it  is 
only  70%. 

Rapid  estimations  of  the  values  of  m  and  q  of  equation 
(1-21)  may  be  obtained  by  selecting  reasonable  terms  such 
that  mq  =  n  (n  is  total  phosphate  per  drug) .   It  has  been 
shown  that  (85)  in  a  medium  containing  0.002M  to  0.005M 
KH»PO.  (pH  5.9)  the  binding  of  3-aminoacridinium  cation  to 
DNA  is  best  expressed  by 

+        3"  2- 

BH   +  (P04)3    v   ^    BH(P04)3  (3-1) 

for  the  surface  mode,  and 

BH+  +  2(P04)2~  ^=±     BH((P04)2)2~  (3-2) 

for  the  intercalative  mode.   In  equation  (3-1) ,  three  DNA 
phosphates  act  as  a  single  moiety,  reacting  with  the  small 
molecular  cation  to  form  the  surface-bound  species.   In 
equation  (3-2) ,  2  phosphates  act  as  a  single  reacting 
entity  with  two  of  these,  in  turn,  interacting  with  one 
cation  to  form  the  intercalated  complex.   The  values  of  m 
and  q  are  3  and  1,  respectively,  for  the  first  case  and  2 
and  2  for  the  second.   Substituting  these  into  equation  (1-21) 
yields 


Ks  =  [BHP]/[BH]([Pt/3]-[BHP]))  (3-3) 

and 
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Kj  =  [BHP]/([BH]  ([Pt/2]-2[BHP])2)   .  (3-4) 

Kg  and  Kj  are  the  apparent  association  constants  for  the 
surface-bound  and  intercalatively-bound  drug-DNA  complexes, 
respectively.   These  equations  were  tested  to  certify  their 
validity  under  our  experimental  conditions.   Recall  that 
plots  of  log( [BHP]/[BH])  vs.  log ( [Pt/ml -q [BHP] )  yield  lines 
whose  slopes  are  equal  to  q  and  whose  intercepts  are  equal 
to  log  K.   Tables  18  through  21  include  calculated  values 
of  log([BHP]/[BH]) ,  log([Pt/3]-[BHP]),  and  log ( [Pt/21 -2 [BHP] ) 
for  representative  titrations.   Figures  7  and  8  are  plots 
of  log([BHP]/[BH])  vs.  log ( [Pt/31 - [BHP] )  for  the  reactions 
in  the  presence  of  0.0025M  Cs+  and  6.3  x  10_4N  Mg++;  and  in 
0.15M  Cs+,  0.010N  Mg++,  and  0.010M  Cs+,  respectively. 
Figures  9  and  10  are  plots  of  log ( [BHP] / [BH] )  vs.  log ( [P  ,  ]- 
2  [BHP])  for  the  same  systems.   Three  general  features  are 
particularly  significant.   One,  the  shapes  and  slopes  of 
the  curves  are  essentially  the  same  whether  log ( [P  .  ]-[BHP]) 
or  log([Pt <2]-2[BHP] )  is  plotted  against  log ( [BHP]/ [BH] ) . 
(Calculated  slopes  are  tabulated  in  Table  2  3.)   Only  in  the 
latter  portions  of  the  titrations  done  in  the  more  dilute 
alkali  metal  electrolytes  does  this  rule  begin  to  break  down. 
Two,  except  for  0.15M  CsH2P04,  all  plots  of  the  data  in- 
volving the  presence  of  monovalent  cations  exhibit  definite 
curvature  —  the  lower  the  concentration,  and  further  down 
a  series  of  the  periodic  chart,  the  greater  the  curvature. 
Conversely,  all  concentrations  of  all  M2+  salts  yield 
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essentially  straight  lines  over  the  whole  titration 
region.   Three,  in  those  cases  in  which  there  is  no 
curvature  throughout  the  titration  range,  the  slopes  of 
the  lines  are  all  nearly  1  and  parallel  the  initial 
regions  of  the  other  titrations. 

That  the  values  of  q  evaluated  from  both  types  of 
plots  are  much  the  same  in  most  media  is  evidence  of  the 
validity  of  this  technique.   It  appears  that  going  from 
q=ltoq=2or3  does  not  substantially  affect  the 
values  of  log ( [Pfc  .  ]-q[BHP] )  except  at  the  extrema  of  the 
titrations.   With  our  procedure,  however,  data  very  early 
and  very  late  in  a  titration  need  not  be  relied  upon  when 
estimating  lines,  unlike,  for  instance,  the  Scatchard 
method.   For  cases  where  the  background  electrolyte  is 
especially  weak  (i.e.,  low  concentrations  of  cations  having 
low  charge  density)  differences  in  the  slopes  of  the  two 
plots  become  appreciable  for  the  later  portions  of  the 
titrations.   For  example,  the  final  slopes  where  0.010M 
Cs  was  employed  as  the  background  electrolyte  are  1.68  and 
1.88  for  q  =  2  and  q  =  1,  respectively,  and  for  0.0025M  Cs+, 
the  values  are  1.71  and  2.16.   Several  explanations  may  be 
advanced  for  this  lack  of  agreement,  including  the  fact 
that  the  relatively  short  overall  range  of  the  amount  of 
DNA  needed  to  complete  the  titrations  magnifies  the  error 
inherent  at  the  extrema.   Also,  the  high  degree  of  curvature 
makes  the  determination  of  the  two  separate  straight  line 
regions  more  subjective.   Perturbations  of  the  reaction 
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system  which  are  significant  only  in  low  ionic  strength 
media  (e.g.,  mild  polymer  denaturation)  would  be  expected 
to  affect  both  types  of  plots  in  the  same  manner.   Re- 
gardless of  the  sources  of  the  discrepancies,  the  pro- 
cedure is  still  useful,  as  both  plots  yield  the  same 
integer  values  of  q.   Final  selection  of  m  and  q  is  made 
on  the  basis  of  the  constancy  of  the  computed  values  of 
the  equilibrium  association  constant  (vide  infra) . 

The  presence  of  at  least  two  distinct  linear  regions 
of  the  plots  for  lower  concentrations  of  alkali  metals 
strongly  indicates  two  binding  stoichiometries.   Numerous 
investigators,  as  outlined  in  the  Introduction,  have  shown 
the  presence  of  two  types  of  binding  processes  which  likely 
account  for  the  different  values  of  q.   What  has  not  been 
previously  reported  is  the  apparent  change  in  the  ratio  of 
phosphate  to  3-aminoacridinium  cation  on  going  from  one 
type  of  binding  to  the  other  as  evidenced  by  these  plots. 
Also,  as  the  ionic  strength  is  increased,  the  extent  of 
intercalative  binding  apparently  decreases.   This  is  mani- 
fested in  Figures  7  through  10  by  a  progressive  trend 
toward  a  wholly  linear  plot  as  the  ionic  strength  increases. 
That  is,  compare  the  shapes  of  the  curves  for  systems  con- 
taining 0.0025M,  0.010M,  and  0.15M  CsH2P04 .   Similarly,  only 
one  stoichiometry  is  evident  from  the  plots  for  M2+-containing 
systems,  which  are  linear  throughout  their  entire  titration 
regions.   Slopes  of  1  +  .  1  for  all  of  these,  including  that 
for  the  reaction  in  the  presence  of  1  x  10~5N  Ba(0  CCH  ) 

•—     s>     z. 
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indicates  that  only  surface  binding  occurs  in  the 
presence  of  the  divalent  metal  cations  over  the  titration 
range  investigated. 

Previously,  it  was  brought  out  that,  in  the  presence 
of  low  concentrations  of  monovalent  cations  as  background 
electrolytes,  plots  of  log ( [BHP]/ [BH] )  vs.  log ( [P  .  ] -q [BHP] ) 
exhibit  two  regions  of  linearity.   During  the  early  part  of 
a  titration  the  slope  of  the  lines  is  approximately  1  while 
in  the  later  region  the  slope  is  about  2.   As  the  ionic 
strength  is  increased,  or  as  the  charge  density  of  the 
cation  of  the  electrolyte  becomes  greater,  the  slope  of 
the  later  region  decreases,  eventually  becoming  equal  to 
about  1.   Slopes  of  the  initial  portions  are  essentially 
unaffected  by  alterations  in  the  supporting  media  in  which 
the  titrations  are  conducted. 

The  ramifications  of  these  relationships  may  be  explored 
by  first  defining  the  binding  interaction  in  terms  of  con- 
ventional metal-ligand  chemistry.   We  shall  analogize  the 
cationic  3-aminoacridinium  to  the  metal  center  and  the  DNA 
phosphates  to  ligands  which  bind  to  the  center.   For  clarity, 
we  shall  further  define  two  types  of  ligands,  having  indepen- 
dent and  variable  affinities  for  the  cation  centers;  a  weakly 
bound  ligand  (corresponding  to  surface  binding)  which  forms 
a  1:1  complex,  and  a  more  strongly  bound  ligand  (corresponding 
to  intercalative  binding)  which  forms  a  1:2  metal: ligand 
complex.   These  definitions  are  based  on  the  values  of  q  for 
surface-bound  and  intercalatively-bound  ligand  types  being  1 
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and  2,  respectively.   Finally,  we  shall  assume  that  the 
ligands  compete  for  metal  binding  sites  such  that  if  a 
cation  is  bound  in  the  surface  mode,  then  it  is  unavailable 
for  intercalative  binding,  and  vice  versa. 

During  the  initial  stages  of  a  titration  there  is  a 
large  excess  of  3-aminoacridinium  over  DNA  phosphate  or, 
by  analogy,  an  excess  of  metal  over  ligand.   In  such  a 
situation,  there  would  be  minimal  competition  for  metal 
sites  and  both  forms  of  binding  would  be  expected  to  occur 
to  their  fullest  extent,  based  on  their  respective  association 
constants.   We  may  view  the  titration  of  the  cation  with  DNA 
ligand  as  a  decrease  in  the  metal  center  concentration  as 
the  process  proceeds.   Addition  of  aliquots  of  DNA  serves 
to  limit  the  number  of  sites  available  for  binding,  causing 
the  more  weakly  bound  surface  species  to  be  affected  first. 
At  this  stage,  the  strongly  bound  intercalative  species  is 
relatively  unaffected.   This  argument  presupposes  that  the 
binding  affinity  of  the  intercalative  mode  is  substantially 
greater  than  that  of  the  surface  type.   Hence,  the  first 
region  of  the  titration  reflects  surface  binding  of  3- 
aminoacridinium  to  DNA.   As  the  titration  progresses,  the 
relative  number  of  available  binding  sites  becomes  limited 
to  such  an  extent  that  the  equilibrium  involving  the  inter- 
calated complex  is  the  one  predominantly  affected.   Now  the 
plot  of  log([BHP]/[BH])  vs.  log ( [Pt/J -q [BHP] )  is  a  reflection 
of  equilibrium  concentrations  governed  by  the  stronger  mode 
which,  in  turn,  results  in  the  slope  of  the  line  being  2 


instead  of  1.   Note  that  the  difference  between  the 
relative  binding  affinities  of  the  two  modes  must  be 
sufficiently  great  to  allow  one  to  predominate  over  the 
other  in  the  initial  and  final  portions  of  the  titrations. 
On  the  other  hand,  if  they  are  too  disparate,  only  one 
would  be  seen  throughout  the  experimentally  accesible 
titration  range.   Increasing  the  ionic  strength  of  a 
reaction  medium  causes  electrostatic  contributions  to  the 
total  free  energy  to  become  less  exothermic  for  a  process  in 
which  oppositely  charged  moieties  coalesce.   Based  on  this, 
decreases  in  the  apparent  association  constants  for  both 
surface  and  intercalative  binding  would  be  anticipated.   It 
is  assumed,  at  this  stage  of  the  argument,  that  changes  in 
electrostatic  free  energies  are  the  predominant  affects 
accompanying   alterations  in  ionic  strength.   Extrapolation 
of  the  straight  lines  of  Figures  7  through  10  to  the  vertical 
axis  shows  that,  as  the  ionic  strength  is  increased,  both 
K   and  K,  do  become  less  positive  in  agreement  with  the  above 
prediction.   However,  such  extrapolations  may  be  misleading 
in  the  intercalative  cases.   Recall  that,  upon  altering  the 
ionic  strength,  the  slopes  of  the  initial,  surface  binding 
regions  of  the  plots  are  essentially  unchanged  and  that  they 
form  a  family  of  lines.   This  constancy  of  slope  indicates 
that  there  are  no  changes  in  the  surface  binding  reaction 
(aside  from  equilibrium  shifts)  as  the  ionic  medium  is 
altered.   However,  two  factors  contribute  to  the  decrease  in 
the  intercept  values  of  intercalative  plots:  (1)  lower  values 
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of  log ( [BHP]/[BH] )  for  any  given  point  along  the  abcissa 

as  ionic  strength  is  increased,  and  (2)  the  decrease  in 

the  slopes  of  the  lines  as  the  ionic  strength  is  increased. 

Two  explanations  for  the  change  in  slope  may  be  advanced. 

The  value  of  q  may  be  indeed  changing  from  2  to  1 ,  for  the 

intercalation  reaction,  as  the  ionic  strength  is  increased. 

If  this  is  true,  then  either  the  overall  stoichiometry  is 

changing  such  that  mq  ^  4  or,  the  stoichiometry  is  remaining 

constant  and  m  is  becoming  equal  to  4  while  q  is  decreasing 

to  1.   A  more  likely  explanation  can  be  advanced  on  the  basis 

of  relative  changes  of  K   and  KT .   If  K   and  KT  do  not  change 

"         S        I         S        I 

in  the  same  manner  as  the  ionic  strength  is  varied,  then 
what  we  may  be  observing  is  an  apparent  change  in  q  for 
intercalative  binding  when,  in  reality,  the  apparent  change 
is  due  to  an  overlapping  of  surface  binding  manifestations 
with  intercalative  binding.   Assume  that  the  values  of  m  and 
q  remain  unchanged  for  each  type  of  binding,  then,  the  real 
slopes  of  the  two  line  segments  will  remain  constant  and 
their  point  of  intersection  will  be  dependent  on  the  ratio 
of  K  to  K  .   For  example,  in  one  extreme  case  there  is 
essentially  no  intercalative  binding,  K  /K   is  very  large, 
and  only  the  line  having  slope  1  will  be  observed.   At  the 
other  extreme,  K  /K  would  be  very  small  when  only  inter- 
calative binding  occurs  and  a  single  line  having  slope  2 
would  be  seen.   If  the  system  is  perturbed  in  a  manner  which 
caused  K   to  decrease  more  rapidly  than  KT  the  point  of 
intersection  of  the  two  lines  will  reflect  the  change  in 
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K  /K  .   In  plots  such  as  Figures  7  through  10,  this  would 

be  manifested  as  a  shift  toward  the  right,  or  to  less 

negative  log ( [BHP] / [BH] )  values.   The  total  number  of 

sites  which  are  involved  in  surface  and  intercalative 

binding  are,  of  course,  directly  related  to  the  magnitudes 

of  the  equilibrium  constants.   In  other  words,  an  increase 

in  K  /K   indicates  that  a  greater  percentage  of  the  total 

occupied  sites  are  occupied  by  surface-bound  species.   The 

ratio  of  K„  to  KT  is  seen  to  increase  for  the  binding 
si 

reactions  between  3-aminoacridinium  and  DNA  in  the  presence 

of  alkali  metal  electrolytes  as  the  ionic  strength  is  raised-. 

5 
The  following  data  are  representative  of  Kg/K  x  10   at 

ionic  strengths  of  0.0025,  0.0050,  0.025,  0.050,  and  0.10(M), 

respectively.   With  (CHO.NH^PO.  as  supporting  electrolyte: 

0.44,  0.56,  1.28,  1.83,  and  2.9;  whereas  with  LiH2P04  as 

supporting  electrolyte;  0.43,  10.7,  1.2,  5.0,  and  10.7, 

-5 
respectively.   Excepting  the  value  of  10.7  x  10    for  0.005 

LiH2P04,the  ratios  show  that  as  the  ionic  strength  is 

increased,  the  total  percentage  of  bound  drug  which  is 

surface-bound  increases.   Moreover,  for  any  given  ionic 

strength,  the  ratio  is  greater  when  tetrame thy 1 ammonium 

dihydrogen  phosphate  is  the  electrolyte  than  when  the 

lithium  salt  is  present.   It  is  obvious  from  the  decreases 

in  both  K   and  KT  that  increasing  the  dielectric  of  the 
si 

medium  does  not  enhance  surface  binding  relative  to  inter- 
calative binding.  Rather,  K  does  not  decrease  as  rapidly 
as  K  .   The  change  in  K  /K  may  result  from  specific  ion 
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competition  of  the  electrolyte  cations  with  the  drug 
cations.   As  the  ionic  strength  increases,  the  concen- 
tration of  cations  which  may  directly  compete  for  the  DNA 
phosphate  ligands  also  increases.   Moreover,  an  increase 
in  the  effectiveness  of  competition  by  the  metals  for  the 
ligands  would  be  expected  as  the  charge  density  on  the 
metal  increases.   Therefore,  specific  ion  interactions  may 
strongly  affect  surface  binding.   It  is  important  to  note 
that  no  evidence  from  the  experiments  conducted  with  the 
alkali  and  alkaline  earth  metals,  as  counterions,  indicated 
any  alterations  in  surface  binding  stoichiometry  resulting 
from  changes  in  the  medium  (aside  from  the  magnitude  of  K  ) . 
That  is,  all  initial  slopes  of  figures  such  as  Figure  7  are 
equal  to  1.   Apparent  changes  in  q  for  the  intercalative 
region  are  not  likely  to  be  due  to  intercalation  of  metal 
cations  because  intercalation  of  the  symmetrically  charged, 
highly  hydrophilic  metal  cations  is  improbable.   Specific 
ion  interferences  occurring  at  the  surface  may,  however, 
affect  the  overall  intercalation  of  the  drug  cations. 

Apparent  Association  Constants  for 
Surface  and  Intercalative  Binding 

Tables  24  through  32  include  the  apparent  association 
constants  for  surface  binding,  K  ,  of  3-aminoacridinium 
cation  to  DNA  in  the  presence  of  various  concentrations  of 
background  electrolytes.   For  a  given  titration,  values  of 

■ 

K  were  computed  at  each  data  point  as  discussed  in  the 
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Experimental  section.   By  definition,  an  equilibrium 
constant  is  independent  of  relative  and  absolute  con- 
centrations (specifically,  activities)  of  reactants  and 
products.   In  fact,  this  requirement  may  be  used  as  a 
test  of  the  validity  of  an  expression  postulated  to 
represent  an  equilibrium  situation.   Inspection  of 
representative  data,  included  in  Tables  18  through  21, 
reveals  that  the  values  of  K  for  points  within  a  specific 
titration  are  reasonably  constant  throughout  the  initial 
regions.   However,  in  many  cases,  notably  when  low  con- 
centrations of  alkali  metal  salts  are  present,  the  values 
calculated  from  the  later  points  of  the  titrations  are 
significantly  larger  than  those  calculated  from  the  initial 
portions.   For  example,  in  the  presence  of  0.0025M  CsH-PO., 

K  computed  using  points  from  the  initial  portion  of  the 

5-1  . 

titration  is  about  2  x  10  M   while  toward  the  end  it  is 

5   -1 
about  7  x  10   M   .   It  may  be  recalled,  however,  that  the 

plots  of  log  ([P,/m]-  q[BHP])  for  this  system,  indicate 

that  the  value  of  q  is  different  for  the  two  regions  and, 

therefore,  the  binding  reactions  are  different.   Thus,  when 

testing  various  values  of  m  and  q  in  search  of  the  pair 

which  results  in  the  most  consistent  K  ,  it  is  necessary 

that  K  be  constant  only  throughout  the  region  over  which 

it  applies.   Conversely,  upon  establishing  that  specific 

values  of  m  and  q  are  optimal,  the  test  of  constancy  may 

be  used  to  define  titration  regions  involving  only  one 

binding  mode.   The  region  of  the  titration  from  which 
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averaged,  reported  K  values  were  computed,  was  the  same 
for  all  titrations.   It  has  been  shown  that,  as  the  ionic 
strength,  I,  is  decreased,  the  ratio  of  bound  to  free 
drug  at  which  intercalative  binding  becomes  significant 
increases.   The  ratio  also  increases  as  the  charge  density 
of  the  countercation  increases.   For  the  most  dilute 
solutions  of  rubidium,  cesium,  and  tetrame thy 1 ammonium 
dihydrogen  phosphate,  intercalative  effects  begin  to 
become  important  when  the  logarithms  of  the  ratios  approach 
+0.2.   In  solutions  having  high  ionic  strengths,  or  which 
contain  divalent  cations,  surface  binding  predominates  even 
at  the  highest  ratios  of  bound  to  free  drug.   To  allow  re- 
liable comparisons  of  calculated  surface  binding  association 
constants  between  high  and  low  ionic  strength  solutions,  only 
data  where  log  ([BHP]/[BH])  was  less  than  +0.2  were  used. 
The  lower  limit  of  the  ratio  which  provides  for  reliable 

values  of  K  was  estimated  to  be  -0.8.   Below  -0.8,  data 
s 

often  deviated  significantly  from  values  obtained  above 

this  limit.   Using  the  above  guidelines,  the  averaged  values 

of  K   listed  in  Tables  24-32  were  computed.   It  is  readily 

acknowledged  that  the  method  of  limiting  the  range  of  data 

is  somewhat  arbitrary  and  may  be  eventually  replaced  with 

more  specific,  but  as  yet  unknown,  guidelines. 

The  number  of  individual  data  points,  n,  used  to 

calculate  K  ,  for  each  titration,  is  given  in  the  tables. 
s  3 

Specific  points  which  were  used  are  indicated  by  double 
asterisks  (**)  in  Tables  2  through  12  and  18  through  21. 
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In  most  instances,  n  was  greater  than,  or  equal  to  3. 
Standard  deviations  from  the  mean,  a,  are  also  included. 
Deviations  were  generally  within  5%  to  10%  of  the  averaged 
values  of  K  ,  which  shows  that  the  constancy  over  the 
intervals  selected  is  within  acceptable  limits  for  these 
types  of  systems.   As  a  rule,  the  lower  the  concentration 
of  background  electrolyte,  the  greater  the  standard  devia- 
tion.  Similarly,  the  lower  the  charge  density  of  the 
cation  of  the  supporting  electrolyte,  the  greater  the 
standard  deviation.   For  divalent  metal  electrolytes  the 
values  of  K  were  constant  over  nearly  the  entire  accessible 

range.   In  0.010N  Mg(02CCH3)  ,  K  varied  from  1.32  x  104  M_1 

4   -1 
to  1.21  x  10   M   over  the  entire  titration  for  an  average 

of  1.27  x  104  M_1  +  5%.   In  6.3  x  10~4N  Mg (02CCH  )  ,  K 

4-1  4 

went  from  a  low  of  4.71  x  10   M    to  a  high  of  5.09  x  10 

-1  4   -1 

M  "  (average:  4.90  x  10   M   +  4%).   The  essentially  in- 
variant apparent  equilibrium  association  constants  through- 
out the  entire  titrations,  along  with  the  single  straight 
line  plots  of  log ( [BHP] / [BH] )  vs.  log ( [P ./m]-  q[BHP]),  is 
convincing  evidence  that  only  the  effects  of  the  surface 
bound  species  are  observed  in  these  systems  and  that, 
either  intercalative  binding  is  not  taking  place,  or  it 

is  being  masked.   On  the  other  hand,  in  0.15M  CsH„PO,,  K 

2   4    s 

varies  from  an  average  over  the  whole  titration  by  +  25% 
and  by  +  55%  in  0.0025M  CsH"2P04 .   The  magnitudes  of  the 
standard  deviations  over  the  range  -0.8  <_  log  (  [BHP]/ [BH]  < 
0.2  are  proportional  to  those  given  above  for  the  entire 
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range.   The  higher  deviations  at  lower  ionic  strengths 
may  be  due  to  the  greater  extent  of  intercalative  binding 
during  the  initial  parts  of  the  titrations  at  the  low 
ionic  strengths.   Also,  there  may  be  a  change  in  the 
nature  of  the  surface  binding  arising  from  mild  denatura- 
tion  of  the  DNA.   In  solutions  of  low  ionic  strength,  the 
configuration  of  the  polyanion  may  be  sufficiently  altered 
from  that  of  its  native  state  to  affect  the  surface 
binding.   Addition  of  even  small  quantities  of  DNA  to 
solutions  having  very  low  initial  ionicities  causes  large 
relative  changes  in  their  ionic  strengths.   Therefore,  if 
DNA  configuration  is  dependent  on  ionic  strength  and  K  is 
dependent  on  DNA  configuration,  a  substantial  variation 
may  be  seen  over  a  narrow  titration  region. 

Agreement  between  apparent  surface  binding  association 
constants  obtained  from  replicate  titrations  was  generally 
within  +  10%  in  the  presence  of  all  concentrations  of 
divalent  metals  and  in  solutions  containing  monovalent 
cations  in  excess  of  about  0.010M.   In  light  of  the  com- 
plexity of  the  systems,  errors  of  this  magnitude  are 
considered  wholly  acceptable.   For  the  more  dilute  solutions 
of  monovalent  electrolytes,  agreement  between  replicate  runs 
was  considerably  poorer,  ranging  from  +  35%  to  +  3%.   It  is 
believed  that  the  lack  of  precision,  under  these  conditions, 
stems  from  uncontrolled  perturbations  of  the  configuration  of 
DNA  resulting  from  the  low  ionic  strength.   Changes  in  solu- 
tion pH  may  be  postulated  as  an  aggravating  factor  contributing 
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to  imprecision  in  the  poorly  buffered,  dilute  solutions 
in  that  Schulman  and  Capomacchia  (85)  have  shown  that 
pH  changes  from  5.9  to  7.4  may  substantially  alter 
binding  of  3-aminoacridinium  to  DNA.   However,  this  is 
not  likely  a  major  contributor  because  of  the  relatively 
small  amount  of  DNA  polyanion  added  to  the  solutions. 
Even  if  DNA  were  a  strong  pH  7  buffer,  the  amount  added 
would  not  affect  the  pH  as  much  as  the  far  greater  volumes 
used  for  titrations  in  the  presence  of  dilute  alkaline 
earths  where  good  precision  was  obtained.   Also,  at  pH 
5.9,  the  H„P04   anion  is  a  far  better  buffer  at  any  given 
concentration  than  is  the  acetate  anion  used  with  the 
alkaline  earth  electrolytes. 

The  limiting  Debye-Hiickel  law  predicts  a  linear 
relationship  between  the  logarithms  of  equilibrium  associa- 
tion constants  and  the  square  roots  of  the  total  ionic 
strengths  of  the  solutions.   Figures  11  through  16  are 
plots  of  the  log  of  the  apparent  surface  binding  associa- 
tion constant,  K  ,  vs.  the  square  root  of  the  total  ionic 

1/2 
strength,  I    .   Calculations  of  K  and  I,  from  which  the 

data  for  these  plots  are  derived,  were  based  on  simple 

molar  concentrations  of  all  species  in  solution.   We  have 

included  this  data  here  so  that  it  may  be  compared  with 

literature  data  and  also  to  allow  comparisons  with  more 

refined  results  presented  later.   The  values  of  K  are  the 

s 

averaged  terms  reported  in  Tables  24  through  32  derived 

1/2 
from  molar  concentrations.   Similarly,  the  I  '       values  are 
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averages  of  the  ionic  strength  at  the  beginning  and  at 

the  end  of  the  regions  over  which  K  was  obtained  and 

s 

include  contributions  of  all  species  toward  the  total 
ionic  strength  (all  concentrations  in  molar  terms) .   It 
is  immediately  apparent,  upon  inspection  of  the  figures, 
that  the  systems  do  not  obey  the  limiting  Debye-Hiickel 
relationship.   Most  significantly,  all  plots  exhibit 
decided  curvature  over  the  ionic  range  investigated.   Such 
non-linearity  is  not,  however,  particularly  surprising, 
nor  is  it  necessarily  unique  to  the  3-aminoacridinium-DNA 
system  under  consideration  since  the  limiting  case  is 
designed  only  for  dilute,  ideal  solutions  of  electrolytes. 
Ideal  systems  are  defined  as  being  infinitely  dilute  solu- 
tions consisting  of  point  charges  -  ions  having  no  dimen- 
sions and  engaging  in  none  other  than  simple  electrostatic 
interactions.   Real  solutions  which  approximate  the  ideal 
generally  consist  of  simple  inorganic  ions  having  high 
charge  densities  with  total  ionic  strengths  less  than  0.01M. 
Even  if  our  systems  were  ideal  in  all  respects  other  than 
concentration,  linearity  would  be  unlikely  over  the  broad 
ionic  strength  range  studied.   The  curvature  should,  if 
simply  a  result  of  deviations  from  electrostatic  ideality, 
be  obviated  using  corrections  included  in  extended  Debye- 
Hiickel  equations.   These  will  be  considered  below. 

Reconsidering  portions  of  Figures  11  through  16  in 
which  the  total  ionic  strength  is  less  than  0.01  (I1/2  <  0.1), 
we  notice  that  there  is  approximate  linearity,  most  easily 
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recognized  in  the  cases  of  divalent  metal  electrolytes. 
Precision  of  these  data  may  be  estimated  from  their  repro- 
ducibility at  a  given  ionic  strength  or  from  the  tendency 
of  data,  at  various  ionic  strengths,  to  form  a  smooth 
curve  when  plotted  against  I1'2 .   Using  either  of  these 
criteria,  it  appears  that  the  use  of  alkaline  earth 
acetates  as  background  electrolytes  provides  more  repro- 
ducible results  than  alkali  phosphates.   This  suggests 
that  pH  control  may  not  be  the  most  critical  factor  in 
maintaining  reproducibility  because,  for  the  same  ionic 
strength,  the  buffer  capacity  of  the  phosphates  is  markedly 
greater  than  that  of  the  acetates  used  for  the  alkaline 
earth  solutions. 

If  the  conformation  of  the  polyanion  is  an  important 
factor  in  surface  binding,  then  perhaps  it  is  maintained 
more  effectively  by  the  divalent  cations.   Their  greater 
charge  and  charge  density  may  provide  for  increased  inter- 
action between  the  polyanion  and  the  metal  cations  thus 
generating  a  disproportionately  higher  local  ionic  field 
around  the  polymer.   Hence,  though  the  bulk  ionicities  of 
the  two  types  of  solutions  are  similar,  the  local  electro- 
static potential  surrounding  the  polyanion  binding  sites 
may  be  much  greater  in  the  presence  of  the  alkaline  earths, 
thereby  stabilizing  the  native  DNA.   If  stabilization  in 
the  divalent  metal  acetate  systems  were  due  to  acetate 
rather  than  phosphate,  then  the  data  for  K02CCH  as  electro- 
lyte would  be  expected  to  be  more  precise  than  that  for 
KH2P04.   Figure  12  shows  that  such  is  not  the  case. 
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The  imprecision  at  low  ionic  strengths  is  sufficiently 
great  as  to  preclude  determining  straight  lines  having 
slopes  agreeing  within,  say,  +  20%  for  replicate  runs. 
The  large  magnitudes  of  the  surface  binding  constants 
for  our  model  compound,  to  DNA,  allow  nearly  optimal 
conditions  of  reactant  concentrations,  compared  to  many 
systems  reported  in  the  literature.   The  lack  of  precision 
evidenced  at  ionic  strengths  less  than  0.010N  for  our 
investigations  indicates  that  conclusions,  based  on  data 
obtained  at  only  a  single  ionic  strength  below  0.01N,  may 
be  suspect. 

Figures  11  through  16  reveal  a  second  manner  in  which 
these  systems  disobey  the  limiting  Debye-Hiickel  law  in 
that,  for  any  given  ionic  strength,  the  magnitude  of  Kg 
is  a  function  of  the  type  of  background  electrolyte.   Again, 
the  actual  case  is  not  surprising  since  we  are  dealing  with 
real  ions  having  finite  and  differing  charge  densities. 
These  differences  mean  that  the  electrostatic  fields  around 
the  ions  will  be  dissimilar  and,  in  turn,  the  dielectric 

properties  of  the  media  will  be  affected.   Extrapolation 

1/2 
of  plots  of  log  K  vs.  I  '       for  dilute  solutions  of  the 

several  different  types  of  electrolytes  would,  in  theory, 

yield  lines  intersecting  at  the  ordinate.   The  y-intercept 

would  correspond  to  the  log  of  the  association  constant  at 

zero  ionic  strength.   To  extrapolate  the  lines,  in  our  case, 

is  not  feasible  due  to  the  imprecision  of  the  data.   However, 

instructive  qualitative  information  may  be  obrained  by 
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comparing  values  of  K   for  solutions  of  intermediate  ionic 
strength.   For  0.010M  solutions  of  the  following  cations, 
the  apparent  association  constant  for  the  reaction  between 
3-aminoacridinium  cation  and  DNA  decreases  in  the  order: 
(CH3)4N+>Rb+>K+>Na+> (Cs+)>Li+.   The  order  is  the  same  in 
0.025M  solutions  (Rb  was  not  evaluated  at  this  concentra- 
tion) and  also  in  0.050M  MH2PO..   In  0.10M  solutions,  the 
order  was  (CH3) 4N+>> (Cs+) >K+   Rb+>Na+>Li+.   In  0.010M 
electrolytes  Kg  ranged  from  17  x  104  M-1  to  9  x  104  M_1, 
in  0.025M  electrolytes,  from  11  x  104  M_1  to  5.4  x  104  M_1, 
in  0.050M  electrolytes,  from  7.3  x  104  M-1  to  3.9  x  104  M-1, 
and  in  0.10M  electrolytes,  from  6.1  x  104  M_1  to  2.9  x  104  M-1 
With  the  exception  of  Cs  ,  the  order  of  the  association 
constants  parallels  the  order  of  the  sizes  of  the  electro- 
lyte cations.   Ionic  radii  for  the  cations  are  (86):  Cs+(1.68 

A),  (CH3)4N+(>1.3  A),  Rb+(1.48  A),  K+(1.33  A),  Na+(0.97  A), 

+      o  , 

and  Li  (0.64  A).   The  exact  size  of  the  hydrated  (CH-J  N 

ion  is  not  known  but  it  is  likely  that  it  is  at  least  as 
large  as  Rb  .   It  has  no  available  empty  orbitals  through 
which  any  form  of  dative  binding  may  occur  and  thereby, 
provides  a  means  of  assessing  possible  effects  of  available 
unfilled  metal  orbitals  on  the  metal's  role  in  the  inter- 
action between  the  drug  and  DNA.   Qualitatively,  the 
behavior  of  the  system  in  the  presence  of  the  totally  satu- 
rated tetramethylammonium  cation  is  the  same  as  in  the 
presence  of  any  of  the  alkali  metal  cations.   Quantitatively, 
the  apparent  surface  binding  association  constants  are 


94 


greatest  for  the  reaction  in  (CH-.)  4N   solutions,  in  accord 
with  its  large  size  and  low  charge  density.   There  is  no 
evidence  that  the  empty  orbitals  impart  any  special 
properties  to  the  alkali  metals  in  terms  of  their  effects 
on  the  binding  of  3-AA  to  DNA. 

The  seemingly  strong  suppression  of  the  binding  process 
by  cesium  ion  is  not  logical  in  view  of  the  effects  of 
(CH-J.N   and  Rb  .   Because  Cs   is  larger  than  Rb   it  would 
be  expected  to  exert  a  weaker  or,  at  least,  similar  effect 
than  rubidium.   Despite  the  fact  that  Cs   is  in  the  ex- 
pected order  in  the  0.10M  series,  the  otherwise  anomalous 
behavior,  in  its  presence,  is  likely  due  to  some  undefined 
variable  such  as  its  original  purity.   The  20%,  by  weight, 
impurity  of  the  solid  CsOH  used  to  prepare  the  phosphate, 
as  determined  by  acid-base  titration,  was  considered  to  be 
entirely  water  of  hydration. 

At  any  given  ionic  strength,  based  on  formal  molarities, 
the  alkaline  earths  more  effectively  suppress  the  binding 
of  3-aminoacridinium  to  DNA  than  do  the  alkali  metals. 
This  behavior  is  not,  of  course,  anticipated  on  the  basis 
of  the  limiting  Debye-Huckel  law  but  may  be  partially 
rationalized  with  extended  equations.   The  order  of  the 
magnitudes  of  log  Kg  for  the  binding  of  3-AA  to  DNA  in 

the  presence  of  the  alkaline  earths,  decreases  in  the 

++   ++   ++   ++ 
order:  Ba   >Sr   >Mg   >Ca    over  the  entire  range  of  ionic 

strength  studied  (for  0.010N  electrolyte,  log  K   is  slightly 

less  with  Mg   than  with  Ca   ) .   In  decreasing  size  of 
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++       °     ++ 
unhydrated  cations  (86)  the  order  is:  Ba   (1.35  A) >Sr 

(1.13  A)>Ca++(0.96  A) >Mg++ ( . 65  A).   Ion  size  parameters 

for  the  metal  ions  are  5,  5,  6,  and  8,  respectively  (87). 

The  order:  Ba++>Sr++>Ca  +  is  logical,  both  in  terms  of 

ion  size  and  in  terms  of  ion  size  parameters.   However, 

calcium  more  strongly  inhibiting  the  binding  of  3-AA  to 

DNA  than  does  Mg++  is  in  opposition  to  charge  density  and 

ion  size  arguments.   This  reversed  order  indicates  that 

the  effect  of  the  metal  ions  encompasses  more  than  simple 

bulk  solution  electrostatic  factors  which  can  be  correlated 

with  ion  sizes  and  charge  densities.   The  most  apparent 

additional  effect  is  specific  ion  competition  in  which 

the  overall  energetics  of  the  competition  for  phosphate 

sites  must  be  evaluated. 

Figures  17  through  2  0  are  plots  of  the  apparent 

intercalative  association  constants,  K  ,  against  the  square 

1/2 
root  of  the  ionic  strength,  I  '  ,  for  the  alkali  metals. 

They  are  similar  to  Figures  11  through  16  as  they  show  the 

1/2 
relationship  between  log  K  and  I  '      based  on  conventionally 

employed,  uncorrected  molar  concentrations.  In  the  presence 

of  the  alkali  metals,  the  relationships  between  the  apparent 

intercalative  association  constants,  K  ,  and  the  ionic 

strength  were  essentially  the  same  as  for  K  and  I  discussed 

s 

above.    That  is,  limiting  Debye-Hiickel  behavior  is  not 
obtained  over  the  entire  concentration  region.   Curvature 
is  especially  apparent  at  higher  concentrations,  and  the 
nature  of  the  background  cation  affects  the  degree  of 
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binding.   Except  in  the  presence  of  0.10M  electrolytes, 
the  order  of  the  magnitudes  of  the  intercalative  binding 
constants  for  3-aminoacridinium  and  DNA  in  the  presence 
of  the  monovalent  cations  is  the  same  as  for  surface 
binding:  (CH3)  N+>Rb+>K+>Na+> (Cs+) >Li  .   Again,  the 
sequence  can  be  rationalized  on  the  basis  of  charge 
density,  excepting  cesium.   In  0.10M  solutions,  the 
order  is  (CH3) .N+>>Rb+  - K+> (Cs+) >Na+>Li+.   This  series 
differs  from  the  one  for  K  only  in  the  position  of  Cs 
which,  again,  is  probably  an  experimental  artifact.   We 
were  unable  to  observe  intercalative  binding  of  3-amino- 
acridinium to  DNA  in  the  presence  of  alkaline  earth  cations 
at  any  concentrations.   Data  were  generated  and  plotted 
for  these  systems  similar  to  Figures  17  through  20  but 
are  not  included  here.   The  order  of  the  "K  "  values  so 
determined  were  entirely  random  with  respect  to  either 
unhydrated  ionic  radii  or  ion  size  parameters.   Possible 
reasons  why  the  intercalative  binding  was  not  seen  in 
these  media  were  presented  previously. 

To  summarize,  the  relationship  between  apparent 
association  constants  for  both  surface  and  intercalative 
binding  of  3-aminoacridinium  to  DNA  over  the  concentration 
range  of  0.10N  to  0.00063M  alkali  metal  background  electro- 
lyte and  0.025N  to  0.00063N  alkaline  earth  electrolytes  do 
not  obey  the  Debye-Hiickel  limiting  law.   In  all  cases,  the 
plots  are  non  linear  and  the  magnitudes  of  the  association 
constants  are  dependent  on  the  nature  of  the  cations.   With 
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the  exception  of  Cs  ,  the  relative  effects  of  the  alkali 
metals  parallel  their  respective  charge  densities.   The 
same  is  true  for  surface  binding  in  the  presence  of  di- 
valent cations,  except  that  Ca    and  Mg++  are  reversed. 
Plausible  comparisons  between  charge  density  and  relative 
values  of  K   in  the  presence  of  divalent  electrolytes 
cannot  be  made  because  of  the  inability  to  evaluate  inter- 
calative  binding  in  solutions  containing  the  divalent 
cations.   The  implications  by  other  investigators  that 
there  is  a  straight-forward  relationship  between  ionic 
composition  of  solutions  in  terms  of  molarities  and  binding 
constants  must  be  considered  with  reservations. 

Association  Constants  Based  on  Activities  — 
Extended  Debye-Hiickel  Considerations 

The  procedures  outlined  above  for  the  determination  of 
binding  constants  of  small  molecule  to  DNA  are  employed 
by  the  majority  of  investigators  in  this  field.   That  is, 
calculations  are  made  with  the  assumption,  either  explicit 
or  implied,  that  molar  (or  molal)  concentrations  of  species 
are  equal  to  their  activities.   Only  when  activity  coeffi- 
cients are  equal  to  1  is  this  assumption  strictly  valid  — 
a  condition  met  only  by  dilute  solutions  of  "well-behaved" 
systems  such  as  fully  ionized,  monovalent,  inorganic  salts. 
Ideally,  evaluation  of  the  relationship  between  ionic 
strength  and  association  constants  should  be  based  on 
solutions  having  total  ionic  strengths  approaching  zero,  and 


98 


including  low  concentrations  of  reactants  and  products. 
Unfortunately,  the  conformational  instability  of  DNA  in 
low  ionic  strength  media  (<  0.0050 (M) )  prevents  such 
investigations.   Furthermore,  limits  of  instrumental 
detectability  and  the  nature  of  the  equilibrium  demand 
moderately  high  concentrations  of  the  biopolymer.   But, 
while  concentrations  of  background  electrolyte  and 
reactants  must  be  higher  than  desired,  a  more  accurate 
reflection  of  the  relationship  between  ionic  strength 
and  the  association  of  small  molecules  to  DNA  may  be 
made  using  species'  activities  as  opposed  to  molar 
concentrations . 

Activities  of  each  electrolyte  at  the  various  molar 
concentrations  were  estimated  by  determining  the  mean 
activity  coefficients  of  the  constituent  ions  from 
equation  (2-18) .   Ion  size  parameters,  d,  were  obtained 
from  the  literature  (87)  and  are  included  in  Table  33. 
The  activities  of  the  electrolytes  were  used,  in  turn, 
in  equation  (2-18)  as  the  values  of  I  when  estimating  the 
activity  coefficients  of  the  species  involved  in  the  binding 
processes.   Ion  size  parameters  for  BHP,  BH,  and  DNA  phos- 
phate were  not  available  and  had  to  be  estimated.   DNA 
phosphates  were  assigned  a  d  value  of  4  x  10   ,  being 
approximately  the  same  as  the  value  for  free  H-PO  ~.   Free 

and  bound  3-aminoacridinium  cation  were  given  d  values  of 

—  8  —8 

3  x  10  "  and  4  x  10   ,  respectively.   The  assignment  of  ion 

size  parameters  for  these  species  is,  admittedly,  rather 
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arbitrary  as  there  are  few  guidelines  by  which  to  estimate 
them.   However,  it  is  the  effects  of  varying  the  back- 
ground electrolyte  which  is  of  principal  concern  here. 
So,  as  long  as  the  values  of  d  for  the  species  involved 
in  the  equilibrium  are  held  constant,  relative  effects 
of  changing  the  medium  can  be  estimated  though  absolute 
constants  may  be  incorrect.   After  calculating  activity 
coefficients  for  reactants  and  products,  the  correction 

factor,  E  #  was  computed  from  equation  (2-17)  where  E   = 
b  s 

^BHp/^BH*  9p^  "   Table  33  lists  the  correction  factors  for 
the  surface  binding  along  with  logarithms  of  the  corrected 
association  constants,  log  K  ,  in  the  presence  of  varying 
concentrations  of  the  electrolytes.   For  any  given  initial 
molarity  of  alkali  metal  electrolyte,  only  minor  variations 
in  Es  are  observed  for  the  differing  types  of  cations. 
These  small  differences  are  reflections  of  the  minor 
effect  d  has  on  the  correction  factor.   For  instance, 
values  for  0.10M  LiH2P04  and  0.10M  CsH2PO   are  4.28  and 
4.20,  respectively,  where  d  is  8  x  10~8  and  2.5  x  10_8  for 
the  respective  cations.   Below  about  0.010M  electrolyte 
the  correction  factor  is  independent  of  the  character  of 
the  monovalent  cation.   For  ionic  strengths  greater  than 
about  0.01(M),  the  order  of  increasing  log  K   is  the  same 
as  for  log  Kg ,  namely:  (CH3) 4N+>Rb+>K+>Na+> (Cs+) >Li+. 
(Cesium's  position  is  not  the  same,  but  data  pertaining 
to  it  are  probably  not  reliable,  vide  supra.)   The  parallel 
order  for  the  corrected  values  with  those  based  on  molar 
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concentrations  is  expected  in  light  of  the  relatively- 
minor  effect  that  the  ion  size  parameter  has  on  E  . 
Below  I  =  O.OIO(M),  the  data  is  too  scattered  to  make 
definitive  statements.   Too,  at  low  ionic  strengths,  the 
assumption  that  DNA  phosphates  can  be  treated  as  if  they 
were  species  independent  of  one  another  becomes  especially 
approximate.   Free  species  may  be  diluted  to  any  degree 
whereas,  below  a  given  value,  further  addition  of  bulk, 
solvent  will  have  little  effect  on  altering  the  ionic 
environment  of  a  site  on  the  polymeric  surface  since  its 
environment  is  predominantly  determined  by  near-neighbor 
charges  along  the  polymer  and  not  by  charged  species  in 
the  bulk  solvent.   Hence,  the  greater  the  degree  of  dilution, 
the  greater  the  disparity  between  the  calculated  ionicity 
and  the  true  ionic  strength  local  to  the  binding  site. 
There  is  a  substantial  decrease  in  the  range  of  K 


i 

s 


compared  to  K  values  for  the  range  of  concentration  (activity) 
studied.   At  higher  ionic  strengths  the  value  of  £   is  maximal 
while  K  is  minimal  while  at  low  ionic  strength  E   is  low 

with  log  K  being  at  its  maximum.   The  product  of  the  two, 

i 

therefore,  tends  to  decrease  the  range  of  K   relative  to  K  . 

1       1/2 
As  a  result,  plots  of  log  K  vs.  I    are  more  linear  as 

evidenced  in  Figures  21  and  22.   Having  corrected  for 

variations  in  activities  of  species  in  solution,  and  also 

having  presumably  eliminated  the  effects  of  differing  elec- 

trolyte  ion  sizes,  theory  predicts  that  plots  of  log  K  vs. 

1/2 
I    would  be  straight  lines  which  are  superimposable . 
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Obviously,  while  data  for  I  >  0.010(M)  for  most  of  the 

systems  do  form  reasonably  straight  lines,  they  do  not 

superimpose.   Moreover,  there  is  still  some  curvature  in 

the  plots  which  is  most  pronounced  for  lithium,  becoming 

less  apparent  as  one  goes  from  Na   to  (CH-.)  .N+.   These 

two  features  likely  result  from  specific-ion  competition 

by  metals  for  the  phosphate  sites  on  the  DNA. 

Figures  23  and  24  are  plots  of  the  corrected  surface 

binding  equilibrium  association  constants,  K  vs.  I  ' 

s 

for  reactions  in  the  presence  of  alkaline  earth  electro- 
lytes.  It  can  be  seen  from  the  figures  and  from  data  in 
Table  33  that  the  order  of  the  magnitudes  of  K   in  the 
presence  of  the  alkaline  earths  is  dependent  on  background 
electrolyte  concentration.   This  is  in  contrast  to  the 
trends  for  Kg  for  alkali  metal  solutions  and  also,  for 
the  uncorrected  values  of  both  alkali  and  alkaline  earth 
electrolytes.   The  orderings  are  not  random  but  are  due 
to  the  unequal  slopes  of  the  distinctly  straight  lines 
for  each  type  of  cation.   The  lack  of  superposition  of 
the  lines  shows  that  specific  ion  competition  probably  is 
operative. 

The  above  treatment  has  shown  that  application  of 
extended  Debye-Huckel  equations  to  surface  binding  data  is 
sufficient  to  rationalize  the  nonlinear  relationship  between 
Ks  and  I*   The  remaining  deviations  are  probably  due  to 
specific  ion  effects  which  will  be  briefly  considered  below. 
Similar  attempts  to  explain  the  intercalative  data  are  not 
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so  fruitful.   The  correction  factor  for  the  intercalative 
binding  expression,  £  ,  is  essentially  a  constant,  equal 
to  1,  for  all  ionic  strengths  investigated  (the  highest 

value,  for  I  =  0.15(M),  is  1.05).   Of  course,  the  wide 

1/2 
range  of  K  values  and  the  nonlinearity  of  log  K  vs .  I 

are  unaltered  by  application  of  the  correction  for  activi- 
ties.  The  implications  here  are  uncertain  at  this  time, 
but  it  appears  that,  whatever  it  is  that  predominantly 
affects  the  intercalative  binding,  it  is  not  explainable 
in  terms  of  simple  concentrations  or  ionicities.   One 
possibility  is  that  the  intercalative  mode  is  highly 
susceptible  to  minor  changes  in  the  pitch  of  the  DNA 
which  is,  in  turn,  a  function  of  ionic  strength. 

Evaluation  of  a  Simple  Competitive  Binding  Model 

A  thorough  investigation  of  specific  ion  competition 
between  electrolyte  metals  and  3-aminoacridinium  cations, 
for  DNA  phosphate  sites,  would  include  complex  equilibria 
involving  mixed  complex  systems  which  are  beyond  the  scope 
of  this  work.   Such  a  study  would  include  the  selection  of 
probable  stoichiometries  of  metal-DNA,  or  metal-drug-DNA 
mixed  complexes  using  procedures  similar  to  those  for  the 
drug-DNA  systems.   These  models  would  then  be  evaluated  to 
determine  the  one  most  probable,  on  the  basis  of  a  "best-fit." 
The  most  elementary  of  these  will  be  briefly  considered  here. 

This  system  would  be  one  in  which  there  are  no  mixed 
complexes  and  in  which  the  binding  stoichiometries  of  the 
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metal,  and  3-AA,  to  DNA  are  the  same.   It  has  been  shown 
that  3  DNA  phosphates,  acting  as  a  single  tridentate 
ligand,  bind  to  each  3-AA,  which  we  shall  assume  is  the 
same  for  the  metal  binding.   Considering  the  3  phosphates 
as  constituting  each  site,  P,  in  equations  (1-4)  and  (1-5) , 
we  may  apply  equation  (1-7)  to  our  system.   Then,  the 

relationship  between  the  observed  association  constant, 

■ 

K  ,  for  the  3-aminoacridinium-DNA  reaction,  the  actual 
s' 

it 

association  constant  for  drug-DNA  binding,  K  ,  and  the 
association  constant  for  the  metal-DNA  interaction,  K  , 
as  developed  by  Passero  and  Gabbay  (64)  would  be 

K  =  K  /(l  +  KM  M)  (3-5) 

S      S  M 

where  [M]  is  the  molar  concentration  of  the  electrolyte 

metal  cation.   Substituting  values  of  K  and  their  cor- 

3  s 

responding  [M]  values  into  equation  (3-5)  allows  for  the 

ii 

simultaneous  solution  of  K   and  K  .   Table  34  lists  calculated 

s      M 
ii 

values  of  log  K  and  K   for  the  competitive  reactions  of 
alkali  and  alkaline  earth  cations  and  3-AA  for  DNA  phosphate 
sites.   The  concentration  ranges  of  electrolytes  listed  in 
Table  34  were  selected  on  the  basis  of  the  straight  line 
regions  of  plots  of  K   against  [M] .   For  the  simple  systems 
assumed  here,  there  would  be  a  linear  relationship  between 
the  two  parameters  because  of  the  first  order  concentration 
terms  for  all  species  in  both  equilibrium  expressions.   Data 
for  lithium  as  an  electrolyte  were  not  included  in  Table  34 
since  there  was  no  region  of  linearity  throughout  the  whole 


104 


range  of  metal  ion  concentration.   This  may  indicate 
that  the  interaction  between  DNA  and  Li   is  substantially 
different  than  between  DNA  and  the  other  metals.   Caution 
must  be  advised  in  drawing  conclusions  about  trends  ob- 
served in  the  values  of  the  metal  association  constants 
for  each  series.   These  values  are  obtained  as  the  result 
of  small  differences  between  large  numbers  (K  )  which,  in 
themselves,  have  large  error  limits.   It  is  evident, 
though,  that  K  values  for  the  divalent  metals  are  approxi- 
mately an  order  of  magnitude  greater  than  the  corresponding 
values  for  the  alkali  metals.   The  logarithm  of  the  equilib- 
rium association  constant  for  the  surface  binding  of  3- 
aminoacridium  cation  to  DNA,  having  been  corrected  for 
ionic  strength,  ion  sizes  of  species  in  solution,  and 
specific  ion  interactions  by  alkali  metals  appears  to  be 
5.41  +  0.10.   While  not  as  large  as  for  the  alkali  metal 

"  II 

value  of  K  ,  the  value  of  4.71  +  0.08  for  K   in  the  presence 
s  —  s 

of  divalent  metal  competitors  appears  reasonable.   Having 

ostensibly  corrected  for  ionic  charges,  along  with  the 

■I 

above-mentioned  parameters,  the  values  of  log  K   for  the 

s 

alkali  and  alkaline  earth  series  should  be  equal.   That 
they  are  substantially  different  proves  that  the  model 
needs  further  refinement.   It  is  logical,  though,  that 
the  binding  stoichiometries  of  divalent  cations  and  mono- 
valent cations  with  DNA  phosphates  would  be  different. 
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Thermodynamics  of  the  Binding  of  3-aminoacridinium 
and  7-aminoquinolinium  to  DNA 

The  standard  free  energies,  AG  ,  enthalpies,  AH°, 
and  entropies,  AS  ,  for  the  two  binding  reactions  of  3- 
aminoacridinium  and  7-aminoquinolinium  to  DNA  were 
evaluated  from  the  temperature  dependencies  of  their 
equilibrium  association  constants.   By  comparing  the 
relative  magnitudes  of  AH°  and  AS°  for  3-aminoacridinium 
with  those  of  its  benolog,  7-aminoquinolinium,  specific 
effects  resulting  from  variations  in  the  aromatic  portion 
of  the  molecule  may  be  evaluated.   Differences  in  chemical 
properties,  especially  acidities,  must  be  borne  in  mind 
when  making  such  comparisons. 

The  general  spectral  characteristics  of  the  3-amino- 
acridinium-DNA  system,  along  with  the  regions  in  the  titra- 
tions corresponding  to  surface  and  intercalative  binding 
and  the  evaluation  of  their  constants  have  all  been 
previously  discussed  for  reactions  at  25.0°C.   Except 
for  small  changes  in  the  magnitudes  of  the  equilibrium 
constants,  the  systems  at  15.0°C  and  35.0°C  behave  in 
essentially  the  same  manner.   These  topics  will  be  briefly 
considered  for  the  reactions  between  7-aminoquinolinium  and 
DNA  at  25.0°C. 

Figure  25  contains  representative  spectra  of  the 
absorptiometric  titration  of  a  2  x  10  M  solution  of  7- 
aminoquinolinium  with  the  potassium  salt  of  DNA  (10~2M  in 
DNA  phosphate) .   The  total  absorbances  of  the  7-AMQ  solutions 
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are  not  as  great  as  for  the  3-AA  solutions  due  to  the 

former's  lower  molar  absorptivities  (free,  7940_  M   cm 

-1   -1 
at  392  nm;  bound,  5332  M   cm   at  407  nm) .   Drug  concen- 
trations could  have  been  used  which  would  have  allowed 
absorbances  greater  than  0.4.   However,  it  was  decided 
to  use  the  same  formal  concentrations  of  7-aminoquinoline 
as  3-aminoacridine.   Then,  any  differences  observed  in 
the  titration  behavior  of  the  two  could  not  be  due  to 
differences  in  their  concentrations.   Also,  the  lower 
binding  affinity  of  7-aminoquinolinium  for  DNA  makes  it 
desirable  to  keep  its  concentration  at  a  minimum  to  allow 
maximum  ratios  of  DNA  to  drug  without  introducing  compli- 
cating concentration  problems. 

The  visible  spectra  of  both  the  free  and  bound  7- 
aminoquinolinium  cation  consist  of  only  one  envelope, 
having  no  fine  structure,  unlike  the  spectra  of  the  3- 
aminoacridinium  species  which  have  well  differentiated 

L   and   L,  bands.   Presumably,  substitution  of  the  amino 
group  in  the  3  position  of  the  molecule  results  in  the   L, 

transition  having  increased  probability  over  the   L  .   Thus, 

a 

the  observed  spectra  consist  predominantly  of  the   L,       A 

transition.   The   L  band  may  be  buried  beneath  the   L. 

though  no  attempts  were  made  to  ascertain  its  position. 

The  lack  of  distinct   L   and   L  bands  and  the  absence  of 

a       b 
fine  structure  in  the  spectra  of  either  the  free  or  bound 

species  preclude  discussions  concerning  possible  molecular 

orientations  upon  binding.   There  is  no  evidence  to  indicate, 
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however,  that  7-aminoquinoliniura  does  not  interact  with 
DNA  in  essentially  the  same  manner  as  does  3-aminoacridinium. 
Similar  reaction  stoichiometries  and  values  of  m  and  q  in 
equation  (1-21)  are  in  agreement  with  this  presumption. 
Significant  shifts  in  the  baseline  of  the  spectra  occur 
with  the  addition  of  moderate  to  large  amounts  of  DNA. 
Figure  25B  is  a  reproduction  of  a  titration  spectrum  in 
which  it  may  be  seen  that  the  absorbance  at  wavelengths 
greater  than  475  nm  increases  by  approximately  0.03  absor- 
bance  units  upon  addition  of  1400  yl  of  10   M  DNA  phosphate 
to  8.04  ml  of  solution  (resultant  P/D  =  105).   Magnitudes 
of  baseline  shifts  for  7-aminoquinolinium  titrations  are 
particularly  significant  in  view  of  the  relatively  small 
total  changes  in  absorbance,  at  the  analytical  wavelength, 
throughout  the  titrations.   Reverse  titrations  in  which 
stock  DNA  solution  was  used  to  titrate  pH  5.90  buffer  in  a 
4.00  cm  cell  showed  that  an  increase  in  the  baseline  absor- 
bance of  about  0.03  units  may  be  anticipated  from  light 
scattering  arising  from  the  presence  of  relatively  high 
concentrations  of  DNA  polymer.   Correcting  the  absorbances 
at  500  nm  of  each  individual  spectra  to  zero  results  in 
absorptiometric  titration  spectra  similar  to  the  one  shown 
in  Figure  25A.   Here  it  becomes  evident  that  what  appeared 
to  be  an  isobestic  point  at  -  420  nm  was,  in  fact,  an  arti- 
fact and,  thus,  the  system  may  involve  more  than  two  species 
which  are  in  equilibrium.   It  is  also  evident  that,  even 
after  addition  of  1400  yl  of  DNA  solution  to  the  8.04  ml  of 
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drug,  (P/D  =  105)  the  titration  was  not  complete.   Adding 
greater  amounts  of  DNA  is  counterproductive  because  of 
increasingly  severe  light  scattering  and  the  dilution  of 
the  sample.   Baseline  absorbance  at  the  analytical  wave- 
length was  determined  at  the  point  of  intersection  of  a 
line  drawn  between  the  absorbances  at  475  nm  and  325  nm,  and 
a  perpendicular  at  the  analytical  wavelength.   The  analytical 
wavelength  was  that  of  the  band  maximum  for  the  total  solu- 
tion. 

The  band  maximum  for  free,  protonated  7-aminoquinolinium 
is  at  392  nm.   As  DNA  is  added  to  the  drug,  there  is  both  a 
decrease  in  the  intensity  of  the  band  and  a  red  shift  of  its 
maximum.   The  hypochromism  and  bathochromism  are  both  quali- 
tatively similar  to  those  seen  in  the  titrations  of  3-amino- 
acridinium.   However,  the  7-AMQ  red  shift  indicates  a 
significantly  greater  energy  difference  between  the  free 
and  bound  7-aminoquinolinium  compared  to  the  difference 
between  free  and  bound  3-aminoacridinium.   In  the  presence 
of  a  100-fold  excess  of  DNA  phosphate  to  drug,  the  band 
maximum  for  bound  7-aminoquinolinium  appears  at  407  nm  cor- 
responding to  a  difference  in  the  transition  energies  between 
the  ground  and  excited  states  of  the  free  vs.  the  bound 
complexes  of  940  cm"   (2.70  Kcal/mole  7-AMQ).   Recall  that 
the  transition  energy  difference  for  free  and  bound  3-amino- 
acridinium amounted  to  only  about  1.3  Kcal/mole  3-AA. 

There  were  no  titrations  of  7-aminoquinolinium  in  which 
an  abrupt  shift  in  the  spectra,  similar  to  that  seen  in  the 
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L,  band  of  the  3-aminoacridinium  system,  was  observed. 
The  cbsc-nce  of  a  dramatic  shift  at  DNA  phosphate  levels 
in  excess  of  those  causing  shifts  in  the  3-aminoacridinium 
systems  shows  that  shifts  in  the  latter  instances  were  due 
to  something  other  than  simple  changes  in  the  bulk  solvent's 
composition.   It  is  not  known  if  greater  concentrations  of 
DNA  would  have  resulted  in  dramatic  shifts  in  the  7-amino- 
quinolinium  spectra. 

Figures  26  and  27  are  plots  of  log ( [ART]/ [AH] )  vs. 
log( [Pt/3]-[AHP] )  and  log ( [P  /2] -2 [AHP] ) ,  respectively  for 
the  titrations  of  7-aminoquinolinium  with  DNA  at  15.0  C, 
25.0  C,  and  35.0  C.   In  all  instances,  two  distinct  regions 
of  linearity  are  seen,  indicative  of  two  separate  modes  of 
binding.   In  Figure  26  the  averaged  slopes  of  the  initial 
portions  of  the  curves,  corresponding  to  the  surface  binding 
process,  is  1.03  +  0.13.   This  value  is  nearly  the  same  as 
for  3-aminoacridinium  and  shows  that  equation  (2-15)  also 
applies  to  7-aminoquinolinium  surface  binding.   The  averaged 
values  for  the  same  region  of  Figure  27  is  1.02  +  0.12.   The 
value  for  q  for  the  intercalative  binding,  obtained  from  the 
slopes  of  the  later  regions  of  the  lines  in  Figure  27,  is 
1.64  +  0.09  (from  Figure  26,  q  =  1.8  +  0.4).   These  are 
somewhat  ambiguous  in  terms  of  the  actual  integral  value  of 
q  but  there  is  little  doubt  that  it  is  greater  than  1  and 
is  presumed  to  be  less  than  3,  based  on  the  lack  of  constancy 
of  the  computed  values  of  K,  when  q  =  3.   Because  of  the 
relatively  weak  binding  of  7-aminoquinolinium  to  DNA,  it  was 
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not  possible  to  titrate  the  system  in  a  region  in  which 
the  surface  binding  did  not  overlap  intercalative  binding. 
As  a  result,  the  later  points  in  Figures  26  and  27  probably 
include  some  contributions  from  surface  binding  which  would 
cause  the  slopes  to  be  less  than  2.   On  the  bases  of  best- 
fit  analyses  and  the  similarity  of  the  binding  with  that  of 
3-aminoacridinium,  values  of  m  =  q  =  2  were  selected  for 
the  intercalative  mode  of  7-aminoquinolinium  binding  to 
DNA. 

Table  35  includes  the  calculated  apparent  association 
constants  for  the  surface  and  intercalative  binding  of  3- 
aminoacridinium  and  7-aminoquinolinium  to  DNA  at  15.0  C, 
25.0°C,  and  35.0°C.   Constants  were  computed  using  the 

estimated  activities  of  the  species  (equations  (1-30)  and 

—  8 
(2-17)).   An  ion  size  parameter  equal  to  3.5  x  10    was 

used  for  the  free  drug  cations  and  bound  complexes.   For 
any  given  temperature,  the  apparent  surface  binding  associa- 
tion constant  of  7-aminoquinolinium  is  approximately  an  order 
of  magnitude  lower  than  for  3-aminoacridinium  while  the 
constants  for  the  intercalative  process  are  about  two  orders 
of  magnitude  lower  for  the  former.   Magnitudes  of  Kg  and  K1 
for  both  compounds  decrease  with  increasing  temperature  in- 
dicating a  negative  free  energy  for  the  binding  process. 
Figures  28  and  29  are  plots  of  log  K  vs.  1/T  (  K)  for  3- 
aminoacridinium  and  7-aminoquinolinium,  respectively,  and 
Figures  30  and  31  are  similar  plots  of  log  K  vs.  1/T  (  K) . 
With  the  possible  exception  of  Figure  31,  the  relationship 
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between  log  K  and  1/T  is  linear  as  predicted  by  the 
van't  Hoff  relationship  (equation  (1-31)).   Values  of 
the  standard  enthalpy  for  the  reactions,  included  in  Table 
36,  were  calculated  from  the  slopes  of  the  lines  of  Figures 
28  through  31.   Standard  entropies  of  binding  were  computed 
from  values  of  AH   and  AG   at  298  K.   The  standard  state 
upon  which  the  thermodynamic  parameters  are  based  is 
related  to  activities  of  the  species  estimated  from  their 
molar  concentrations  using  equation  (1-30) .   While  the 
total  free  energy  of  binding  of  3-aminoacridinium  is  more 
negative  than  for  7-aminoquinolinium  (AG  =  -7.35  Kcal/mole 
3-AA  vs.  -6.07  Kcal/mole  7-AMQ)  the  standard  enthalpy  change 
or  heat  of  reaction  involving  3-aminoacridinium  is  less  than 
that  involving  7-aminoquinolinium  (AH   =  -4.90  Kcal/mole  3-AA 
vs.  -5.13  Kcal/mole  7-AMQ).   Enthalpy  is  predominantly  a 
measure  of  the  overall  energy  changes  resulting  from  the 
disruption  of  bonds  associated  with  reactants  and  the 
formation  of  new  bonds  among  product  species.   Considered 
in  this  vein,  it  is  evident  that  the  greater  binding  affinity 
of  3-aminoacridinium  over  7-aminoquinolinium  for  DNA  cannot 
be  attributed  to  electrostatic  bonding  energies  because, 
based  solely  on  AH°,  the  binding  of  7-AMQ  to  DNA  would  be 
greater  than  that  of  3-AA.   Two  types  of  electrostatic  inter- 
actions are  probably  most  significant  here:  the  ion-dipole 
forces  between  drug  cations  and  water  molecules  of  the 
solvent,  and  the  ion-ion  interactions  of  the  charged  ring 
nitrogens  with  the  DNA  phosphate  anions.   Ion-dipole  forces 
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between  solvent  molecules  and  the  phosphate  anions  would 
be  less  important  than  the  above  two.   The  two  drug  cations 
are  similar  in  terms  of  the  structural  characteristics  which 
would  affect  the  attractive  forces  between  them  and  water 
molecules.   Both  have  hydrogen-bond  acceptor  amino  groups 
and  protonated  hydrogen-bond  ring  nitrogens.   Thus,  to  a 
first  approximation,  the  total  energy  necessary  to  disrupt 
the  ion-dipole  interactions  between  water  and  3-amino- 
acridinium  cations  would  be  expected  to  be  about  the  same 
as  that  for  water  and  7-aminoquinolinium.   More  important, 
though,  is  a  consideration  of  the  relative  acidities  and 
basicities  of  the  groups.   The  extent  of  hydrogen  bonding 
between  the  solvent  and  a  hydrogen-bond  acceptor  is  related 
to  the  acceptor's  basicity  while  the  extent  between  the 
solvent  and  a  hydrogen-bond  donor  can  be  correlated  with 
the  donor's  acidity.   The  protonated  heterocyclic  ring 
nitrogens  are  hydrogen-bond  donors.   Therefore,  7-amino- 
quinoline  (pK  =  6.65),  being  a  stronger  acid  than  3-amino- 
acridine  (pK  =  8.04),  will  form  stronger  hydrogen  bonds 
with  water  making  the  AH°  for  its  binding  reaction  less 
positive  than  that  of  3-AA.   Likewise,  the  basicity  of  the 
substituted  amino  group  of  7-aminoquinoline  is  somewhat 
greater  (pK  =  0)  than  the  amino  group  on  3-aminoacridinium 
(PKa  -  -1.5),  meaning  the  hydrogen-bond  acceptor  of  the  former 
would  interact  more  strongly  with  water  than  the  latter,  thus, 
generating  a  less  negative  AH°  for  reactions  involving  7-AMQ 
as  a  reactant.   The  above  factors  tend  to  make  the  change  in 


113 


enthalpy  for  the  reaction  of  7-aminoquinolinium  with  DNA 
less  negative  than  that  of  3-aminoacridinium  with  DNA. 
It  is  clear  that  simple  arguments  based  solely  on  ion- 
dipole  energies  of  the  reactant  cations  are  not  sufficient, 
in  themselves,  to  explain  the  more  negative  AH   for  the 
3-AA  system. 

The  amount  of  energy  evolved  upon  formation  of  the 
electrostatic  ion-ion  bonds  between  the  cationic  drug 
molecules  and  the  anionic  phosphates  of  the  DNA  may  explain 
the  slightly  more  negative  AH   for  7-aminoquinolinium 
binding  over  that  for  3-aminoacridinium.   The  cation  having 
the  greater  charge  density  at  the  heterocyclic  nitrogen 
would  bind  more  strongly.   7-aminoquinolinium,  being  the 
stronger  acid,  would  generate  the  stronger  electrostatic 
interaction  of  the  two,  in  agreement  with  the  experimentally 
observed  values.   The  effect  of  the  third  ring  of  3-amino- 
acridine  as  a  steric  barrier  to  hydrogen  bonding  of  the  ring 
nitrogen  would  serve  to  make  AH°  for  the  reaction  involving 
3-AA  more  negative  than  for  7-AMQ,  but  this  effect  would 
likely  be  small. 

Since  the  binding  sites  are  presumed  to  be  identical 
for  the  two  cations,  the  energies  involved  in  disrupting 
bonds  associated  with  the  solvent  and  the  DNA  binding  sites 
would  be  the  same.   Likewise,  enthalpic  contributions  arising 
from  the  formation  of  the  solvent  cages  surrounding  the  bound 
species  would  be  expected  to  be  similar  if  the  binding  geo- 
metries of  the  two  are  alike. 
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Intercalation  of  7-aminoquinolinium  and  3-amino- 
acridinium  may  occur  by  their  "slipping"  between  DNA 
base  pairs.   Such  a  mechanism  would  allow  for  interactions 
between  the  aromatic  base  pairs  and  the  ring  systems  of 
the  cations  and  thus  the  3-AA,  having  a  larger  aromatic 
surface  area,  might  be  expected  to  interact  more  strongly. 

Changes  in  the  standard  enthalpies  for  the  intercala- 
tive  binding  processes  of  the  two  cations  to  DNA  are  more 
negative  than  for  the  surface  binding,  but  are  qualitatively 
similar  in  that  AH°  for  7-aminoquinolinium  is  slightly  more 
negative  (-15.1  Kcal/mole)  than  for  3-aminoacridinium  (-13.7 
Kcal/mole) .   The  absolute  magnitude  of  AH°  for  7-AMQ  is  1.06 
times  greater  than  AH°  for  3-AA  while  the  absolute  magnitude 
of  AH°  for  7-AMQ  is  1.10  times  greater  than  AH°  for  3-AA. 
These  similar  ratios  show  that  in  terms  of  the  enthalpic 
contributions  to  the  free  energy  of  binding,  intercalation 
of  3-aminoacridinium  and  7-aminoquinolinium  are  the  same. 
That  is,  the  larger,  more  lipophilic  aromatic  portion  of 
3-AA  does  not  impart  any  particular  advantage  to  the  com- 
pound in  terms  of  the  enthalpy  of  the  intercalation  process. 

The  entropy  of  surface  binding  of  3-aminoacridinium  is 
considerably  more  positive  than  for  the  similar  process  with 
7-AMQ.   For  3-AA,  TAS°  at  298°K  equals  2.45  Kcal/mole  and 
for  7-AMQ,  TAS°  at  298°K  is  0.94  Kcal/mole.   In  fact,  it  is 
this  relatively  large  difference  in  entropy  which  obviates 
the  enthalpic  advantage  7-AMQ  has  over  3-AA,  resulting  in 
the  more  negative  free  energy  of  binding  of  3-aminoacridinium. 
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Entropy  change  is  a  measure  of  the  extent  to  which  the 
randomness  of  a  system  is  altered  during  the  course  of  a 
reaction,  being  positive  when  the  final  state  is  less 
ordered  than  the  initial  state.   In  terms  of  the  stoichio- 
metrics of  reactants  and  products,  AS°  for  the  two  would 
be  expected  to  be  the  same.   Obviously,  the  solvent  and 
the  solute's  effect  on  its  structure  are  of  great  importance 
in  the  present  instance.   The  role  of  solvent  structure 
may  be  understood  by  envisaging  the  bulk  solvent  as  having 
a  relatively  random,  fluid  structure.   Dispersed  throughout 
the  solvent  are  the  organic  cations,  3-AA  and  7-AMQ,  which 
both  disrupt  the  bulk  solvent  structure  by  their  presence. 
But  the  tricyclic  3-aminoacridinium,  being  more  lipophilic 
and  larger  in  size,  will  demand  a  greater  extent  of  solvent 
structuring  to  accommodate  it.   We  have  already  considered 
the  effects  that  the  polar  regions  of  the  cations  would  have 
on  the  water  via  hydrogen  bonding.   Therefore,  since  the 
DNA  sites  of  the  two  species  are  presumed  to  be  the  same, 
the  binding  of  3-aminoacridinium  will  result  in  a  more 
positive  entropy  change  as  the  cation  is  "removed"  from 
the  aqueous  solvent.   This  would  be  true  even  if  surface 
binding  effected  only  partial  removal  of  the  cationic 
species  from  the  water. 

Large  error  limits  for  the  values  of  AS   for  the  inter- 
calative  binding  of  the  two  cations  to  DNA  mandate  caution 
in  comparing  them,  though  it  appears  that  AS   for  the  binding 
of  7-aminoacridinium  is  considerably  less  positive  than  for 
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3-aminoacridinium.   This  relationship  is  in  parallel  with 
the  surface  binding  AS°  values  and  may  be  due  to  the  same 
factors  as  discussed  above.   The  entropy  changes  given 
in  Table  36  are  per  mole  of  cation  but,  since  the  surface 
and  intercalative  modes  involve  different  ratios  of  DNA 
phosphate  to  cation,  comparisons  of  the  entropies  of  the 
two  processes  can  be  misleading.   However,  AS   terms  for 
surface  binding  are  significantly  positive  whereas  for 
the  intercalative  mode  they  are,  within  experimental  error, 
nearly  zero  or  less  than  zero.   Intercalative  binding  entropy 
changes  include  a  positive  contribution  arising  from  the 
release  of  water  molecules  as  the  cations  are  removed  from 
the  bulk  solvent.   These  would  likely  be  about  the  same  as 
the  AS°  terms  for  surface  binding  meaning  that  the  entropy 
associated  with  the  insertion  of  a  3-AA  cation  already  on 
the  surface  of  the  polymer  would  be  approximately  -10.4 
(-2.2  -  8.2)  e.u.  or  -3.10  Kcal/mole  at  298°K.   For  7-AMQ 
the  value  would  be  about  -17  e.u.  (-14  -  3)  or  -5.07  Kcal/ 
mole  at  298°K.   These  negen tropic  terms  indicate  that  it  is  the 
highly  negative  enthalpies  of  the  intercalative  reactions 
which  result  in  the  overall  negative  free  energy  of  binding. 
Surface  bound  cations  will  apparently  spontaneously  bind 
intercalatively  due  to  the  small  enthalpic  advantage  afforded 
the  latter  process.   It  appears  that,  while  the  enthalpy  . 
arising  from  intercalation  is  decisive,  it  is  only  slightly 
more  negative  for  3-aminoacridinium  than  for  7-aminoquinolinium. 
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Earlier,  we  discussed  the  binding  reactions  in  terms 
of  ligand-raetal  models.   Thus,  it  is  instructive  to  continue 
the  analogy  to  specifically  include  the  chelate  effect. 
The  well  known  chelate  effect  (88)  refers  to  the  more 
negative  free  energy  of  formation  of  polycoordinate  metal 
complexes  formed  from  polydentate  chelates  over  similar 

complexes  formed  from  monodentate  ligands  of  chelates  having 

2+ 

fewer  coordination  sites.   For  example,  log  34  for  Cu(NH-.)  . 

is  11.9  while  log  $2  for  Cu(H2NCH  CH  NH  )  2+  is  20.0  (89)  in 
which  the  former  is  formed  from  4  monodentate  ammonia  ligands 
and  the  latter  from  2  bidentate  ethylenediamine  chelates. 
The  chelate  effect  for  non-transition  metals  is  almost 
entirely  an  entropy  effect  (90).   For  transition  metals, 
greater  crystal  field  stabilization  of  chelates  over  mono- 
dentate complexes  introduces  an  enthalpy  contribution  (91) . 
The  more  positive  entropy  for  chelate  binding  means,  in 
terms  of  probability,  that  it  is  easier  to  form  an  aggregate 
from  a  few  molecules  than  from  many  molecules.   The  three 
DNA  phosphates  involved  in  surface  binding  could  conceivably 
act  as  three  monodentate  ligands,  a  monodentate  ligand  and  a 
bidentate  chelate,  or  as  a  single  tridentate  chelate.   The 
case  in  which  the  phosphates  act  as  a  single  chelate  gives 

the  most  consistent  calculated  values  of  K   and  is  also  the 

s 

one  which  would  be  expected  based  on  the  chelate  effect.   On 
the  other  hand,  the  4  phosphates  involved  in  intercalative 
binding  could  react  as  4  monodentates;  1  monodentate  and  1 
tridentate;  2  monodentates  and  1  bidentate;  2  bidentates;  or 
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as  1  tetradentate  ligand.   Clearly,  the  chelate  effect 

predicts  that  the  tetradentate  ligand  would  generate  the 

4- 
most  positive  entropy.   However,  the  use  of  a  (PO.).    term 

in  the  equilibrium  association  expression  yields  values  of 

KT  which  are  not  nearly  as  consistent  throughout  the  inter- 

2-  . 

calative  titration  region  as  those  when  2(P04)2    is 

employed.   Thus,  an  antichelate  effect  is  seen  for  the 
intercalative  process  in  which  the  greater  AH   for  the  2- 
bidentate-ligand  mechanism  appears  to  be  sufficient  to 
compensate  for  any  positive  entropy  gain  which  might  arise 
from  a  reaction  involving  a  single  tetradentate  chelate. 

When  discussing  the  forces  which  account  for  the 
behavior  of  polymers  and  biologically  important  systems, 
such  as  lipid  membranes  in  aqueous  media,  the  "hydrophobic 
bond"  is  often  invoked.   The  avoidance  here  of  any  reference 
to  this  term  as  an  aid  to  explaining  the  energetics  of  the 
binding  reactions  has  been  intentional.   In  most  cases, 
"hydrophobic  forces"  are  merely  manifestations  of  well 
established  fundamental  forces  arising  from  energies  of 
adhesion  and  cohesion  and  from  entropic  factors.   It  perhaps 
could  be  argued  that  the  "driving  force"  for  the  intercalative 
binding  of  the  compounds  to  DNA  is  the  formation  of  hydrophobic 
bonds  between  the  aromatic  rings  of  the  cations  and  the 
aromatic  portions  of  the  DNA  bases.   But  as  Hildebrand  (92) 
has  pointed  out,  it  is  rare  where  the  energy  involved  in  the 
interaction  of  two  organic  species  in  an  aqueous  system  is 
any  greater  than  it  is  in  a  corresponding  non-aqueous 
environment.   For  example,  the  alkyl  chains  in  soap  micelles 
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are  not  bonded  together  by  a  phobia  for  water  —  they  are 
just  as  adhesive  in  the  absence  of  water.   In  other  words, 
there  is  not  necessarily  any  phobic  response  of  the 
species  toward  water,  but  rather  the  water  molecules 
would  "prefer"  to  be  hydrogen  bonded  together  rather  than 
separate  to  admit  the  immiscible  molecules. 


APPENDIX  I 
FIGURES 


Figure  1.   The  structures  of  selected  substituted 
acridines  and  related  compounds 
(Structures  I  through  XVI) 


I.       3-AMINOACRIDINIUM    CATION       (3-AA) 


II.       3,6- Dl  AMINO- IO-METHYLACRIDINIUM 
CATION     (ACRIFLAVIN) 


III.     3,6-DIAMINOACRIDINIUM    CATION 
(PROFLAVINE) 


IV.    BENZ  [c]ACRIDINIUM     CATION 
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V.  DIBENZ  |ac]  ACRIDINIUM   CATION 

NHCH(CH3KCH2)3NEt2 
OCH, 

VI.  QUINACRINE     CATION    (ATEBRIN, MEPACRIN) 


*NHCH2CH(CH3)(CH2)2N  Et2 


VII.  CHLOROQUIN     CATION 

NH(Y)NR,R2 


VIII.  GENERAL    STRUCTURE    OF   SOME 

ANTIMALARIALS 
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IX.       9-AMINOACRIDINIUM     CATION 


NHCH2CH(OH)CH2NEt2 
OCH, 


X.        ACRANIL     CATION 


XI.       NEOMONACRIN     CATION 


XII.     9-AMIN0-l,2,3,4-TETRAHYDR0ACRIDINIUM 
CATION 
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N02 


NHR 


© 


a)     R"  (CH2)2NH(CH3)2 


© 


XIII. 


b)     R«  (CH2)2N(CH3)2CH2N(CH3)3 

REPORTER    MOLECULES   FROM    REFERENCE 
(3  1) 


(CH3)2N 


N(CH3)2 


XIV.     3,6-DIMETHYLAMINOACRIDINIUM 
CATION     (ACRIDINE    ORANGE) 
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XV.      ETHIDIUM     BROMIDE     CATION 


H2N 


XVI.    7-AMIN0QUIN0LINIUM    CATION     (7-AMQ) 
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Figure  3.   Absorptiometric  titration  of  3-aminoacridinium 
cation  with  DNA,  cesium  salt,  pH  5.90,  25.0°C. 
Background  electrolyte:  C  H„PO  ,  initial 
concentration,  0.15M 

=  0) 


a 

P/D  =  —  ([P] 

b 

P/D  =2.89 

c 

P/D  =5.97 

d 

P/D  =9.65 

e 

P/D  =17.0 

f 

P/D  =26.7 

g 

P/D  =65.3 

P/D  is  the  ratio  of  total  DNA  phosphate  to  total 
3-aminoacridinium  cation,  D  =  1.35  x  10"'  moles. 
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Figure  4.   Absorptiometry  titration  of  3-aminoacridinium 
cation  with  DNA,  cesium  salt,  pH  5.90,  Zb.U   c. 
Background  electrolyte:  CsH2P04,  initial 
concentration,  0.0025M 


a 

P/D  =  —  ([Pi  : 

=  0) 

b 

P/D  =  0.652 

c 

P/D  =  1.63 

d 

P/D  =  2.94 

e 

P/D  =4.57 

f 

P/D  =  5.87 

g 

P/D  =13.7 

P/D  is  the  ratio  of  total  DNA  phosphate  to  total 
3-aminoacridinium  cation,  D  =  1.33  x  10   mOles. 
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Figure  5.   Absorptiometric  titration  of  3-aminoacridinium 

cation  with  DNA,  magnesium  salt,  pH  5.90,  25.0°C, 
Background  electrolyte:  Mg(02CCH3)2,  initial 
concentration,  0.010N 

a  P/D  =  —  ( [P]  =  0) 

b  P/D  =3.14 

c  P/D  =9.57 

d  P/D  =18.9 

e  P/D  =37.7 

f  P/D  =69.1 

g  P/D  =  104 

P/D  is  the  ratio  of  total  DNA  phosphate  to  total 
3-aminoacridinium  cation,  D  =  1.31  x  10" '  moles. 
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Figure  6.   Absorptiometry  titration  of  3-aminoacndiniun^ 

cation  with  DNA,  magnesium  salt,  pH  5.90,  25.0  C. 
Background  electrolyte:  Mg(02CCH3)2,  initial 
concentration,  0.00063N 

=  0) 


a 
b 

P/D  =  —U 
P/D  =  1.63 

c 

P/D  =3.59 

d 

P/D  =  6.4  9 

e 

P/D  =  13.0 

f 

P/D  =  2  6.0 

g 

P/D  =64.9 

P/D  is  the  ratio  of  total  DNA  phosphate  to  total 
3-aminoacridinium  cation,  D  =  1.27  x  10"  moles. 


.8- 


136 


300 


4"3o" 
WAVELENGTH  (nm) 


Figure  7.  Log ( [BHP]/ [BH] )  vs.  log ( [Pt/3] " tBHp3 )  with  0-0025 
CSH2P04  and  6.3  x  10"4N  Mg(02CCH3)2,  as  supportin 
electrolytes,  pH  5.90,  25.0°C 

A   0.0025M  C  H2P04 

O   6.3  x  10~4N  Mg(02CCH3)2 
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Figure  8.   Log ( [BHP] / [BH] )  vs.  log ( [Pt/3] " tBHp] )  with  0.15M 
CSH2P04,  0.010M  CSH2P04,  and  0.010N  Mg(02CCH3)2 
as  supporting  electrolytes,  pH  5.90,  25.0°C 

A   0.15M  CsH2P04 

•   0.01M  C  H2P04 

O   0.010N  Mg(02CCH3)2 
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Figure  9.  Log ( [BHP] / [BH] )  vs.  log ( [Pt/2l -2 [BHP] )  with  0.00 
CSH2P04  and  6.3  x  10  4N  Mg(02CCH3)2  as  supportin 
electrolytes,  pH  5.90,  25.0°C 

A   0.0025M  C  H2P04 

0   6.3  x  10_4N  Mg(02CCH3) 2 
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Figure  10.   Log ( [BHP] / [BH] )  vs.  log ( [Pt/2] "2 [BHP] )  with  0.15 
CSH2P04,  0.01M  CSH2P04,  and  0.010N  Mg(02CCH3)2 
as  supporting  electrolytes,  pH  5.90,  25.0°C 

A   0.15M  CsH"2P04 

•   0.010M  CsH2P04 

O   0.010N  Mg(02CCH3)2 
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Table  22.   Percent  total  3-aminoacridinium  bound  to  DNA 
in  the  presence  of  various  electrolytes  after 
addition  of  excess  DNAa 


Electrolyte 


0.15M  C  H2P04 

0.10M  C  H2P04 

0.010M  C  H2P04 

0.0025M  C  H„PO. 
s  2   4 


0.10M  LiH  PO 

0.010M  LiH  PO. 

6.3  x  10~4M  LiH2P04 


0.010N  Ba(0  CCH  ) 

6.3  x  10"4N  Ba(0  CCH  ) 


0.010N  Mg(02CCH  ) 

6.3  x  10"4N  Mg(02CCH3)2 


Ionic  strength 
(Dc 

P/Dd 

%  of  total 
drug  bound 

0.14 

65 

90 

0.094 

55 

94 

0.010 

17 

97 

0.0028 

14 

97 

0.094 

57 

89 

0.010 

15 

99 

0.0012 

11 

98 

0.017 

103 

77 

2    0.0027 

47 

90 

0.017 

104 

70 

0.0034 

65 

90 

apH  5.90,  T  =  25.0°C. 


Concentrations  before  any  DNA  added. 

c  • 

In  equilvalents  per  liter  for  the  solution  containing  the 

volume  of  DNA  necessary  for  the  P/D  ratio  listed. 
Total  moles  of  DNA  phosphate  per  total  moles  of  3-AA. 
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Table  23.  Initial  and  final  slopes  of  plots  of  log([BHP]/ 
[BH])  vs.  log ( [Pt/m] -q[BHP] )  in  various  concen- 
trations of  CsH  PO   and  Mg(02CCH3)2 


Electrolyte  Initial  slope   Final  slope 


0.15M  C  H„PO. 
s  2   4 

0.010M  C  H-PO. 
s  2   4 

0.0025M  C  H_PO. 
s  2   4 


0.010N  Mg(02CCH  )2 

6.3  x  10~4N  Mg(02CCH3)2 


£ 

IIC 

1^ 

IIC 

1.14 

1.16 

1.14 

1.16 

1.07 

1.23 

1.68 

1.88 

1.10 

1.14 

1.71 

2.16 

0.96 

0.97 

0.96 

0.97 

0.98 

1.0 

0.98 

1.0 

titrations   done  at   25.0°C,    pH   5.90. 
blog( [BHP]/[BH])    vs.    log( [Pt/21 -2 [BHP] ) 
Clog([BHP]/[BH])    vs.    log([Pt/3]-[BHP]) . 
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Table  34.   Equilibrium  association  constants  for  the  binding 
of  3-aminoacridinium  to  DNA,  Ks ,  and  for  the 
binding  of  alkali  and  alkaline  earth  metal  cations 


to  DNA, 

K  a 

M 

Na+ 

Alkali  metals 

Cs+ 

metal 

K+ 

Rb+ 

(CH3)4N+ 

metal 

0.10- 

0.050- 

0.050- 

0.050- 

0.10- 

cone. (M) 

0.025 

0.010 

0.010 

0.025 

0.025 

log  Kg 

5.158- 
5.264 

5.271- 
5.393 

5.272- 
5.420 

5.198- 
5.420 

5.410- 
5.407 

K 


M 


4.06 


8.82 


11.3 


14.6 


-0.09 


log  K 


5.31 


5.43 


5.47 


5.44 


5.41 


metal 


Alkaline  earth  metals 


Mg 


++ 


Ca 


++ 


Sr 


++ 


Ba 


++ 


metal 
cone. (N) 


0.010- 
0.0025 


0.0050- 
0.0013 


0.010- 
0.0025 


0.010- 
0.0025 


log  K 


4.409- 
4.608 


4.501- 
4.594 


4.410- 
4.631 


4.532- 
4.685 


K 


M 


96.1 


70.5 


114 


65.5 


log  K 


4.70 


4.63 


4.74 


4.75 


Calculated  using  equation  (3-5).   All  reactions  at  25.0  + 
0.1°C,  pH  5.9.   For  DNA  and  3-AA  reactant  concentrations 
and  other  experimental  parameters,  refer  to  the  appropriate 
tables  for  the  specific  metals. 
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Table  34  (Continued) 


Counteranion,  H-PO. 
cCounteranion,  CH  COO~ 
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APPENDIX  III 
COMPUTER  PROGRAMS 


These  programs  were  designed  for  the  calculations 
of  apparent  association  constants  and  related  data 
employing  a  Litton-Monroe  model  326  Scientist  computer. 
Terms  not  defined  here  are  either  defined  in  the  text 
or  are  explained  in  the  operator ' s  manual  for  the 
computer.   Brackets  signify  molar  concentrations,  un- 
bracketed  terms  indicate  moles.   V  denotes  volume, 
subscripts  o,  t,  and  °°  refer  to  values  before  any  DNA 
has  been  added,  values  at  a  given  point  in  the  titration, 
and  values  after  excess  DNA  has  been  added,  respectively. 

Program  I 

This  program  has  been  designed  to  compute   the 
following:  [BH] ,  [BHP] ,  Kj,  log ( [P  / 2] -2 [BHP] ) ,  log([BHP]/ 
[BH]  )  ,  K  ,  and  log  (  [P  /_.]  -  [BHP]  )  .   For  the  program  as 
written,  the  total  initial  volume  must  be  10.00  ml  with 
a  sample  pathlength  of  4.00  cm.   Parameters  which  remain 
constant  throughout  a  given  titration  are  stored  in 
registers  1,  2,  4,  and  9  (their  contents  are  listed  in  the 
table  below) .   All  data  is  entered  as  written,  using 
scientific  notation,  and  all  results  are  displayed  in 
their  final  form,  that  is,  no  correction  factors  are 
necessary.   To  run  a  program,  enter  total  volume  of  DNA 
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solution  (yl) ,  press  START,  enter  total  absorbance  of 
the  solution,  press  START.   The  program  will  be  executed 
automatically,  pausing  for  3  to  6  seconds  to  display  each 
of  the  above  values  in  the  sequence  in  which  they  are 
listed.   To  initiate  another  computational  cycle,  repeat 
the  above  procedure.   At  the  completion  of  a  computation, 
the  registers  will  have  stored  within  them  the  following 
information  which  may  be  recalled,  if  necessary. 


Register Information 

0  A. 


4 
5 
6 


(A  -A  )  """BH 
o  °°     o 


[BHP] 


BH 
o 

[BH] 


EVa't 


Register 

Information 

7 

V   in  liters 

8 

Total  1  of  DNA 
added  up  to 
point  t 

9 

P^/p  (moles  of 

DNA  phosphate 
per  1) 

CHG 

SGN 

KI 

log([Pt/2]-2[BHP]t 


The  individual  steps  of  the  program,  along  with  their  key- 
strokes and  Munroe  326  code  numbers  are  presented  below. 
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STEP 

KEY 
STROKES 

CODE 

STEP 

KEY 
STROKES 

CODE 

1 

ST 

300 

51 

) 

027 

2 

8 

010 

52 

f 

024 

3 

STOP 

033 

53 

RCL 

310 

4 

ST 

300 

54 

7 

007 

5 

0 

000 

55 

= 

020 

6 

( 

026 

56 

ST 

300 

7 

RCL 

310 

57 

5 

005 

8 

8 

010 

58 

fEXP 

114 

9 

X 

023 

59 

RCL 

310 

10 

EXP 

014 

60 

3 

003 

11 

CHSN 

013 

61 

T 

024 

12 

6 

006 

62 

RCL 

310 

13 

) 

027 

63 

7 

007 

14 

+ 

021 

64 

= 

020 

15 

• 

012 

65 

ST 

300 

16 

0 

000 

66 

3 

003 

17 

1 

001 

67 

fEXP 

114 

18 

= 

020 

68 

RCL 

310 

19 

ST 

300 

69 

3 

003 

20 

7 

007 

70 

! 

024 

21 

RCL 

310 

71 

( 

026 

22 

9 

Oil 

72 

( 

026 

23 

X 

023 

73 

RCL 

310 

24 

RCL 

310 

74 

6 

006 

25 

8 

010 

75 

- 

022 

26 

T 

024 

76 

( 

026 

27 

RCL 

310 

77 

RCL 

310 

28 

7 

007 

78 

3 

003 

29 

s 

020 

79 

X 

023 

30 

ST 

300 

80 

2 

002 

31 

6 

006 

81 

) 

027 

32 

( 

026 

82 

) 

027 

33 

RCL 

310 

83 

Ax 

025 

34 

1 

001 

84 

2 

002 

35 

- 

022 

85 

X 

023 

36 

RCL 

310 

86 

RCL 

310 

37 

0 

000 

87 

5 

005 

38 

) 

027 

88 

) 

027 

39 

X 

023 

39 

=  • 

020 

40 

RCL 

310 

90 

ST 

300 

41 

2 

002 

91 

CHSN 

013 

42 

= 

020 

92 

fEXP 

114 

43 

ST 

300 

93 

( 

026 

44 

3 

003 

94 

RCL 

310 

45 

( 

026 

95 

6 

006 

46 

RCL 

310 

96 

- 

022 

47 

4 

004 

97 

( 

026 

48 

- 

022 

93 

RCL 

310 

49 

RCL 

310 

99 

3 

003 

50 

3 

003 

100 

X 

023 
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STEP 

KEY 
STROKES 

CODE 

STEP 

KEY 
STROKES 

CODE 

101 

2 

002 

129 

3 

003 

102 

) 

027 

130 

r 

024 

103 

) 

027 

131 

( 

026 

104 

ST 

300 

132 

RCL 

310 

105 

» 

012 

133 

5 

005 

106 

LOG 

061 

134 

X 

023 

107 

fEXP 

114 

135 

( 

026 

108 

fEXP 

114 

136 

RCL 

310 

109 

( 

026 

137 

6 

006 

110 

RCL 

310 

138 

- 

022 

111 

3 

003 

139 

RCL 

310 

112 

f 

024 

140 

003 

003 

113 

RCL 

310 

141 

) 

027 

114 

5 

005 

142 

) 

027 

115 

) 

027 

143 

)   . 

027 

116 

LOG 

061 

144 

fEXP 

114 

117 

fEXP 

114 

145 

fEXP 

114 

118 

RCL 

310 

146 

( 

026 

119 

6 

006 

147 

RCL 

310 

120 

T 

024 

148 

6 

006 

121 

1 

001 

149 

- 

022 

122 

• 

012 

150 

RCL 

310 

123 

5 

005 

151 

3 

003 

124 

= 

020 

152 

) 

027 

125 

ST 

300 

153 

LOG 

061 

126 

6 

006 

154 

STOP 

033 

127 

( 

026 

155 

JUMP 

350 

128 

RCL 

310 

156 
Program  2 

START 

033 

This  commercially  available  program  (Litton-Monroe, 
Orange,  N.J.,  number  9042W)  solves  the  simultaneous 
equations  (2-8),  (2-9),  and  (2-10)  for  the  concentrations 
of  7-aminoquinoline-containing  species,  [AHP] ,  [AH],  and 
[A].   Program  execution  is  outlined  in  the  operator's 
manual.   The  following  table  correlates  the  terminology 
of  the  program  with  the  quantities  used  for  this  specific 
problem  based  on  the  general  equations 
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A1X  +  B^Y  +  C,Z  -   P 


A„X  +  B0Y  +  C0Z  =  Q 
2      2      3 


A3X  +  B3Y  +  C  Z  ■  R 


Program  Variables  or 

designation  Coefficient 

X  [AHP] 

Y  [AH] 

Z  [A] 

Ax  ?AHp  (5332  M-1™"1) 

B1  eAH  (7940  M-1cm_1) 

C1  eA  (347  M~1cm~1) 

P  At/4 

A2  1 

B2  1 

c2  1 

Q  ct  =  co/vt 

A3  0 

B  Ka  (2.239  x  10~7) 


3 
:3 


C_  -[H+]  (-1.259  x  10~6) 


Program  3 

This  program  calculates  association  constants  for  the 
binding  of  7-aminoquinolinium  to  DNA.   The  following  are 
computed  and  displayed  sequentially,  as  written:  K  ,  log([AHP]/ 
[AH]),  log([Pt/2]-2[AHP]),  Ks ,  and  log ( [Pt/3] - [AHP] ) .   A  total 
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initial  volume  of  8.03  ml  and  a  solution  pathlength  of 
4.00  cm  is  required.   Concentrations  of  AHP  and  AH  are 
obtained  from  equations  (2-8),  (2-9),  and  (2-10),  using 
Program  2.   Before  running  a  program  (moles  of  DNA 
phosphate  per  yl) • 1/2  is  loaded  into  register  9.   To 

execute  a  program,  enter  yl  of  DNA  solution,  press  START, 

5  R 

enter  [AHP]  x  10  ,  press  START,  enter  [AH]  x  10  ,  press 

START.   The  concentrations  of  AHP  and  AH  are  converted  to 

the  proper  orders  of  magnitude  in  the  program.   At  the 

completion  of  a  computation  the  registers  will  contain 

the  following  information  for  retrieval,  if  necessary. 


Register Information 

0  [AHP],. 


1 
2 

3 

4 


[AH]t 


log([AHP]t/[AH]t) 


Register 

Information 

7 

V   in  liters 

8 

yl  of  DNA 

9 

Pfc/2   (1/2  moles 

DNA  phosphate 
per  yl) 

CHG 
SGN 


K 


K. 


"Wt 


The  individual  steps  of  the  program,  along  with  their 
keystrokes  and  Munroe  326  code  numbers  are  presented  below. 


STEP 

KEY 

STROKES 

CODE 

STEP 

KEY 
STROKES 

CODE 

1 

ST 

300 

51 

CHSN 

013 

2 

8 

010 

52 

5 

005 

3 

STOP 

033 

53 

= 

020 

4 

ST 

300 

54 

ST 

300 

5 

0 

000 

55 

1 

001 

6 

STOP 

033 

56 

RCL 

310 

7 

ST 

300 

57 

0 

000 

8 

1 

001 

58 

■=■ 

024 

9 

( 

026 

59 

( 

026 

10 

RCL 

310 

60 

( 

026 

11 

8 

010 

61 

RCL 

310 

12 

X 

023 

62 

6 

006 

13 

EXP 

014 

63 

- 

022 

14 

CHSN 

013 

64 

( 

026 

15 

6 

006 

65 

RCL 

310 

16 

) 

027 

66 

0 

000 

17 

+ 

021 

67 

X 

023 

18 

• 

012 

68 

2 

002 

19 

0 

000 

69 

) 

027 

20 

0 

000 

70 

) 

027 

21 

8 

010 

71 

Ax 

025 

22 

= 

020 

72 

2 

002 

23 

ST 

300 

73 

X 

023 

24 

7 

007 

74 

RCL 

310 

25 

( 

026 

75 

1 

001 

26 

RCL 

310 

76 

) 

027 

27 

9 

Oil 

77 

= 

020 

28 

X 

023 

78 

ST 

300 

29 

RCL 

310 

79 

CHSN 

013 

30 

8 

010 

80 

fEXP 

114 

31 

) 

027 

81 

( 

026 

32 

■f 

024 

82 

RCL 

310 

33 

RCL 

310 

83 

0 

000 

34 

7 

007 

84 

f 

024 

35 

= 

020 

85 

RCL 

310 

36 

ST 

300 

86 

1 

001 

37 

6 

006 

87 

) 

027 

38 

RCL 

310 

88 

LOG 

061 

39 

0 

000 

89 

ST 

300 

40 

X 

023 

90 

4 

004 

41 

EXP 

014 

91 

fEXP 

114 

42 

CHSN 

013 

92 

( 

026 

43 

5 

005 

93 

RCL 

310 

44 

= 

020 

94 

6 

006 

45 

ST 

300 

95 

- 

022 

46 

0 

000 

96 

( 

026 

47 

RCL 

310 

97 

RCL 

310 

48 

1 

001 

98 

0 

000 

49 

X 

023 

99 

X 

023 

50 

EXP 

014 

100 

2 

002 

288 


289 


STEP 

KEY 
STROKES 

CODE 

STEP 

KEY 
STROKES 

CODE 

101 

) 

027 

124 

RCL 

310 

102 

) 

027 

125 

0 

000 

103 

LOG 

061 

126 

) 

027 

104 

STOP 

033 

127 

X 

023 

105 

RCL 

310 

128 

RCL 

310 

106 

6 

006 

129 

1 

001 

107 

X 

023 

130 

) 

027 

108 

a 

012 

131 

= 

020 

109 

6 

006 

132 

ST 

300 

110 

6 

006 

133 

# 

012 

111 

6 

006 

134 

fEXP 

114 

112 

6 

006 

135 

( 

026 

113 

= 

020 

136 

RCL 

310 

114 

ST 

300 

137 

2 

002 

115 

2 

002 

138 

- 

022 

116 

RCL 

310 

139 

RCL 

310 

117 

0 

000 

140 

0 

000 

118 

% 

024 

141 

) 

027 

119 

( 

026 

142 

LOG 

061 

120 

( 

026 

143 

STOP 

033 

121 

RCL 

310 

144 

JUMP 

350 

122 

2 

022 

145 

START 

033 

123 

022 

Program  4 

This  program  calculates  the  activity  coefficients  of 
species  at  a  given  ionic  strength  based  on  equation  (1-30). 
The  coefficients  may  be  stored  for  subsequent  determination 
of  K   from  K  [equation  (2-17)]  .  The  program  is  designed  for 
aqueous  systems  at  25.0°C  for  which  B  is  equal  to  3.3  x  10  . 
To  execute  a  program,  enter  I,  press  ST  7,  enter  d  for  the 
i    species,  press  ST  8,  enter  Z.,  press  ST  6,  press  JUMP, 
CHG  SGN,  START.   The  value  of  a.  will  be  displayed.   Store 
the  values  of  a ____.  a   ,  and  an   or  the  respective  values 

rSnF    rsn         F 


for  AHP,  AH,  and  A  in  registers  4,  2,  and  5,  respectively. 
Upon  storing  these  values,  the  correction  term  for 


290 


intercalative  binding  [n  =  2  of  equation  (2-17)  may  be 
computed  by  pressing  JUMP,  9,  START],.  For  the  surface 
binding  correction  factor  press  JUMP,  .,  START. 

At  the  completion  of  a  series  of  computations,  the 
registers  will  contain  the  following  information  which 
may  be  retrieved,  if  necessary. 


Register 


Information 


0 
1 

2 
3 

4 


W/^BH-0^ 
aBHP/(aBH'aPt2) 


a 


BH 


Register Information 

7  I 


8 


CHG 
SGN 


d.  (the  last  one 
entered) 


BHP 

6        Z .  (the  last  one 
i 

entered) 

The  individual  steps  of  the  program,  along  with  their 
keystrokes  and  Munroe  32  6  code  numbers  are  presented  below. 


STEP 

KEY 

STROKES 

CODE 

STEP 

KEY 
STROKES 

CODE 

1 

f .LABEL 

213 

11 

7 

007 

2 

( 

026 

12 

Ax 

025 

3 

( 

026 

13 

• 

012 

4 

RCL 

310 

14 

5 

005 

5 

6 

006 

15 

X 

023 

6 

AX 

025 

16 

, 

012 

7 

2 

002 

17 

5 

005 

3 

X 

023 

18 

0 

000 

9 

( 

026 

19 

9 

011 

10 

RCL 

310 

20 

) 

027 
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STEP 

KEY 

CODE 

STEP 

KEY 

CODE 

STROKES 

STROKES 

21 

) 

027 

51 

RCL 

310 

22 

T 

024 

52 

4 

004 

23 

( 

026 

53 

T 

024 

24 

( 

026 

54 

( 

026 

25 

RCL 

310 

55 

RCL 

310 

26 

7 

007 

56 

5 

005 

27 

Ax 

025 

57 

X 

023 

28 

* 

012 

58 

RCL 

310 

29 

5 

005 

59 

2 

002 

30 

X 

023 

60 

) 

027 

31 

RCL 

310 

61 

) 

027 

32 

8 

010 

62 

ST 

300 

33 

X 

023 

63 

0 

000 

34 

3 

003 

64 

STOP 

033 

35 

• 

012 

65 

JUMP 

350 

36 

3 

003 

66 

. 

012 

37 

EXP 

014 

67 

fLABEL9 

211 

38 

7 

007 

68 

( 

026 

39 

) 

027 

69 

RCL 

310 

40 

+ 

021 

70 

0 

000 

41 

1 

001 

71 

T 

024 

42 

) 

027 

72 

RCL 

310 

43 

) 

027 

73 

5 

005 

44 

CHSN 

013 

74 

) 

027 

45 

ALOG 

161 

75 

ST 

300 

46 

STOP 

033 

76 

1 

001 

47 

JUMP 

350 

77 

STOP 

033 

48 

CHSN 

013 

78 

JUMP 

350 

49 

f LABEL. 

212 

79 

9 

011 

50 

( 

026 

Program  5 

This  program  calculates  the  ionic  strength  of  a 
solution  resulting  from  the  addition  of  a  given  volume  of 
DNA  to  10.00  ml  of  solution  having  a  known  initial  ionic 
strength.   The  contribution  of  background  electrolyte  in 
the  DNA  solution  is  included.   An  initial  and  final  volume 
of  DNA  solution  may  be  entered  to  obtain  an  averaged  ionic 
strength  over  a  range.   Calculations  for  solutions  containing 
monovalent  or  divalent  cations  may  be  made  depending  on  the 
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constants  used  for  steps  51  and  67  of  the  program.   Before 
running  the  program,  the  normality  of  the  10.00  ml  solution 
containing  the  BH  or  AH  is  stored  in  register  0  and  1/2 
(moles  of  DNA  phosphate  +  equivalents  of  supporting  elec- 
trolyte in  the  DNA  solution)  per  yl  is  stored  in  register 
2.   To  execute  a  program,  enter  V,  of  DNA  solution  (yl) , 
press  START,  enter  V„  of  DNA  solution  (yl) ,  press  START. 
The  averaged  ionic  strength  will  be  displayed.   If  the 
total  final  ionic  strength  of  a  solution  containing  a 
given,  single,  volume  of  DNA  is  desired,  execute  the 
program  by  entering  that  value  (yl)  and  pressing  START, 
START.   At  the  completion  of  a  calculation  the  registers 
will  contain  the  following  information  which  may  be 
retrieved,  if  necessary. 

Register Information 

0  Normality,  N,  of  the  10.00  ml  solution  before 
any  DNA  has  been  added. 

1  Average  volume  of  DNA  added  (yl)  =  (V^NA  + 
V2DNA)/2. 

2  1/2 (moles  of  DNA  phosphate  +  equivalents  of 
supporting  electrolyte  in  the  DNA  solution) . 

3  Total  final  average  volume  of  solution  = 
0.01000  +  average  volume  of  DNA  added,  in 
liters . 

2 

4  Total  final  volume  of  solution  x  10  . 

5  Ionic  strength,  I,  of  the  final  solution  if 
register  0  contains  0. 

6  I  of  the  final  solution  if  register  2  contains 
0  (dilution  effect  only) . 
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Register 


Information 


7 
8 
9 


V-jDNA  (yl)  . 

V2DNA  (yl) . 

I  of  final  solution  (averaged  volume  and 
averaged  ionic  strength) . 


CHG 

SGN 


The  individual  steps  of  the  program,  along  with  their 
keystrokes  and  Munroe  32  6  code  numbers  are  presented  below. 


STEP 

KEY 
STROKES 

CODE 

STEP 

KEY 
STROKES 

CODE 

1 

ST 

300 

31 

2 

002 

2 

7 

007 

32 

) 

027 

3 

STOP 

033 

33 

ST 

300 

4 

ST 

300 

34 

3 

003 

5 

8 

010 

35 

( 

026 

6 

( 

026 

36 

RCL 

310 

7 

( 

026 

37 

3 

003 

8 

RCL 

310 

38 

X 

023 

9 

7 

007 

39 

EXP 

014 

10 

+ 

021 

40 

2 

002 

11 

RCL 

310 

41 

) 

027 

12 

8 

010 

42 

ST 

300 

13 

) 

027 

43 

4 

004 

14 

T 

024 

44 

( 

026 

15 

2 

002 

45 

( 

026 

16 

) 

027 

46 

RCL 

310 

17 

ST 

300 

47 

1 

001 

18 

1 

001 

48 

X 

023 

19 

( 

026 

49 

RCL 

310 

20 

( 

026 

50 

2 

002 

21 

RCL 

310 

51 

X 

023 

22 

1 

001 

52 

3* 

003* 

23 

X 

023 

53 

) 

027 

24 

EXP 

014 

54 

-r 

024 

25 

CHSN 

013 

55 

RCL 

310 

26 

6 

006 

56 

3 

003 

27 

) 

027 

57 

) 

027 

28 

+ 

021 

58 

ST 

300 

29 

EXP 

014 

59 

5 

005 

30 

CHSN 

013 

60 

( 

026 
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STEP 

KEY 
STROKES 

CODE 

61 

( 

026 

62 

RCL 

310 

63 

0 

000 

64 

X 

023 

65 

1 

001 

66 

. 

012 

67 

5** 

005** 

68 

) 

027 

69 

T 

024 

70 

RCL 

310 

71 

4 

004 

72 

) 

027 

73 

ST 

300 

74 

6 

006 

75 

( 

026 

76 

RCL 

310 

77 

5 

005 

78 

+ 

021 

STEP    KEY        CODE 
STROKES 

79  RCL  310 

80  6  006 

81  )  027 

82  ST  300 

83  9  Oil 
8  4  STOP  0  33 
8  5  JUMP  350 
86  START  033 

*  51  3  003  for  M2+ 
51  2  002  for  M+1 


**  67  5  005  for  M24 


67  0  000  for  M 


1+ 
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